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ABSTRACT 

In  the  first  portion  of  this  Technical  Memorandum,  the 
general  approach  to  the  complete  digital  computer  simulation  of 
Anti -Submarine  Warfare  (ASW)  Combat  is  formulated  and  discussed. 
The  discussion  centers  upon  the  three  basic  components  of  combat- 
command  and  control  functions,  machines  and  the  environment-  and 
identifies  the  general  parameters  and  dynamic  variables  required 
to  simulate  these  components.  (While  the  primary  concern  is 
ASWT  combat,  the  general  principles  are  applicable  to  any  type  of 
combat).  The  second  portion  of  this  memorandum  presents  a de- 
tailed description  of  the  simulation  of  marine  combat  between 
two  submarines.  In  this  simulation  one  submarine  is  deployed 
on  a Forward  Area  Patrol  Mission.  The  mission  of  the  second 
submarine  is  to  transit  the  patrol  area  of  the  first  submarine. 
This  digital  computer  simulation  contains  many  features  discussed 
in  the  general  approach  and  demonstrates  the  feasibility  of  con- 
structing a general  ASW  Combat  Simulation  for  a digital  computer. 
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1 . INTRODUCTION 

Technological  progress  in  recent  years  has  made  possible 
the  creation  of  highly  sophisticated  military  systems.  Continuing 
technological  progress  during  the  foreseeable  future  will  make  even 
more  sophisticated  systems  possible.  The  experimental  evaluation 
of  existing  systems  and  their  role  in  combat  is  both  expensive  and 
time  consuming.  The  experimental  evaluation  of  proposed  systems 
and  their  role  in  combat  is  not  possible.  Thus,  accurate  analytical 
techniques  for  the  rapid  evaluation  of  the  effectiveness  of  military 
systems,  and  their  employment,  are  required. 

To  accurately  analyze  military  systems  and  their  role  in 
combat,  it  is  necessary  to  address  the  complete  problem.  When 
modern  complex  military  systems  are  considered,  digital  computer 
simulation  is  the  most  promising  method  of  rapidly  addressing  the 
complete  problem.  In  the  past,  analytical  studies  without  the  aid 
of  a digital  computer  have  been  useful  in  the  analysis  of  effective- 
ness and  the  determination  of  optimum  employment  of  military  systems. 
As  these  systems  have  become  more  complex,  however,  the  results  ob- 
tained by  these  studies  have  become  questionable.  A single  such 
study  can  consider  only  a few  facets,  frequently  only  one,  of  a 
realistic  engagement.  Furthermore,  the  few  aspects  of  the  problem 
that  are  considered  may  not  conform  to  reality  because  of  simp- 
lifying assumptions  required  to  yield  a tractable  problem. 

It  is  now  possible  to  construct  a dynamic,  three- 
dimensional,  digital  computer  simulation  of  combat  which  addresses 
the  complete  problem.  In  this  simulation,  the  dynamic  in-situ 
performance  of  all  relevant  subsystems  and  weapons  possessed  by 
each  participant  can  be  reproduced  with  an  accuracy  commensurate 
with  the  current  state  of  understanding  of  these  elements.  Each 
participant  in  the  simulation  can  be  provided  with  a command  and 
control  function  which  accurately  reproduces  the  performance  of 
that  participant ' s decision  maker. 
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Once  constructed,  such  a model  would  provide  a tool 
for  realistically  assessing  the  impact  on  mission  effectiveness 
of  both  subsystem  and  weapon  performance  and  tactics.  Further- 
more, the  model  can  be  constructed  so  as  to  permit  the  rapid 
assessment,  in  terms  of  mission  effectiveness,  of  the  effects 
of  changes  in  system  performance  parameters,  changes  in  the 
tactics  employed  by  the  participants,  and/or  changes  in  environ- 
mental conditions. 

The  purpose  of  this  technical  memorandum  is  two- fold. 
First,  the  general  problem  of  the  development  and  use  of  a 
digital  computer  simulation  of  ASW  combat  is  discussed.  This 
discussion  is  contained  in  Chapters  2 and  3 where  the  general 
principles  for  combat  simulation  are  set  forth,  the  general 
problems  associated  with  sonar  system  simulation  are  identified, 
and  an  approach  to  the  solution  of  these  problems  is  outlined. 

The  second  purpose  of  this  technical  memorandum  is 
to  present  a detailed  description  of  a digital  computer  program 
prepared  to  simulate  the  engagement  between  two  submarines. 

This  particular  simulation  was  prepared  to  assist  in  evaluating 
the  effectiveness  of  certain  sonar  subsystems  proposed  for 
installation  in  the  Post  World  War  II  diesel-electric  attack 
submarines  and  in  the  pre-SXIPJACK  nuclear  attack  submarines. 
This  computer  simulation  is  important  for  two  reasons.  It 
demonstrates  the  feasibility  of  simulating  a combat  situation, 
and  the  results  of  the  simulation  can  materially  aid  in  the 
evaluation  of  the  effectiveness  of  the  subsystems  in  question. 
Chapter  4 contains  a detailed  discussion  of  the  factors  con- 
sidered in  the  two-submarine  combat  simulation.  Sample  re- 
sults obtained  by  exercising  this  model  are  presented  in 
Chapter  5.  Chapter  6 discusses  problem  areas  encountered  in 
the  current  simulation  and  those  anticipated  to  occur  in  the 
expansion  of  the  current  model  into  a general  ASW  combat  simu- 
lation model. 
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2.  THE  PURPOSE  OF  ASW  COMBAT  SIMULATION 


The  purpose  of  combat  simulation  is  to  provide  a 
tool  for  the  evaluation  of  mission  effectiveness  in  terms  of 
the  relevant  system,  tactical,  and  operational  parameters. 

The  question  of  mission  effectiveness  is  fundamental  in  each 
of  the  following  types  of  analyses: 

(a)  The  analysis  of  real  or  hypothetical  military 
systems  in  given  situations  for  specific  tactical  doctrines; 


(b)  The  analysis  of  real  or  hypothetical  tactical 
doctrines  in  given  situations  for  specific  military  systems; 


1 


(c)  The  analysis  of  both  tactical  doctrine  and 
military  system  requirements  necessary  to  accomplish  particular 
types  of  missions. 

A single  digital  computer  simulation,  properly  con- 
structed, can  provide  the  means  of  performing  any  of  these 
analyses.  A detailed  discussion  of  the  principles  involved 
in  the  construction  of  such  a simulation  is  the  subject  of 
the  next  chapter  am  only  the  basic  features  will  be  given 
here.  Basically,  the  simulation  of  combat  involves  the  simula- 
tion of  the  actions  of  all  participants  involved,  both  friendly 
and  enemy.  The  simulated  action  of  each  participant  is  governed 
bv  the  simulated  performance  of  all  relevant  information  sub- 
systems and  weapons  possessed  by  that  participant,  by  the  simu- 
lated performance  of  the  participant's  vehicle,  and  by  a simu- 
lated command  and  control  function.  The  tactics  employed  by 
each  participant  are  incorporated  into  that  participant's 
command  and  control  function.  The  performance  of  all  sub- 
systems and  weapons  is  simulated  in  terms  of  the  subsystem  and 
weapon  performance  parameters  and  the  dynamic  and  environmental 
variables  affecting  performance.  Finally,  to  be  useful,  the 
model  must  be  amenable  to  rapid  alteration  of  number  of  parti- 
cipants, tactics  employed,  subsystem  and  weapon  performance,  and 
character  of  the  dynamic  and  environmental  variables. 
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To  use  Che  simulation  model  in  any  one  of  the  three 
types  of  analysis,  one  specifies: 

1.  the  mission  to  be  simulated; 

2.  the  participants,  both  friendly  and  enemy; 

3.  the  initial  positions  of  each  participant; 

4.  the  subsystems  and  weapons  possessed  by  each 
participant; 

5.  the  tactics  to  be  employed  by  each  participant, 
and 

6.  the  environment. 

The  model  is  then  exercised  to  simulate  the  course 
of  the  engagement  from  initial  encounter  to  conclusion.  The 
exercise  of  the  simulation  reproduces  the  evolution  of  the 
engagement  as  well  as  the  ultimate  result.  Thus,  the  effects 
of  tactics,  subsystem  and  weapon  performance  parameters,  etc., 
upon  mission  effectiveness  of  all  of  the  participants  can  be 
ascertained  and  analyzed. 

The  simulation  of  one  engagement  would  seldom  provide 
sufficient  information  for  a system  or  tactical  analysis.  One 
would  generally  exercise  the  model  many  times,  holding  constant 
those  parameters  to  be  evaluated  while  varying  the  other  param- 
eters of  the  engagement  over  a reasonable  range  of  values. 

A specific  example  of  the  use  of  the  model  will  now 
be  considered.  A frequent  measure  of  mission  effectiveness  is 
kill  probability.  If  it  were  desired  to  evaluate  the  effective- 
ness of  several  different  sonar  systems,  and  if  kill  probability 
were  the  measure  of  mission  effectiveness,  one  would  proceed 
as  follows: 

(a)  Provide  a typical  vessel  with  the  simulated 
performance  of  one  of  the  candidate  systems; 

(b)  Repeatedly  exercise  the  combat  model  for 
relevant  combinations  of: 
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(1)  missions, 

(2)  number  of  participants, 

(3)  subsystems  of  other  participants, 

(4)  tactics  of  all  participants, 

(5)  geographic  regions  of  encounter, 

(6)  geometries  of  encounter, 

(7)  other  dynamic  variables; 

(c)  Compute  a kill  probability  for  each  case 
considered ; 


(d)  Compute  an  average  kill  probability  by  appro- 
priately weighing  and  averaging  the  results 
obtained  in  each  case; 

(e)  Repeat  st-eps  (b)  , (c),  and  (d)  for  each  of 
the  other  candidate  systems. 


In  this  manner,  the  measure  of  mission  effectiveness 
for  each  candidate  system  would  be  determined. 


The  computer  simulation  in  this  example  provides 
important  by-products.  If  several  sets  of  tactics  were  employed 
in  conjuntion  with  the  systems  to  be  evaluated,  the  results  of 
the  simulation  would  automatically  suggest  the  tactics  most 
appropriate  for  use  with  each  system. 


A second  important  by-product  would  be  contained  in 
the  display  of  the  status  of  the  simulation  during  the  course 
of  the  engagement.  This  status  display  would  permit  the  evalu- 
ation of  the  worth  of  each  component  subsystem  of  the  total 
system  during  all  stages  of  an  engagement.  This  information 
should  be  of  significant  value  in  trade-off  studies  of  various 
components  of  the  total  system. 


CONFIDENTIAL 

This  page  is  unclassified. 


CONFIDENTIAL 


an  6500  TRACOR  LANE.  AUSTIN.  TEXAS  78721 


Other  types  of  analyses,  such  as  the  evaluation  of 
the  variation  of  performance  parameters  of  existing  subsystems 
or  the  evaluation  of  tactics,  would  be  carried  out  in  a manner 
similar  to  the  example  just  described. 

The  results  of  a single  exercise  of  a combat  simula- 
tion model  are  illustrated  in  Fig.  1.  An  engagement  between 
a diesel-electric  submarine  patrolling  a barrier  and  another 
diesel-electric  submarine  attempting  to  transit  the  barrier  is 
illustrated.  The  engagement  was  simulated  using  the  Forward 
Area  Patrol  Mission  Computer  Simulation  Model  to  be  described 
later  in  this  memorandum. 

In  this  example,  the  patrol  submarine  possesses  a 
detection  advantage.  The  patrol  submarine  detects  the  transi- 
tor  while  the  latter  is  snorkeling  and  attempts  to  estimate  the 
transitor's  course  and  speed.  Once  this  course  and  speed  have 
been  estimated,  the  patrol  submarine  estimates  the  position  at 
which  the  transitor  will  next  snorkel,  and  proceeds  to  that 
point  to  intercept. 

The  transitor,  however,  is  following  an  evasive  course. 
As  a result,  the  patrol  submarine  computes  an  intercept  position 
which  is  some  distance  from  the  transitor's  base  course.  Never- 
theless, the  patrol  submarine  re-detects  the  transitor  when  he 
snorkels,  and  launches  a weapon.  At  the  point  of  weapon  launch, 
the  kill  probability  is  computed  without  displaying  continuing 
positions  of  the  boats. 
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3.  A GENERAL  ANTI-SUBMARINE  WARFARE 
COMBAT  SIMULATION  MODEL 

For  the  purposes  of  this  Technical  Memorandum,  the 
term  "digital  computer  simulation"  is  defined  to  be  a digital 
model  of  a real  world  situation,  incremented  through  a simu- 
lated time  environment.  The  central  and  most  difficult  tasks 
in  simulation  are  the  definition  of  the  problem  and  the  reduct- 
ion of  the  real  world  situation  into  a meaningful  and  useful 
set  of  abstractions.  In  the  following  section,  the  problem 
of  combat  simulation  will  be  explicitly  defined.  The  remainder 
of  the  chapter  will  be  devoted  to  reduction  of  the  problem  to 
a set  of  abstractions  which  may  be  used  to  adequately  simulate 
ASW  combat  with  a digital  computer. 

3.1  PROBLEM  DEFINITION 

The  general  problem  of  simulating  combat  is  to  simu- 
late the  course  of  an  engagement  from  initial  encounter  to 
ultimate  conclusion  for  the  arbitrary  confrontation  of  an 
arbitrary  number  of  participants  each  of  whom  may  possess: 

1.  Certain  specified  information  sensing  devices, 

2.  Certain  specified  information  processing  devices, 

3.  Certain  communication  devices, 

4.  Certain  specified  weapons, 

5.  A specified  amount  of  previously  acquired 
information. 

6.  A specified  tactical  doctrine, 

7.  The  ability  to  assess  information  available 
to  him  and  to  make  decisions  based  upon  his 
assessment,  and 

8.  The  ability  to  control  his  own  or  his  vehicle's 
motion,  his  own  subsystems  and  his  own  weapons. 
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In  addition,  the  simulation  model  is  required  to  be 
sufficiently  general  so  that  different  systems,  different  tactics, 
different  environmental  conditions,  and  different  situations  can 
be  incorporated  without  undue  difficulty. 

3.2  BASIC  CONSIDERATIONS  IN  ASW  COMBAT  SIMULATION 

Combat  is,  in  reality,  a very  complex  interaction 
among  participants  each  of  whom  possesses  certain  combinations 
of  the  capabilities  listed  in  the  problem  definition.  To  simu- 
late any  complex  process  occurring  in  the  real  world,  it  is  first 
necessary  to  identify  each  essential  component  of  the  process. 

Once  this  has  been  done,  the  simulation  of  the  process  consists 
of  the  simulation  of  the  behavior  of,  and  the  interaction  be- 
tween, each  essential  component  during  the  course  of  time. 

In  the  most  general  terms,  the  essential  components 
of  combat  are: 

1.  The  command  and  control  function  of  each 
participant , 

2.  The  relevant  machines  involved  in  the  combat,  and 

3.  The  medium  (environment)  in  which  the  combat 
occurs . 

The  command  and  control  function  consists  of  the 
processes  of  assessment  and  decision  making  and  the  act  of 
implementing  these  decisions.  Capabilities  5 through  8 of 
each  participant  as  specified  in  the  problem  definition  are 
the  basic  elements  of  the  command  and  control  function.  The 
relevant  machines  involved  in  the  combat  are  the  various 
participating  vehicles,  information  sensors,  processors,  and 
communications  subsystems  and  weapons.  It  should  be  noted 
that  upon  launch  certain  weapons,  such  as  homing  torpedoes, 
should  be  regarded  as  participants  in  the  combat.  Capabilities 
1 through  4 as  specified  in  the  problem  definition  involve 
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the  relevant  machines.  The  medium  is  an  essential  component 
of  combat  because  it  affects  both  the  performance  of  machines 
and  the  command  and  control  function. 


Each  component  of  combat  interacts  with  every  other 
component,  either  directly  or  indirectly.  Examples  of  direct 
interactions  are  (1)  the  radiation  of  information  by  certain 
machines  (such  as  ships),  the  flow  of  this  information  through 
the  medium,  and  the  reception  of  this  information  by  other 
machines  (such  as  sonar  sensors),  (2)  information  flow  from 
sensing  and  processing  machines  to  command  and  control  functions, 
and  command  and  control's  direction  of  certain  machines  such 
as  own  ship  motion.  The  interaction  between  two  command  and 
control  functions  is  an  example  of  an  indirect  interaction.  The 
alteration  of  one  participant's  course  of  action  by  his  command 
and  control  function  may  alter  the  second  participant's  assess- 
ment of  the  situation.  This  new  assessment  could  require  a 
different  course  of  action  on  the  part  of  the  second  partici- 
pant. All  of  these  interactions  will  be  discussed  in  detail 
later. 

The  type  of  combat  considered  places  more  or  less 
obvious  restrictions  upon  the  locale,  the  participants  and 
their  capabilities.  Thus,  ASW  combat  must  occur  in  some  region 
of  the  ocean.  The  simulation  must  be  capable,  however,  of 
treating  arbitrary  oceanographic  conditions  and  arbitrary 
geometries  of  encounter. 

The  participants  in  ASW  combat  as  considered  in  the 
model  being  discussed  here  are  restricted  to  the  following: 

1.  Submarines, 

2.  Surface  Vessels, 

3.  Airplanes  and  Helicopters, 

4.  Certain  Detached  Sonars, 
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5.  Certain  Countermeasure  Devices, 

6.  Certain  Weapons. 


The  information  sensors  and  processors  are  limited 


to : 


following : 


3.3 


1. 

Communication  Equipment, 

2. 

Optical  Equipment, 

3. 

Acoustic  Devices, 

4. 

Electro-Magnetic  Devices, 

5. 

Fire  Control  and  Tracking 

Information  Processors. 

The 

participant's  weapons  are 

restricted  to  the 

1. 

Hedgehog, 

2. 

SUBROC, 

3. 

ASROC, 

4. 

Torpedoes , 

5. 

Mines , 

6. 

Aerial  Bombs, 

7. 

Artillery. 

THE 

ROLE  OF  STATISTICS  IN  ASW 

COMBAT  SIMULATION 

Before  discussing  in  detail  the  construction  of  a 
general  ASW  combat  simulation  model,  the  role  played  by 
statistics  in  such  a model  will  be  considered.  The  description 
of  the  effects  of  the  medium  and  the  performance  of  certain 
machines  is  in  part  statistical  in  nature.  For  example,  the 
propagation  loss  in  the  ocean  from  a fixed  source  to  a fixed 
receiver  with  stationary  oceanographic  conditions  is  found  to 
vary  about  some  fixed  value.  The  variation  in  propagation  loss 
is  apparently  due  to  variation  in  the  microstructure  of  the 
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ocean.  Range  information  provided  by  ranging  subsystems  is 
found  to  vary  about  some  value  even  when  all  known  relevant 
parameters  are  constant.  The  performance  of  detection  sub- 
systems is  generally  expressed  in  terms  of  detection  proba- 
bility. 

The  accurate  simulation  of  all  of  these  processes 
must  include  the  proper  statistical  description.  The  effects 
of  these  statistics  must  be  accounted  for  in  the  evaluation 
of  the  results  produced  by  the  model.  A method  of  doing  this 
would  be  to  average  out  the  statistics  by  repeatedly  exercising 
the  model. 

Tactics,  however,  should  not  be  treated  on  a statis- 
tical basis  in  the  model.  The  command  and  control  function 
for  each  participant  should  consist  of  a specified  plan  of 
action  for  all  contingencies.  To  be  useful,  the  simulation 
of  each  command  and  control  function  must  be  amenable  to  con- 
venient alteration  so  that  logical  tactical  alternatives  can 
be  considered.  Only  a non-statistical  treatment  will  enable 
the  tactics  to  be  studied  in  the  detail  required  by  analyses 
addressed  in  this  memorandum. 

3.4  BASIC  ORGANIZATION  OF  THE  COMBAT  SIMULATION  MODEL 

Proper  organization  of  the  computer  simulation  is 
required  to  accurately  treat  the  interactions  between  the 
various  essential  components  of  the  combat.  This  organiza- 
tion will  be  discussed  in  this  subsection,  and  a detailed 
discussion  of  the  individual  components  will  comprise  the 
remainder  of  the  chapter. 

Figure  2 is  a schematic  representation  of  most  of 
the  direct  interactions  between  the  essential  components  of 
combat.  The  only  interaction  not  shown  is  a collision  or 
near  collision  of  participants.  The  figure  depicts  one  command 
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and  control  function,  the  components  providing  information  to  it, 
the  components  controlled  by  it,  the  interaction  of  all  these 
components,  and  the  influence  of  the  actions  of  the  command 
and  cont  ol  function  upon  the  sensors  possessed  by  the  other 
command  and  control  functions  involved  in  the  combat. 

The  Command  and  Control  functions  are  the  central 
components  of  combat.  Just  as  the  decision  makers  drive  the 
actual  combat,  so  should  the  simulation  of  their  actions  drive 
the  computer  simulation.  If  a computer  program  is  to  accurately 
simulate  combat,  it  must  be  organized  around  the  command  and 
control  function  of  each  participant. 

The  command  and  control  function  may  be  thought  of  as 
continuous  repetition  of  the  following  steps: 

1.  Assessment, 

2.  Decision,  and 

3.  Action. 

The  first  step  is  the  assessment  of  the  situation.  The  second 
step  is  the  decision  to  act  in  accordance  both  with  the  assess- 
ment of  the  situation  and  with  own  objectives.  The  third  step 
is  the  action  to  implement  the  decision. 

Figure  3 is  a schematic  representation  of  the  cyclic 
repetition  of  the  command  and  control  process  and  its  inter- 
action with  the  subsystems  and  vehicles.  The  dashed  lines  in 
Fig.  3 represent  information  flow  between  participants.  The 
lines  are  dashed  to  indicate  that  one  participant  may  or  may 
not  be  aware  of  the  other  participant's  actions. 

The  concept  that  combat  is  the  orderly  repetition  by 
each  participant  of  certain  steps,  as  illustrated  in  Figure  3, 
is  to  some  extent  an  idealization.  Combat  is  in  fact  the 
interaction  of  a number  of  time -dependent  processes,  occurring 
simultaneously.  Nevertheless,  the  concept  of  combat  as  illus- 
trated can  serve  as  a basis  for  the  organization  of  a digital 
computer  simulation,  provided  that  the  interaction  between  the 
various  time  dependent  processes  is  properly  accounted  for. 
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The  time-dependent  processes  in  combat  are: 

1.  Vehicular  Motion, 

2.  Information  Flow, 

3.  Subsystem  and  Weapon  Performance, 


4.  Variations  of  the  Medium, 

5.  The  Assessment  Process, 

6.  The  Decision  Making  Process,  and 

7.  The  Action  Steps. 


All  of  these  processes  proceed  simultaneously,  and  each  process 
influences  one  or  more  of  the  others.  Proper  treatment  of  the 
simulation  of  the  flow  of  time  and  proper  organization  of  the 
computer  program  are  required  to  adequately  account  for  the 
coupling  between  processes. 


In  digital  computer  simulation,  the  continuous  passage 
of  time  must  be  represented  by  finite  time  steps.  The  validity 
of  the  use  of  finite  time  increments  to  simulate  continuous  or 
discrete,  coupled  processes  is  well  known,  and  will  be  only 
briefly  discussed  here. 

The  finite  time  increment  technique  proceeds  as  follows: 

At  the  end  of  each  time  increment  the  changes  in  status  of  all 
time- dependent  factors  occurring  during  that  time  increment  are 
evaluated.  All  coupling  between  the  various  time- dependent 
factors  can  be  considered,  provided  that  the  character  of  the 
coupling  is  known.  It  is  only  necessary  to  organize  the  flow 
of  the  computer  program  in  such  a way  that  the  time-dependent 
factors  and  couplings  are  computed  in  the  proper  order.  Figure  4 
shows  the  basic  program  flow  for  combat  simulation. 

The  "Action  Steps"  require  special  consideration.  An 
Action  Step  that  alters  the  course  of  action  of  any  participant 
can  alter  the  time  dependence  of  all  other  time-dependent  processes. 
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The  time  of  occurrence  of  an  Action  Step  cannot  be  predicted  from 
initial  conditions.  Thus,  the  Action  Steps  couple  all  of  the 
time- dependent  factors  in  an  a priori  unpredictable  way. 

If  the  finite  time  increment  technique  is  used  in  the 
simulation,  an  Action  Step  can  only  occur  at  the  end  of  each 
time  increment.  Therefore,  the  selection  of  a proper  time  incre- 
ment must  be  governed  by  accuracy  considerations  of  the  individual 
time-dependent  factors  and,  in  addition,  the  magnitude  of  the 
time  increment  must  be  selected  so  that  the  results  of  the  simula- 
tion will  not  be  materially  affected  should  the  decision  point 
have  occurred  at  any  time  during  the  time  step. 

3.5  THE  COMMAND  AND  CONTROL  FUNCTION 

Figure  5 is  a schematic  representation  of  the  Command 
and  Control  Function  for  each  participant  and  its  interfaces 
with  the  subsystems  and  vehicles  which  it  controls.  One  of  the 
factors  contributing  to  the  assessment  step  of  Command  and  Control, 
the  Information  Bank,  has  been  subdivided  in  some  detail  because 
each  of  the  indicated  subdivisions  represents  an  important  aspect 
of  the  computer  simulation. 

The  first  step  in  the  Command  and  Control  Function  is 
the  Assessment  Process.  As  indicated  in  Fig.  5,  three  factors 
enter  into  an  assessment.  The  first  factor  is  the  "Information 
Bank".  The  Information  Bank  is  the  cumulative  total  of  all  infor- 
mation received  from  all  participants  in  the  combat.  It  is 
processed  information.  It  may  be  either  processed  and  displayed 
by  the  subsystems  themselves,  or  it  may  be  processed  by  personnel 
subordinate  to  Command  and  Control.  Two  inputs  are  shown  to  the 
Information  Bank,  "Information  Sensors  and  Processors"  and 
"Information  From  Other  Sources".  As  all  information  received 
from  external  sources  is  received  via  sensors,  the  indicated 
division  is  for  simulation  purposes.  The  Information  Sensors 
and  Processors  division  is  information  received  via  subsystems 
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whose  performance  is  to  be  simulated  in  the  computer  simulation 
(for  example,  information  received  via  sonar).  The  "Information 
From  Other  Sources"  division  is  information  received  via  sub- 
systems whose  performance  normally  is  not  to  be  simulated  (for 
example,  radio  messages).  The  simulat  on  of  the  subsystems 
providing  information  to  the  information  bank  will  be  the  sub- 
ject of  a later  section  of  this  chapter. 

The  second  element  entering  into  the  Assessment  Process 
is  "Knowledge  and  Experience".  This  is  information,  either 
correct  or  erroneous,  possessed  by  the  decision  maker,  but  not 
derived  from  the  other  participants.  For  example,  the  decision 
maker  in  a submarine  will  have  some  idea  of  the  existing  ocean- 
ographic conditions  and  of  the  detection  range  of  his  detection 
subsystems  as  a function  of  these  oceanographic  conditions.  Thus, 
at  the  time  a detection  is  achieved,  his  assessment  would  nec- 
essarily be  dependent  upon  his  estimated  maximum  detection  range. 

The  third  factor  entering  into  the  Assessment  Process 
is  the  "Tactical  Situation".  The  decision  maker's  assessment 
is  based  not  only  on  the  information  available  to  him,  but  also 
upon  what  he  is  and  has  been  trying  to  accomplish.  For  example, 
a submarine's  loss  of  contact  of  a surface  vessel  might  be  inter- 
preted differently  by  the  decision  maker  in  the  submarine  if  he 
were  trying  to  evade  the  surface  ship  or  if  he  were  attempting 
to  close  with  the  surface  ship  and  launch  a weapon. 

The  result  of  the  assessment  process  may  be  thought  of 
as  the  decision  maker's  total  current  evaluation  of  the  situation. 
This  result  is  one  of  the  three  factors  entering  into  the  de- 
cision process.  The  second  factor  entering  into  this  process 
is  the  "Tactical  Doctrine".  This  is  the  mission  of  the  partici- 
pant, the  guide  lines  by  which  he  should  effect  this  mission  and 
the  constraints  under  which  he  should  operate. 
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The  third  factor  entering  into  the  decision  process 
is  the  "Knowledge  and  Experience"  of  the  decision  maker.  This 
controls  the  implementation  of  Tactical  Doctrine  in  accordance 
with  the  decision  maker's  assessment  of  the  situation.  It 
might  also  include  reasonable  predictions  of  future  action  by 
the  other  participants. 

The  objective  of  the  decision  process  is  the  deter- 
mination of  that  course  of  action  most  favorable  to  success 
of  the  decision  maker's  mission.  The  action  step  in  Fig.  5 
denotes  the  act  of  implementing  this  course  of  action  once 
determined.  Since  a possible  result  of  the  decision  process 
is  to  continue  on  present  course  of  action,  the  action  step 
may  consist  of  doing  nothing  to  alter  the  present  course  of 
action. 

There  are  only  five  broad  classes  of  action  available 
to  the  decision  maker  in  combat.  They  are: 

1.  Change  own  ship's  velocity, 

2.  Change  own  ship's  activity, 

3.  Alter  an  Information  Channel, 

4.  Fire  a weapon, 

5.  Create  a new  participant  (e.g.  launch  a torpedo, 
activitate  a countermeasure  device). 

Ship's  velocity  is  used  here  in  the  mathematical  sense  and 
includes  such  things  as  changes  in  course,  speed  and  rate  of 
descent.  A change  in  own  ship's  activity  is  any  change  that 
alters  the  information  received  or  radiated  by  that  participant. 
The  commencement  of  snorkeling  by  a submarine  is  an  example. 

The  last  three  classes  of  action  include  all  controls  that  a 
decision  maker  has  upon  his  own  weapons  and  subsystems. 
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The  action  step  of  each  decision  maker  affects  some 
vehicles  and  subsystems  directly  and  others  indirectly.  The 
effects  of  the  action  step  will  be  discussed  more  fully  in 
the  next  chapter. 

It  is  probably  impossible  to  simulate  the  general 
command  and  control  function  because  of  the  enormous  variety 
of  situations  in  which  a decision  maker  might  find  himself. 

It  is  possible  to  simulate  this  function  for  any  given  type 
of  engagement;  the  method  of  achieving  this  will  now  be  des- 
cribed. It  will  be  clear  that  many  portions  of  the  simula- 
tion for  one  type  of  engagement  will  apply  to  other  types 
of  engagements. 

The  key  to  the  simulation  of  command  and  control  is 
the  tactical  situation  of  each  participant.  A participant's 
tactical  situation  is  his  immediate  goal  and  his  method  of 
accomplishing  the  goal.  Examples  of  tactical  situations  are: 
patrol  according  to  a preassigned  patrol  pattern,  close  with 
a contact  in  order  to  gather  information  about  it,  close  with 
a contact  on  a lead  intercept  course,  evade  a contact  by 
diving  to  maximum  depth  and  running  away  from  the  contact  at 
maximum  speed.  It  should  be  emphasized  that  a tactical 
situation  applies  to  each  participant  and  not  to  the  engage- 
ment as  a whole. 

It  is  necessary  to  consider  each  tactical  situation 
in  which  each  participant  may  find  himself.  It  is  also  nec- 
essary to  ascertain  all  information  that  the  decision  maker 
may  possess  exclusive  of  that  provided  by  his  information 
bank.  (It  is  assumed  at  this  point  that  each  participant's 
information  bank  has  been  simulated.) 

For  each  tactical  situation,  the  following  must  be 
determined : 
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1.  The  participant's  course  of  action; 

2.  The  alteration  in  this  course  of  action  that  may 
be  required  while  in  this  tactical  situation, 
and  the  information  apprising  the  decision  maker 
that  he  should  take  this  course  of  action; 

3.  The  tactical  situations  which  may  logically 
follow ; 

4.  The  information  notifying  the  decision  maker  to 
go  to  another  tactical  situation; 

5.  The  order  of  precedence  of  the  tactical  situations 
that  may  follow.  (This  is  required  in  case 
available  information  suggests  to  the  decision 
maker  that  he  should  go  to  two  or  more  tactical 
situations ) . 

Once  the  tactical  situations  have  been  organized  as  described, 
it  is  a straightf orward  matter  to  write  the  computer  program 
to  simulate  the  command  and  control  function.  An  approach  to 
this  simulation  will  now  be  submitted. 

After  each  participant's  information  bank  has  been 
updated  at  the  end  of  a time  step,  the  program  branches  to  the 
current  tactical  situation  of  participant  number  one.  The 
program  checks  the  information  available  to  that  participant, 
both  from  the  information  bank  and  from  knowledge  and  experience, 
against  the  information  required  to  inform  that  decision  maker 
to  go  to  another  tactical  situation.  This  check  is  done  in 
order  of  precedence  of  possible  tactical  situations.  Should 
the  information  available  indicate  that  the  decision  maker 
should  go  to  another  tactical  situation,  the  necessary  action 
to  effect  this  transfer  is  simulated,  and  a tag  is  set  so  that 
the  program  will  branch  to  the  new  tactical  situation  at  the 
end  of  the  next  time  step. 


I 


24 

UNCLASSIFIED 


UNCLASSIFIED 


LANE. 


AUSTIN  TEXAS 


78721 


Should  the  information  available  to  the  decision  maker 
indicate  that  he  should  not  go  to  another  tactical  situation, 
the  information  is  reexamined  to  determine  if  a modification 
of  course  of  action  is  required  while  in  this  same  tactical 
situation.  If  a modification  is  required  the  appropriate 
action  to  accomplish  this  is  simulated. 

The  above  procedures  are  repeated  for  each  partici- 
pant in  the  engagement  at  the  end  of  each  time  increment. 

It  should  be  possible  to  simulate  the  command  and 
control  function  quite  adequately,  provided  guidance  from 
those  experienced  in  actual  ASW  operations  can  be  obtained. 

3.6  MACHINES  AND  THE  MEDIUM 

The  performance  of  the  machines  of  combat  and  the 
characteristics  of  the  medium  in  which  the  combat  occurs  are 
intimately  related.  Hence,  it  is  convenient  to  combine  the 
discussion  of  these  two  essential  components.  It  is  also 
convenient  to  subdivide  the  machines  into  two  groups  and  dis- 
cuss each  group  separately.  The  first  subdivision  consists 
of  those  machines  that  directly  interface  with  the  assessment 
process  of  command  and  control  and  the  second  consists  of 
those  which  JitcvCly  interface  with  the  action  step  of  command 
and  control. 

Before  discussing  these  machines,  the  relation  be- 
tween subsystem  performance  prediction  and  the  ASW  Combat 
Simulation  Model  should  be  set  forth.  It  is  not  the  function 
of  the  combat  simulation  to  predict  performance  of  subsystems. 
Rather,  a major  purpose  of  the  combat  simulation  is  to  assess 
the  influence  of  subsystem  performance  upon  the  mission  effect- 
iveness of  the  vehicles  upon  which  the  subsystems  are  installed. 
The  combat  simulation  requires  that  the  performance  of  each 
subsystem  be  specified  as  a function,  if  appropriate,  of: 
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1.  Performance  Variables, 

2.  Dynamic  Variables, 

3.  Environmental  Variables,  and 

4.  Time. 


The  performance,  dynamic,  and  environmental  variables  are 
themselves  time-dependent.  Thus,  the  combat  simulation  is 
required  to  treat  these  variables  as  a function  of  time. 

There  is  no  intent  to  minimize  the  importance  of  subsystem 
performance  analysis.  Indeed,  the  realism  of  the  simulation 
model  rests  more  upon  the  adequacy  of  subsystem  performance 
simulation  than  upon  any  other  factor.  The  purpose  here  is 
to  point  out  that  subsystem  performance  analysis  and  combat 
simulation  are  separate  problems.  It  is  anticipated  that 
the  analysis  of  relevant  ASW  subsystems  will  proceed  in 
parallel  with  the  development  of  the  combat  simulation  model. 

3.6.1  Machines  That  Interface  With  the  Assessment  Process 

The  machines  that  interface  with  the  assessment 
process  of  command  and  control  are  the  information  sensors 
and  processors.  These  sensors  and  processors  are  most  logically 
discussed  in  terms  of  information  channels.  An  information 
channel  consists  of  information  sources,  noise  sources,  the 
medium  through  which  the  information  and  noise  is  transmitted 
to  a sensor,  the  sensor,  the  processor,  possibly  additional 
information  required  by  the  processor,  and  the  presentation 
of  the  processed  information  to  command  and  control.  The 
performance  of  human  operators  or  observers  is  considered 
to  be  part  of  the  performance  of  the  processor.  A typical 
information  channel  is  illustrated  in  Fig.  6. 

Each  information  channel  can  be  characterized  by  the 
particular  sensor  used  for  reception  of  the  information  from 
the  medium  and  by  the  processor  which  processes  this  informa- 
tion. Information  sensors  include  communication  equipment. 
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radar  equipment,  optical  sensing  equipment  (pelorus,  periscope, 
etc.),  sonar  equipment,  etc.  In  ASW  combat,  the  primary  sensors 
are  the  sonar  devices.  Thus,  this  section  will  discuss  princi- 
pally the  sonar  information  channels.  Should  it  be  necessary 
to  simulate  other  information  channels,  they  can  be  included 
in  the  simulation,  with  no  significant  modification  of  the 
computer  program  organization.  It  is  noted  that  the  effect 
of  another  channel  can  be  included  without  actually  simulating 
the  channel.  For  example,  the  information  received  by  radio 
can  be  included  without  actually  simulating  the  radio. 

The  simulation  of  an  information  channel  can  be 
divided  into  two  parts.  This  division  enables  the  essential 
factor  to  be  logically  isolated  and  discussed,  and  permits 
the  simulation  of  the  channel  to  be  carried  out  in  a general 
manner.  The  first  part  is  the  simulation  of  the  information 
and  noise  field  at  the  sensor  as  a function  of  the  appropriate 
dynamic  and  environmental  variables  and  time.  The  second  part 
is  the  simulation  of  the  performance  of  the  sensor  and  processor 
as  a function  of  the  information  and  noise  field  at  the  sensor, 
dynamic  variables,  the  performance  variables  of  the  sensor  and 
processor,  time,  and  other  information  required  by  the  processor. 
The  latter  division  will  be  discussed  first. 

3.6. 1.1  Sonar  information  sensors  and  processors.  Sonar 

information  channels  can  be  divided  into  four  basic  categories 
according  to  the  function  performed.  Each  category  is  composed 
of  both  active  and  passive  sonars.  The  categories  are: 

1.  Detection, 

2.  Localization, 

3.  Track  Determination,  and 

4.  Classification, 

and  the  specific  function  performed  by  each  is  shown  in  Table 
I as  follows: 
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TABLE  I.  CATEGORIES  OF  SONAR  SYSTEMS 


Information  Channel 

Function 

Category 

Detection 

Detection 

Contact  Maintenance 

Localization 

Contact's  Range  Determination 
Contact's  Bearing  Determination 

Track  Determination 

Contact's  Course  and 
Speed  Determination 

Classification 

Contact  Classification 

A particular  sonar  can  be  a part  of  more  than  one 
information  channel;  for  example,  a sonar  can  be  a part  of  a 
detection  information  channel  and  a localization  information 
c hanne 1 . 

The  performance  to  be  simulated  as  a function  of  time 
for  each  category  is  shown  in  Table  II. 


TABLE  II.  PERFORMANCE  TO  BE  SIMULATED 
BY  INFORMATION  CHANNELS 


Information  Channel 
Category 

Performance  to  be  Simulated 

Detection 

Probability  of  Detection 
Contact  Maintenance 

Localization 

Contact's  Range  Error 
Contact's  Bearing  Error 

Track  Determination 

Contact's  Course  Error 
Contact's  Speed  Error 

Classification 

Probability  of  Classification 

In  the  simulation,  the  exact  velocity  vector  of  each 
participant,  the  exact  range  between  participants,  and  the  exact 
bearing  of  each  participant  with  respect  to  every  other  partici- 


pant are  available.  Therefore,  the  simulation  of  the  perfor- 
mance of  the  localization  and  tracking  information  channels  is 
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required  to  produce  only  range,  bearing,  course,  and  speed 
errors.  These  errors  are  then  combined  with  exact  values  and 
the  resulting  approximate  values  are  presented  to  command 
and  control. 

The  performance  variables  to  be  considered  include 
such  things  as: 

1.  Operating  Mode, 

2.  Integration  Time, 

3.  Processor  Threshold, 

4.  Operator  Alert  Condition, 

5.  Pinging  Cycles, 

6.  Spoke  Suppression,  and 

7.  Performance  degradation  resulting  from  failure 
or  degradation  of  subsystem  components. 

The  performance  of  the  tracking  information  channels 
is  a function  of  own  ship  motion  information  in  addition  to 
other  variables. 

To  simulate  the  sonar  information  sensor  and  processor 
portions  of  the  information  channels  it  is  necessary  to  know 
the  performance  of  each  channel  as  a function  of  the  informa- 
tion and  noise  field  at  the  sensor,  the  appropriate  performance 
variables,  time  and  other  information  as  required.  Once  this 
is  known,  a computer  program  can  be  constructed  which  simulates 
this  performance.  The  performance  simulation  program  must  be 
linked  to  the  appropriate  information  and  noise  field  simula- 
tions (to  be  described  later)  and  to  the  command  and  control 
function.  The  link  to  command  and  control  must  be  to  the 
assessment  both  from  and  to  the  action  step.  The  link  to  the 
action  step  is  required  so  command  and  control  can  select 
those  performance  parameters  subject  to  its  control. 

The  knowledge  of  the  performance  of  the  sensors  and 
processors  may  be  obtained  from  experimental  measurements  or 
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from  theoretical  performance  predictions.  The  exact  form  of 
expression  of  this  performance  is  immaterial,  so  long  as  it  is 
amenable  to  reduction  to  analytical  formulas  or  numerical  tables. 

If  it  is  desired  to  determine  the  sonar  performance 
required  to  provide  a specified  capability  to  one  of  the  partici- 
pants, completely  hypothetical  performance  may  be  inserted  in 
the  simulation  program. 

The  simulation  of  sonar  performance  must  be  considered 
on  an  individual  basis.  However,  the  computer  simulation  can  be 
constructed  so  that  the  simulated  performance  of  one  sonar  can 
be  replaced  by  the  simulated  performance  of  another  sonar  in  the 
same  category  with  relative  ease.  This  can  be  done  because  all 
sonars  in  the  same  category  furnish  the  same  information  to  command 
and  control.  The  only  modifications  necessary  are  the  replace- 
ment of  one  performance  simulation  by  the  other,  the  insertion 
of  the  appropriate  information  and  noise  field  calculations  for 
the  new  sonar,  and  possibly  certain  minor  modifications  of  command 
and  control  should  the  performance  variables  of  the  two  sonars 
differ. 

3. 6. 1.2  The  information  and  noise  field  at  the  sonar  sensor. 

The  information  and  background  noise  received  by  a sonar 
sensor  are  contained  in  the  variations  of  intensity  of  the 
acoustic  field  in  the  immediate  vicinity  of  the  sensor.  In 
principle,  this  variation  of  intensity  can  be  found  by  construct- 
ing the  appropriate  partial  differential  equation  and  then  inte- 
grating this  differential  equation  subject  to  appropriate  boundary 
conditions.  In  practice,  this  procedure  cannot  be  carried  out 
exactly,  but  it  can  be  carried  out  to  various  degrees  of  approxi- 
mation. The  precise  degree  of  approximation  required  by  the  com- 
bat simulation  model  is  dictated  by  the  performance  simulation 
of  the  sensor  and  processor  of  each  information  channel.  The  in- 
formation and  noise  field  at  each  sensor  must  be  specified  in 
sufficient  detail  to  permit  adequate  simulation  of  the  processed 
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information  delivered  by  the  information  channels  associated  with 
that  sensor.  In  the  remainder  of  this  subsection,  the  basic  factors 
which  must  be  considered  in  the  simulation  will  be  listed.  In  the 
next  subsection  a method  of  simulating  the  information  and  noise 
field  to  first  order  will  be  presented. 

To  simulate  the  information  and  noise  fields  at  sonar 
sensors  in  a dynamic  computer  model,  it  is  necessary  to  consider: 

1.  The  sources  of  information, 

2.  The  sources  of  background  noise, 

3.  The  influence  of  the  medium  in  transmitting  the 
information  and  noise  to  the  sensor. 

Each  of  these  factors  must  be  considered  as  a function  of  the 
dynamic  and  environmental  variables  of  the  combat. 

Insofar  as  passive  sonar  is  concerned,  most  vehicles 
in  ASW  combat  serve  as  information  sources  in  that  they  radiate 
acoustical  energy  into  the  medium.  The  intensity  of  this  acousti- 
cal energy  at  unit  distance  is  a function  of: 

1.  The  vehicle, 

2.  The  vehicle's  speed,  and  depth, 

3.  Other  noise  making  activities  on  board  the  vehicle, 

4.  Direction  in  which  the  sound  is  radiated, 

5.  Time,  and 

6.  Frequency. 

The  "other  noise  making  activities"  factor  includes  such  things 
as  active  sonar  transmission  and  acoustic  countermeasures  as  well 
as  other  noise  making  activities  aboard  the  vehicle. 
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Insofar  as  active  sonar  is  concerned,  the  vehicles  in 
ASW  combat  may  serve  as  information  sources  in  that  they  may  re- 
flect energy  which  originates  as  active  sonar  pulses.  The  target 
reflectivity  is  referred  to  as  target  strength  and  is  a function 
of : 


1. 

Target  vehicle, 

2. 

Pulse  length, 

3. 

Frequency , 

4. 

Target  aspect, 

5. 

Range , 

6. 

Depth,  and 

7. 

Active  sonar  processor 

resolution  time. 

In 

addition  to  information, 

the  sonar  sensors 

receive  extraneous  signals.  These  extraneous  signals  are  called 
sonar  background  noise.  Fig.  7 shows  the  major  sources  of  sonar 
background  noise. 

The  first  type  of  background  noise  concerns  only  an 
active  sonar  signal  and  is  termed  reverberation.  Reverberation 
is  the  transmitted  energy  which  is  reflected  back  to  a receiver 
by  the  medium  and  its  boundaries.  Reverberation  is  a function  of 


1. 

Source  level, 

2. 

Transmit  beam  pattern. 

3. 

Ping  length, 

4. 

Time  measured  from  end  of 

Ping, 

5. 

The  propagation  path, 

6. 

Oceanographic  parameters. 

7. 

Scattering  cross  sections 

for  reverberation,  and 

8. 

Frequency. 
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The  second  type  of  extraneous  signal  is  termed  under- 
water ambient  noise.  The  underwater  ambient  noise  consists  of 
the  ever  present  sea  noise  plus,  possibly,  noise  from  other 
sources  such  as  biological  noise,  rain  noise,  and  so  forth.  Thus, 
the  underwater  ambient  noise  is  a function  of: 

1.  The  source  levels  qf  the  contributor, 

2.  The  path  from  contributor  to  sensor, 

3.  Oceanographic  parameters, 

4.  Time,  and 

5.  Frequency. 

The  third  type  of  extraneous  information  is  termed 
sonar  self-noise  and  consists  of  platform  noise  and  equipment 
noise.  Sonar  self-noise  is  a function  of: 

1.  The  source 'level  of  the  noise  contributor  which 
is  in  turn  a function  of  the  vehicle,  vehicle's 
speed,  vehicle's  position,  and  any  other  noise 
making  activity, 

2.  The  path  from  contributor  to  sensor, 

3 . Time , and 

4.  Frequency. 

The  ocean  medium  exerts  an  influence  on  each  sonar 
information  channel.  The  ocean  serves  as  the  means  by  which 
information  and  noise  is  physically  propagated  to  the  sensor. 

The  intensities  of  the  information  and  noise  decrease,  in 
general,  as  they  propagate  from  the  source.  This  decrease  in 
intensity  is  termed  propagation  loss,  and  is  the  major  medium 
influence  which  must  be  accounted  for  in  the  combat  simulation. 

It  may  also  be  necessary  to  account  for  alterations  in  wave  form 
as  the  energy  propagates.  Propagation  loss  is  a function  of: 
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1.  The  position  of  the  source, 

2.  The  position  of  the  receiver, 

3.  The  paths  over  which  the  sound  travels, 

4.  Oceanographic  parameters,  and 

5.  Frequency. 

Thus,  to  simulate  the  information  and  noise  field 
at  each  sensor  it  is  necessary  to  know,  as  a function  of  the 
listed  time -dependent  parameters,  the  radiated  noise  level,  the 
target  strength,  and,  the  self-noise  level  at  each  participant. 

In  addition  a method  of  calculating  the  reverberation  level 
for  each  active  sonar,  the  propagation  loss  between  each  pair 
of  participants,  and  the  ambient  noise  levels  must  be  included. 
Finally,  coupling,  if  any,  between  these  factors  must  be  included. 

As  in  the  case. of  sonar  performance,  the  knowledge  of 
the  factors  involved  in  the  information  and  noise  field  calcula- 
tions may  be  obtained  from  experimental  measurements  or  theore- 
tical predictions.  At  present,  experimental  data  is  the  best 
source  for  radiated  noise  levels,  self-noise  levels,  ambient 
noise  levels,  and  target  strengths.  Theoretical  calculations 
are  the  best  source  of  reverberation  and  propagation  loss  due 
to  the  lack  of  experimental  data  for  the  wide  range  of  ocean- 
ographic parameters  necessary  to  generalize  the  simulation. 

3 . 6 . 1 . 3 A first  order  approximation  to  the  information  and 
noise  field  at  the  sonar  sensor.  In  many  cases  sonar  performance 
can  be  expressed  in  terms  of  an  information  and  noise  field  which 
is  specified  by  the  source's  range  and  bearing  and  the  short- 
time  average  signal-to-noise  level  at  the  output  of  the  sonar 
beamformer.  This  representation  in  effect  shifts  the  informa- 
tion and  noise  field-sonar  interface  from  the  sensor  to  the 
output  of  the  beamformer. 


36 


UNCLASSIFIED 


UNCLASSIFIED 


6500 


TRACOR 


LANE 


AUSTIN,  TEXAS  78721 


In  the  model  to  be  described  in  the  next  chapter, 
sonar  performance  is  simulated  in  terms  of  short-time  average 
signal-to-noise  levels.  These  signal-to-noise  levels  are 
computed  utilizing  equations  (1)  and  (2),  below. 

Active  Sonar 

S/N  = [S0  - PL  - H?  - HR  + TS]  -‘[10  log  (Na  + Ng  + R)  1 (1) 

where 


SQ  = the  source  level  on  the  axis  of  the  main  lobe  (dB) ; 

PL  = the  two-way  propagation  loss  (dB) ; 

H,p  = transmission  pattern  losses  to  account  for  the  target 
not  lying  on  the  axis  of  the  main  transmit  lobe  (dB) ; 

HR  = receiving  pattern  losses  to  account  for  the  target  not 
lying  on  the  axis  of  the  main  receive  lobe  (dB) ; 

TS  = target  strength  (dB) ; 

= intensity  of  ambient  noise  integrated  over  the  receive 
beam  pattern; 

Ng  = intensity  of  self-noise  integrated  over  the  receive 
beam  pattern; 

R = intensity  of  reverberation  integrated  over  the  receive 
beam  pattern. 

Passive  Sonar 

S/N  = [SR  - PL  - Hr]  - [10  log  (Na  + Ng) 1 (2) 

where 

SR  = the  radiated  signal  of  the  target  integrated 
over  the  receiving  beam  pattern  (dB) ; 
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PL  = the  one-way  propagation  loss  (dB) ; 

= receiving  pattern  losses  to  account  for  the 

target  not  lying  on  the  axis  of  the  main  receive 
lobe  (dB) ; 

= intensity  of  ambient  noise  integrated  over  the 
receive  beam  pattern; 

Ng  = intensity  of  self -noise  integrated  over  the 
receive  beam  pattern; 

It  should  be  noted  that  the  above  equations  are  general  provided 
that  the  indicated  separation  of  effects  is  valid,  and  that  each 
term  in  the  equation  can  require  a very  sophisticated  analysis 
before  it  can  be  evaluated.  For  example,  the  evaluation  of  the 
propagation  loss  (PL)  term  can  be  obtained  in  many  different  ways 
ranging  from  a hand  calculation  to  a complex  simulation  including 
all  oceanographic  factors ' affecting  propagation  loss. 

3.6.2  Machines  That  Interface  with  the  Action  Step 

Each  participant  in  an  ASW  Combat  Simulation  has  five 
possible  tactical  decisions  available  to  him,  as  has  been  dis- 
cussed previously.  These  decisions  concern  change  of  the  parti- 
cipant's velocity  vector  (i.e.,  change  course,  change  speed), 
change  in  the  participant's  activity  (i.e.,  start  or  stop  snorkel 
ing) , weapon  launch,  and  control  of  information  channels.  An 
approach  to  the  simulation  of  the  factors  directly  affected  by 
this  implementation  will  now  be  presented. 

The  motion  of  each  participant  in  the  simulation  is 
described  by  a differential  equation,  which  can  be  integrated 
at  each  time  step  to  produce  the  new  position  and  orientation 
of  the  participant  at  that  time. 

An  activity  change  is  anything  done  by  a participant 
other  than  change  of  velocity  which  alters  the  nature  of  the 
information  either  received  or  radiated,  or  both.  For  example, 
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if  a submarine  starts  snorkeling,  the  radiated  noise  as  well  as 
the  self-noise  of  that  participant  are  markedly  increased. 

To  simulate  the  change  of  activity,  radiated  noise  and 
self-noise  information  must  be  provided  as  a function  of  the 
possible  activity  for  each  participant. 

Each  participant  in  the  simulation  may  have  the  capa- 
bility of  launching  a weapon.  The  particular  types  of  weapons 
generally  employed  in  ASW  combat  are  the  following: 

1.  Hedgehog, 

2.  SUB ROC, 

3.  ASROC, 

4.  Torpedoes, 

5.  Aerial  bombs, 

6.  Artillery,  and 

7.  Mines. 


Once  the  weapon  has  been  launched  it  becomes  a partici- 
pant in  the  simulation,  possibly  possessing  a command  and  control 
function,  information  sensors,  and  decision  implementation  hard- 
ware, as  well  as  radiating  information  to  other  participants. 

The  simulation  of  the  performance  of  these  weapons  can 
be  divided  into  three  areas  once  the  command  to  fire  the  weapon 
has  been  given.  The  initial  area  of  simulation  is  the  launch 
phase,  during  which  the  weapon  accelerates  or  turns  to  reach 
a predetermined  trajectory.  The  second  area  of  simulation 
is  the  midcourse  phase,  during  which  new  trajectories  may  be 
achieved  if  necessary  (and  if  the  weapon  has  this  capability). 

The  third  area  is  the  terminal  (homing)  phase  in  which  the  weapon 
closes  with  the  target  and  explodes  if  triggered. 

Information  channels  may  be  altered  by  the  decision 
maker,  thus  changing  the  types  of  information  stored  in  the 
information  bank. 
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In  addition,  the  opening  of  an  information  channel  may 
change  the  information  radiated  by  the  participant  to  the  medium. 
For  example,  if  a participant  uses  his  active  sonar,  he  radiates 
a signal  which  can  be  received  by  other  participants,  thus  be- 
traying the  presence,  and  possibly  the  position,  of  the  first 
participant.  Therefore,  the  simulation  must  be  provided  with 
radiated  noise  as  a function  of  the  opening  or  closing  of  informa- 
tion channels,  when  applicable. 

As  was  the  case  with  sonar  performance,  the  information 
necessary  to  simulate  hardware  performance  in  the  implementation 
of  decisions  is  derived  from  data  analysis  or  from  a performance 
prediction  model. 
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4.  THE  FORWARD  AREA  PATROL  SIMULATION  MODEL 

In  this  chapter,  a detailed  description  of  the  factors 
considered  in  the  preparation  of  a digital  computer  program  to 
simulate  the  forward  area  patrol  mission  is  presented.  A pro- 
gramming level  description  of  the  computer  program  is  given  in 
the  Appendix.  This  program  was  prepared  to  assist  in  the  evalua 
t ion  of  the  effectiveness  of  certain  sonar  subsystems  proposed 
for  the  SS/SSN  phase  of  the  SISS  program.  (A  program  phase  to 
evaluate  the  effectivenss  of  certain  sonar  subsystems  proposed 
for  installation  on  the  Post  World  War  II  diesel-electric  and 
the  pre- SKIPJACK  nuclear  attack,  submarines). 

4.1  INTRODUCTION 

Briefly,  the  forward  area  patrol  mission  involves  the 
encounter  between  a patrol  submarine  (sometimes  called  the  Bar- 
rier Submarine)  assigned  to  patrol  a certain  geographic  region 
(the  barrier),  and  a transitor  (sometimes  called  the  Penetrator) 
The  mission  of  the  patrol  submarine  is  to  kill  any  submarine 
which  enters  his  assigned  patrol  region.  The  patrol  submarine 
must  not  leave  his  assigned  patrol  region.  The  mission  of  the 
transitor  is  to  traverse  the  barrier  region.  The  first  concern 
of  both  submarines  is  to  avoid  being  killed  themselves.  The 
patrol  submarine  is  provided  with  the  subsystems  under  considera 
tion  in  the  SS/SSN  program.  The  transitor  is  provided  with  only 
detection  and  bearing  determination  capability. 

Time  limitations  imposed  by  the  SS/SSN  study  made  it 
impossible  to  thoroughly  treat  every  aspect  of  the  problem  as 
discussed  in  Chapter  3.  The  simulation,  however,  addresses  all 
of  the  salient  features  of  the  SS/SSN  problem. 

For  example,  the  transitor  submarine  has  not  been 
provided  with  a weapon.  The  mission  of  the  transitor  is  to  pass 
through  the  barrier  and  hence  he  would  only  attack  the  patrol 
submarine  in  rare  instances.  Also,  reverberation  has  not  been 
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included  in  the  program.  The  high  sea  state  specified  for  the 
problem,  however,  implies  that  active  sonar  is  more  likely  to 
be  noise  limited  than  reverberation  limited. 

In  addition  to  the  aid  in  evaluating  the  effective- 
ness of  the  SS/SSN  subsystems,  this  model  illustrates  almost 
every  portion  of  the  general  ASW  combat  simulation  discussed 
in  Chapter  3.  Thus,  it  should  be  a valuable  foundation  upon 
which  to  construct  a more  general  and  more  complete  ASW  combat 
s imulation. 

The  organization  of  this  chapter  conforms  to  the 
organization  of  the  computer  program,  and  consists  of  four  main 
sections : 

1 . Program  Control  (CONTOL) , 

2.  Position  (POSIT), 

3.  Subsystems  (SUBSYS),  and 

4.  Command  and  Control  (TACTIC). 

CONTOL  is  the  input-output  portion  of  the  program  and  also  per- 
forms various  housekeeping  functions.  POSIT  simulates  the  motion 
of  the  submarines  and  corresponds  to  the  uppermost  box  in  Fig.  4. 
SUBSYS  updates  the  information  provided  by  each  subsystem  and 
corresponds  to  the  second  and  third  boxes  of  Fig.  4.  The  TACTIC 
portion  performs  all  other  functions  indicated  in  Fig.  4.  Thus, 
TACTIC  simulates  those  assessment  and  decision  processes  indi- 
cated schematically  in  Fig.  5.  A discussion  of  each  of  these 
four  sections  follows. 

4.2  CONTROL  (CONTOL) 

This  routine  has  three  major  functions,  namely, 

1.  Input  of  variables  and  initialization, 

2.  Control  of  the  simulation, 

3.  Output  of  variables  and  terminal  prints. 
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The  detailed  input  of  variables  is  described  in  the  Appen- 
dix. The  initialization  sets  all  constants  to  initial  values, 
and  stores  input  data  in  the  working  arrays. 

The  control  section  of  the  simulation  is  best  under- 
stood by  the  following  diagram  (Fig.  8)  of  the  program  flow 
directed  by  CONTOL: 


FIG.  8 PROGRAM  FLOW  DIRECTED  BY  CONTOL 

This  flow  cycle  is  continued  until  the  simulation  is 
completed,  then  the  output  of  variables  and  terminal  prints  (as 
desired)  are  effected.  The  various  nrint  options  and  sample 
output  are  described  under  CONTOL  in  the  Appendix  (A.l). 

It  should  be  noted  that  once  a weapon  has  been  ’aunched, 
control  of  the  simulation  is  currently  handled  by  the  Choice  II  B 
torpedo  model.  Control  is  returned  to  CONTOL  by  the  torpedo  model 
only  after  a probability  of  kill  has  been  calculated  and  the  simu- 
lation is  complete,  except  for  termination  prints. 

4.3  POSITION  (POSIT) 

The  subroutine  POSIT  is  entered  to  increment  the  motion 
of  each  submarine  in  the  simulation  model,  and  compute  the  new 
positional  information  for  each  submarine.  Before  discussing  the 
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motion  options  available  to  the  submarines  POSIT,  it  is  necessary 
to  discuss  the  coordinate  system  in  which  this  motion  is  des- 
cribed. 

4.3.1  Coordinate  System 

The  coordinate  system  used  in  the  model  is  an  earth- 
fixed  system  (Fig.  9).  The  origin  is  at  any  convenient  point 
on  the  ocean's  surface.  The  Z-axis  is  directed  radially  out- 
ward from  the  center  of  the  earth.  The  Y-axis  is  in  the  plane 
containing  the  Z-axis  and  the  earth's  axis  of  rotation,  is  normal 
to  the  Z-axis  and  is  directed  in  a northward  direction.  The 
X-axis  is  such  that  X,  Y,  and  Z form  a right  handed  coordinate 
system.  The  maximum  boundary  of  the  current  problems  considered 
is  small  compared  to  the  radius  of  the  earth;  thus,  effects  due 
to  the  curvature  of  the  earth  are  ignored.  The  X-Y  plane  is 
taken  to  be  the  horizontal  plane  and  the  Y-Z  plane  is  the  vertical 
p lane . 

It  should  be  noted  that  in  this  coordinate  system, 
depths  are  negative. 

4.3.2  Ship  Motion 

The  computer  model  is  capable  of  simulating  three 
types  of  ship  motion  at  present.  The  types  of  motion  are: 

(1)  constant  velocity,  (2)  constant  acceleration  and  (3)  circular 
horizontal  motion,  constant  speed.  The  equations  necessary  to 
simulate  these  types  of  motion  are: 

4. 3.2.1  Constant  velocity 

-4  -4  -4 

R (t.  + x)  = R (ti)  + R • [tt  + 1 - ti]  , (3) 

where 

-4 

R (t^  + is  the  position  vector  of  the  ship 
at  time  + p 
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R (t^)  is  the  position  vector  of  the  ship 
at  the  previous  time  point  t.  , and 

• X 

-4 

R is  the  constant  velocity. 

The  position  vectors  are  expressed  in  terms  of  their  cartesian 


/X(tl)\ 

R (tp  « y(t1)  ] • 


components 


\z(ti)/ 

/x(ti  + i}\ 

r (ti  + p -( Y(tt  + p y (5) 

W'i  + p/ 

The  velocity  vector  is  expressed  in  terms  of  its  vertical 
magnitude,  horizontal  magnitude,  and  the  angle  from  the  Y-axis 
to  the  horizontal  component  of  the  velocity  (see  Fig.  10), 

- /C  sin  cp  \ 

R = I C cos  co  | , (6) 


where 


is  the  magnitude  of  the  horizontal  com- 
ponent of  the  velocity, 

is  the  magnitude  of  the  vertical  com- 
ponent of  the  velocity,  and 

is  the  ship's  heading  (see  Fig.  10). 


4. 3. 2. 2 Constant  acceleration 


R (t.  + L)  = R (t.)  + R (t±)  • Cti  + - t.l 


+ h A ? [tt  + x - ttl2 
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FIG.  10  QUANTITIES  ASSOCIATED  WITH  THE 

POSITION  AND  MOTION  OF  THE  SHIPS 


vi-£  _j_  ^ ) ciiiu  are  as  ueriueu  j-u  uie  cu 

velocity  case, 

R (t*  , . ) is  the  velocity  of  the  ship  at  time 


R (t.)  is  the  velocity  of  the  ship  at  time  t.,  and 


A is  the  constant  acceleration 
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The  acceleration  vector  is  expressed  in  terms  of  the  vertical 
magnitude,  horizontal  magnitude,  and  the  angle  from  the  Y-axis 
to  the  horizontal  component  of  the  acceleration, 


a sin  cp\ 
a cos  cp  ] , 

Z / 


where 


a is  the  magnitude  of  the  horizontal  component  of 
the  acceleration, 

M 

Z is  the  magnitude  of  the  vertical  component  of 
the  acceleration,  and 

cp  is  the  ship's  heading. 

4. 3.2.3  Circular  horizontal  motion,  constant  speed 

, xc  + rc cos  ?7(ti  + 1 - ti)  + * (ti>  '■ 

- I C I 

R (ti  + 1>  = ' Yc  + rc  sin  ^(ti  + 1 " ti)  + * (ti}  ! ’ (10) 

' C I 

\ Z(t.)  + z{t.  + x - t.]  / 


where 


X , Y are  the  coordinates  of  the  center  of  the 
c c 

circle  of  motion, 

r is  the  radius  of  the  circle  of  motion, 

c 

Cc  is  the  speed,  Cc<0  produces  clockwise  motion, 

C > produces  counter  clockwise  motion, 
t.  is  the  time  at  some  time  point  of  the 

simulation, 
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+ | is  the  time  at  the  next  time  point, 

Z(t^)  is  the  depth  at  time  t^, 

Z is  the  vertical  speed, 

<lr  (t^)  = tan 


R (t  • ) - Y 
y v x' c 

R (t.)  - X 
x v l'  c 


and 


(ID 


R^(t^)  and  R^(t^)  are  the  X and  Y components  of 
R(t^)  respectively. 


The  ship's  heading  at  any  time  point,  t.,  is  given  by 

r_C  “■  1 


* (ti)  = 


tan 


-1 


~ [ R - X \ 
r V x c) 


C 

( y 

rc  1 c 


R 


(12) 


To  insure  that  co  (t . ) is  in  the  correct  quadrant,  it 
is  necessary  to  retain  the  quantity  (C  /r  ) in  both  the  numerator 
and  denominator  of  the  arc  tangent  argument. 


4.3.3  Quantities  Associated  with  Ship's  Motion  and  Position. 

Several  other  quantities  associated  with  ship  position 
and  ship  motion  are  required  in  the  simulation.  These  quantities 
are  illustrated  in  Figure  10,  and  the  formulas  for  their  computa- 
tion are: 

The  Vector  from  Ship  n to  Ship  m 

pnm  (ti)  = Rfti  (t^  - Rii  (t^,  (13) 

where 

-4 

Rm  (ti)  and  Rn  (t^)  are  the  position  vectors  of 
ship  m and  ship  n respectively. 

t 1 * 
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The  Horizontal  Vector  from  Ship  n to  Ship  m 


p Hnm  ( t ^ ) = 


pnm  (t^x" 
3 rim  (ti)y 

0 


(14) 


where 


onm  (t^)x  and  pnm  (t^)  y are  respectively  the 
magnitudes  of  the  X and  Y components  of  "pnm  (t^)  . 

The  Horizontal  Range  from  Ship  n to  Ship  m 

"*  r--*  "*  ~]h 

1 pHnm  (t^)  I = | pHnm  (t^)  • pHnm  (t^)J  • (15) 

The  True  Bearing  of  Ship  m with  respect  to  Ship  n 


cpTnm  (t^)  = tan 


-1 


pnm  (ti)x 
pnm  (ti)y 


(16) 


The  Relative  Bearing  of  Ship  m with  Respect  to  ship  n 

onm  (t^)  = r-oTnm  - cpn  , if  (mTnm-con)  ^ 0 
= 2-  + (mTnm-cpn) , if  (cpTnm-con)  < 0 

where 


(17) 


ccn  is  the  heading  of  ship  n. 

A detailed  discussion  of  the  variables  computed  by 
POSIT  is  contained  in  the  Appendix  (A. 2).  Once  all  positional 
data  is  computed,  the  subroutine  associated  with  updating  the 
subsystems  is  entered. 

4.4  SUBSYSTEM  (SUBSYS) 

The  simulation  model  requires  an  evaluation  of  the 
detection,  localization,  and  classification  performance  asso- 
ciated with  the  appropriate  subsystems.  The  particular  perfor- 
mance parameters  required  for  each  subsystem  are  listed  in 
Table  III.  Detection  performance  for  the  active  and  passive 
ranging  subsystems  has  not  been  determined  for  this  study. 
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TABLE  HI  PERFORMANCE  PARAMETERS  REQUIRED  FOR  EACH  SUBSYSTEM 


SUBSYSTEMS 

REQUIRED  PERFORMANCE  PARAMETERS 

AN/BQR-2B 

PROBABILITY  OF  DETECTION 
BEARING  ERROR 

AN/BQR-2(  ) DIMUS 

PROBABILITY  OF  DETECTI ON 
BEARING  ERROR 

PUFFS  (AN/BQG-4) 

RANGE  AND  BEARING  ERROR 

an/bqs-4 

RANGE  AND  BEARING  ERROR 

an/bqs-4  MODIFIED 

RANGE  AND  BEARING  ERROR 

AN/BQH-2C 

PROBABILITY  OF  CLASSIFICATION 
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The  environmental  conditions  for  which  the  performance 
data  apply  are  those  stated  in  SOR  23-29.  They  are  based  upon 
statistical  data  from  the  Norwegian  Sea,  Summer  Season,  and  are 
as  follows: 


Summer  Season; 


Sea  State  Four; 

2000  Fathoms  of  Water; 

Average  Bottom  Loss  of  15  dB. 

4.4.1  Passive  Detection  Performance 

Detection  performance  data  in  terms  of  probability 
of  detection  in  both  the  convergence  zone  and  the  direct  path 
zone  are  required  for  the  AN/BQR-2B  and  AN/BQR-2  (DIMUS) 
subsystems.  These  data  have  been  generated  from  a passive  per- 
formance prediction  model  in  which  a certain  assumption  has  been 
made  in  the  absence  of  directly  applicable  data  on  display- 
operator  performance.  The  assumption  integrates  the  results  of 
two  special  purpose  experiments  to  yield  a more  general  model  of 
display-operator  performance  through  which  probability  of  detection 
can  be  determined  for  a specified  signal-to-noise  ratio  and  number 
of  target  lines  on  the  bearing- time  recorder  (BTR)  display.  The 
validity  of  the  model  cannot  be  determined  quantitatively  because 
of  the  lack  of  data  by  which  to  check  the  display-operator  model. 

It  is  thought  to  result  in  a reasonable  assessment  of  absolute 
performance  for  passive  detection  systems  having  a BTR  display 
and  to  be  quite  good  for  determining  the  relative  performance 
of  two  or  more  passive  detection  systems. 

4.4. 1.1  The  Static  performance  prediction  model.  The  perfor- 
mance prediction  model  applies  to  a static  situation  in  which 
probability  of  detection  is  computed  for  a fixed  target  signal- 
to-noise  level  and  a specified  number  of  target  information  lines 
on  the  BTR.  To  apply  this  static  model  to  the  dynamic  situation 
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in  which  target  range  can  change  continuously,  and  in  which  no 
meaningful  fixed  observation  time  corresponding  to  a fixed  number 
of  lines  on  the  display  can  be  established,  a procedure  has  been 
adopted  through  which  probability  of  detection  is  determined  as 
a function  of  time.  This  procedure  will  be  discussed  after  the 
static  model  has  been  described. 

The  static  prediction  model  includes  two  basic  steps: 

(1)  Prediction  of  the  signal  and  noise  levels  at 
the  beamformer  output  for  the  assumed  environmental  conditions; 
and 

(2)  Prediction  of  probability  of  detection  for  a 
specified  number  of  target  information  lines  on  the  BTR  at  the 
beamformer  output  signal-to-noise  ratio  determined  in  step  (1). 

A block  diagram  of  the  process  is  shown  in  Fig.  11. 

The  prediction  of  all  signal-to-noise  levels  (step  1) 
in  this  program  is  handled  by  a separate  routine  termed  STON.  As 
this  routine  serves  localization  and  classification  as  well  as 
detection  subsystems,  it  will  be  described  at  the  end  of  the 
subsystem  descriptions. 

Step  (2)  accounts  for  subsystem  performance  from 
Point  A to  Point  C.  Point  B is  of  interest  in  the  subsequent 
discussion  of  step  (2). 

The  performance  of  each  of  the  passive  subsystems  can 
be  presented  as  a set  of  probability  of  detection  curves,  one 
set  for  each  of  the  passive  detection  subsystems.  The  sets  of 
curves  which  have  been  generated  for  the  AN/BQR-2B  and  AN/BQR-2 
(DIMUS)  subsystems  are  presented  in  Figures  12  and  13  res- 
pectively. Once  the  signal-to-noise  ratio  of  the  input  to  the 
signal  processor  has  been  determined,  and  the  number  of  lines 
on  the  BTR  at  this  signal-to-noise  ratio  has  been  specified, 
then  the  probability  of  detection  can  be  obtained  directly  from 
the  set  of  curves. 
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In  order  to  determine  a set  of  probability  of  de- 
tection curves,  the  results  of  two  display-operator  studies,  one 
by  Deatherage  and  one  reported  in  TDP  S23-19,  are  used  in  con- 
junction with  the  proper  treatment  of  the  signal  processor. 

Deatherage  investigated  a 50-line  display.  He  found 
that  for  fixed  Gaussian  signal  in  Gaussian  noise  at  a particular 
bearing  (1)  probability  of  detection  is  maximum  when  the  noise 
marking  density  is  between  0.5  and  0.6,  and  (2)  if  the  noise 
marking  density  is  set  at  0.5,  the  probability  of  detection  de- 
pends upon  the  signal-plus-noise  marking  density  according  to 
the  curve  given  in  Fig.  14.  The  study  reported  in  TDP  S23-19 
resulted  in  a graph  of  time  to  detect  the  target  with  0.5  prob- 
ability versus  required  signal- to-noise  ratio  at  the  input  to 
the  AN/BQR-2B  DIMUS  signal  processor. 

In  order  to  utilize  the  results  of  these  display- 
operator  studies,  it  is  necessary  to  determine  the  processing 
gain  of  the  output  signal-plus-noise  statistics  of  the  signal 
processor.  A summary  of  techniques  used  to  determine  these  sig- 
nal processor  characteristics  can  be  found  in  P.  B.  Brown's 
summary  on  signal  processors+.  The  procedure  carried  out  consists 
of  two  steps  : 

(1)  Determine  the  processing  gain  of  the  signal 
processor;  this  is  a plot  of  the  processor  output  signal-to- 
noise  ratio,  (S/N)g,  plotted  as  a function  of  the  input  signal- 
to-noise  ratio,  (S/N)^.  (S/N)^  is  an  average  signal  power  to 

average  noise  power  ratio.  (S/N)g  is  the  square  of  (p  - u)/a, 


Deatherage,  B.  H. , "Summary  Report  of  Results,  Conclusions  and 
Recommendations  from  a Psychophysical  Study  of  the  Relative  De- 
tectability of  Target  Tracks  in  Simulated  Passive  Sonar  Dis- 
plays", TRACOR  Document  63-231-U  prepared  under  Contract  NObsr- 
89265  (Sept. 6,  1963). 

4* 

P.  B.  Brown,  "A  Comparison  of  the  Performance  of  Several  Signal 
Processors",  TRACOR  Document  66-203-U,  prepared  under  Contract 
N123 (953)  53354A-NPOLA,  (March  1966). 
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where  p is  the  output  signal-plus-noise  peak,  u is  the  noise- 
only  mean  at  the  output,  and  a is  the  standard  deviation  of  the 
noise-only  output. 

(2)  Prepare  plots  of  the  signal-plus-noise  proba- 
bility of  exceeding  threshold  as  a function  of  threshold,  using 
output  signal-to-noise  ratio  as  a parameter.  These  curves  (one 
for  each  signal-to-noise  ratio)  are  then  relabeled  using  the  in- 
put signal-to-noise  ratio,  obtained  from  the  processing  gain 
curve  determined  in  step  (1),  corresponding  to  specific  output 
signal-to-noise  ratios  used  in  preparing  the  curves.  The  above 
procedure  has  been  carried  out  for  the  AN/BQR-2B  and  the  AN/BQR-2B 
DIMUS  signal  processors. 

At  this  point,  it  is  possible  to  determine  probability 
of  detection  for  a 50-line  display  as  a function  of  signal-to- 
noise  ratio  into  the  corresponding  signal  processor.  Using  the 
results  of  Deatherage's  study,  it  is  assumed  that  the  display 
threshold  is  set  for  a 0.5  noise  only  marking  density.  Since 
the  probability  that  the  signal  out  of  the  processor  exceeds  a 
specified  threshold  (zero  in  this  case)  is  equivalent  to  marking 
density  on  the  display,  the  curve  of  Figure  14  can  be  replotted 
as  probability  of  detection  versus  signal-to-noise  ratio  into 
the  signal  processor.  This  is  done  by  using  the  signal  processor 
characteristics  to  determine  the  (S/N)^  at  zero  threshold  which 
corresponds  to  a given  signal-plus-noise  marking  density.  The 
results  of  following  this  procedure  for  both  the  AN/BQR-2B  and 
the  AN/BQR-2B  DIMUS  signal  processors  are  shown  in  Figure  15. 

These  curves  yield  the  required  information  for  a 50-line  dis- 
play. 

The  results  of  the  study  presented  in  TDP  S23-19  are 
used  in  the  prediction  of  probability  of  detection  for  numbers 
of  display  lines  other  than  50.  The  graph  of  time  to  detect 
with  0.5  probability  versus  required  signal-to-noise  ratio  into 
the  AN/BQR-2  (DIMUS)  signal  processor  can  be  replotted  to  reflect 
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operator  display  performance  only  by  properly  accounting  for  the 
signal  processor  characteristics.  The  "time-to-detect"  axis  can 
be  converted  to  number  of  lines  on  the  display,  and  the  "signal- 
to-noise  ratios  into  the  signal  processor"  axis  can  be  converted 
to  signal-to-noise  ratio  into  the  display.  This  procedure  re- 
sults in  the  graph  shown  in  Figure  16.  Although  this  graph  was 
measured  for  0.5  probability  of  detection,  it  is  assumed  to  re- 
present the  functional  relationship  between  required  signal-to- 
noise  ratio  and  number  of  lines  on  the  display  for  probability 
of  detection  between  0.1  and  0.9.  The  sets  of  probability  of 
detection  curves  given  in  Figures  12  and  13  are  established  by 

(1)  determining  the  required  input  signal-to-noise  ratio  at 

0.1,  0.2,  0.3,....,  0.9  probability  of  detection  from  Figure  15; 

(2)  relabeling  the  signal-to-noise  ratio  axis  of  Figure  16  to 
give  signal-to-noise  ratio  into  the  signal  processor;  and  (3) 
shifting  the  50-line  point  on  the  curve  of  Figure  15  to  coin- 
cide with  these  input  signal-to-noise  ratios. 

4.4. 1.2.  Procedures  for  obtaining  Detection  Performance  in 
the  Dynamic  Operational  Condition.  The  static  performance  data 
discussed  in  the  previous  subsection  do  not  supply  the  perfor- 
mance information  required  for  the  combat  simulation  model.  In 
the  model  it  must  be  ascertained  at  all  times  whether  each  parti- 
cipant has  contact  with  the  other  participant. 

Furthermore,  in  the  dynamic  condition,  the  signal-to- 
noise  ratio  out  of  the  beamformer  changes  continuously  with  the 
target  range.  This  varying  signal-to-noise  ratio  is  reflected 
on  the  sonar  display  as  a non-uniform  signal-plus-noise  marking 
density.  To  apply  the  static  prediction  model  to  this  dynamic 
condition  a procedure  has  been  adopted  through  which  a signal- 
to-noise  field  is  calculated  at  each  time  step  of  the  program. 

The  values  in  the  signal-to-noise  field  are  then  tested  to  yield 
the  greatest  probability  of  detection  for  that  time  step. 
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This  procedure  is  based  on  the  assumption  that  the 
operator  will  consider  in  effect  any  recent  increase  in  marking 
density  to  the  exclusion  of  what  occurred  before  the  increase. 

This  assumption  is  implemented  in  the  following 
manner.  At  each  time  step,  the  signal-to-noise  ratio  at  the 
output  of  the  beamformer  is  computed  for  the  detection  subsystem. 
If  the  signal-to-noise  ratio  is  less  than  a threshold  value 
(the  signal-to-noise  ratio  corresponding  to  a 0.10  probability 
of  detection  for  an  infinite  observation  time),  a value  of  -« 
on  the  log  scale  is  stored  in  the  signal-to-noise  table  for  that 
time  step.  This  signal-to-noise  table  represents  the  display 
screen  seen  by  the  operator.  If  the  table  is  full,  the  oldest 
sigi al-to-noise  ratio  is  dropped  to  make  room  for  the  most  re- 
cent data. 

Two  possible  situations  must  be  considered  at  this 
point.  The  first  is  that  the  target  has  been  detected,  and  it 
must  be  determined  if  detection  has  been  lost.  The  second  situa- 
tion is  that  the  target  has  not  been  detected,  and  it  must  be 
determined  whether  or  not  detection  occurs  during  this  time  step. 
These  situations  will  be  discussed  separately  in  the  following 
text . 

If  the  signal-to-noise  ratio  is  below  the  threshold 
for  five  consecutive  time  steps,  detection  is  assumed  lost,  the 
table  of  signal-to-noise  ratios  is  cleared,  and  the  detection 
process  restarted. 

A conditional  probability  of  detection  is  required 
to  determine  if  detection  has  occurred  at  any  particular  instant 
of  time.  This  conditional  probability  of  detection  is  defined 
to  be  the  probability  of  detection  at  that  instant  multiplied  by 
the  probability  of  failure  to  detect  prior  to  the  instant  in 
question.  Once  the  conditional  probability  of  detection  is 
determined  the  question  of  whether  detection  has  occurred  can 


be  treated  with  standard  Monte  Carlo  techniques.  In  the  procedure 
used  here,  the  conditional  probability  of  detection  is  not  ex- 
plicitly determined;  rather,  a Monte  Carlo  technique  is  applied 
to  the  static  detection  probability  in  such  a way  that  the  static 
detection  probability  is  effectively  transformed  into  the  condi- 
tional probability  of  detection.  The  specific  procedure  is  as 
follows.  After  a value  has  been  stored  in  the  signal-to-noise 
table,  the  values  in  the  table  are  averaged  in  the  following 
manner.  Let  the  total  number  of  values  in  the  Table  be  K,  the 
number  of  time  steps  for  which  information  has  been  displayed 
to  the  operator.  First,  an  average  signal-to-noise  is  computed 
for  the  K values  in  the  Table.  (Recall  that  K corresponds  to 
the  length  of  time  over  which  the  operator  observes  the  target, 
which  is  equivalent  to  some  number  of  lines  of  information  on 
the  BTR  display).  The  average  value  of  the  signal-to-noise  and 
the  number  N are  then  used  to  compute  the  probability  of  detec- 
tion from  the  curves  diagrammed  in  this  section.  (The  computer 
program  actually  uses  a table  look-up  and  linear  interpolation). 
The  averaging  process  is  repeated  for  the  latest  (K-l)  values 
in  the  signal-to-noise  table.  The  probability  of  detection  is 
again  computed  from  the  curves  using  the  new  average  signal-to- 
noise  and  the  number  (K-l).  If  the  latest  probability  of  de- 
tection is  greater  than  that  previously  computed,  it  is  stored 
in  the  probability  of  detection  location.  The  process  is  re- 
peated for  (K-2) , (K-3),  . . .,  1,  until  the  maximum  value  of 
the  probability  of  detection  has  been  determined.  We  refer  to 
this  maximum  value  as  the  static  probability  of  detection  at 
that  time  step. 

The  probability  of  detection  is  then  compared  to  a 
detection  limit  which  is  initially  0.1.  If  it  is  less  than  this 
limit,  detection  is  assumed  impossible,  and  no  test  for  detection 
is  made. 
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occurred  at  this  time  step  is  to  select  a random  number  between 
0.0  and  2.0.  If  the  number  chosen  is  less  than  the  probability 
of  detection  for  that  time  step,  detection  is  said  to  have 
occurred.  If  not,  the  detection  limit  is  increased  by  0.05,  to 
a maximum  of  0.50. 

The  selection  of  a random  number  between  0.0  and  2.0 
rather  than  between  0.0  and  1.0  is  the  process  by  which  the 
static  probability  of  detection  is  transformed  into  conditional 
probability  of  detection.  The  selection  of  the  interval  0.0  to 
2.0  was  determined  by  numerical  experimentation  to  be  the  most 
appropriate  sampling  interval  based  upon  physical  intuition. 

The  validity  of  this  technique  must  be  established  by  experi- 
mentation. 

The  increase  in  the  detection  limit  is  instigated  in 
the  belief  that  if  an  observer  cannot  detect  a target  at  one 
time  step  for  a given  probability  of  detection,  he  will  detect 
nothing  in  succeeding  time  steps  until  the  probability  increases 
by  some  amount  (say  0.05).  Once  the  limit  reaches  0.50,  no 
further  increase  in  the  limit  occurs. 

4.4.2  Localization  Subsystem  Performance 

Any  localization  subsystem  can  be  characterized  by 
the  following  diagram  (Fig.  17): 


FIG.  17.  CHARACTERIZATION  OF  LOCALIZATION  SUBSYSTEM 
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The  input  data  is  the  actual  range  or  bearing  of  the  target.  The 
error  is  approximated  as  the  sum  of  a bias  error  and  an  error 
randomly  selected  from  a normal  distribution  of  standard  deviation  a 
and  mean  0.  It  is  therefore  necessary  to  specify  the  bias (mean) 
error  and  the  standard  deviation  for  each  localization  subsystem, 
either  from  a theoretical  formula  or  by  a fit  to  available  data. 


We  distinguish  between  ranging  and  bearing  systems 
in  the  following  discussion,  rather  than  distinguishing  on  the 
basis  of  equipment.  Hence,  PUFFS  performs  as  both  a bearing 
and  ranging  system.  The  bearing  data  has  one  characteristic 
error,  the  ranging  data  another.  PUFFS  is  therefore  treated 
in  the  ranging  section  and  in  the  bearing  section  of  this  program. 


4.4.2. 1 Ranging  Systems 

PUFFS . According  to  the  MK48  Torpedo  Simulation  Model 
CHOICE  II-B , the  theoretical  range  error  of  the  PUFFS  subsystem 
sonar  is 


: 3 VT~  (RT)2 
R Sf3  sin  B 


(18) 


where 


= range  error  in  yards ; 

RT  = range  in  yards  to  target ; 

Sg  = receiver  spacing  in  feet; 

B = target  bearing  in  degrees; 
aB  = bearing  error  in  radians,  due  to  noise. 


Greenstone,  R. , and  Tully,  J.  J. , The  Effect  of  Torpedo  Radiated 
Noise  on  the  MK48  MOD  0 Torpedo  Weapon  System,  Part  5,  TM 
412.2711-^7,  6 "May  1966', (Contract'  NOw'6’5- 0123-37: 
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°B  = 


0.113 

sr t 


Q + 0.19  Q" 


°B 


error  in  degrees,  due  to  system  noise, 


T = integration  time; 

= 20.6667  + 311111  Q*; 

Q = Q in  dB; 

Q = noise-to-signal  ratio  in  power. 

An  attempt  was  made  to  compare  the  results  of  this 
formula  with  certain  measured  range  errors  obtained  at  Dabob 
Bay  during  TECHEVAL  and  0PEVAL  on  the  AN/BQG-4.  The  submarine 
equipped  with  PUFFS  was  moving  at  a very  slow  speed. 

It  was  first  assumed  that  the  noise  source  was 
equivalent  to  a 15  knot  ship  and  that  the  background  noise 
corresponded  to  Sea  State  2.  Using  the  generalized  radiated 
noise  curve,  the  following  errors  were  computed  with  the  above 
theoretical  formula. 


Range 


Error 


5000  yards 
8000  yards 
10,000  yards 


44  yards 
450  yards 
1350  yards 


If  the  substituted  conditions  are  assumed  to  be 
similar  to  those  existing  during  the  OPEVAL,  then  the  data  above 
can  be  compared  to  the  data  taken  on  the  USS  VOLADOR.  Between 
the  ranges  of  7000  - 9000  yards,  the  standard  deviation  of  the 
VOLADOR  data,  based  on  417  samples,  was  500  yards.  The  standard 
deviation  of  32  samples  taken  at  8000  yards  was  471  yards.  The 
subsurface  target  used  in  the  OPEVAL  was  the  RONQUIL,  snorkeling, 
cavitating,  and  circling  VOLADOR  at  a range  of  5000  yards.  Al- 
though many  unknowns  were  present,  indications  were  that  the  PUFFS 
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equation  in  the  CHOICE  II-B  model  was  in  reasonable  agreement 
with  measured  data  at  ranges  out  to  8000  yards. 

Further  study  indicated  that  the  theoretical  range 
error  was  strongly  dependent  on  the  value  of  signal-to-noise 
ratio.  For  example,  when  signal  is  assumed  to  correspond  to 
a '3  knot  ship,  the  following  errors  are  calculated. 

5000  yards:  oR  = 166  yards; 

8000  yards:  oR  = 1390  yards; 

10,000  yards:  aR  = 4180  yards. 

The  errors  are  both  large  and  in  considerable  disagreement  with 
the  experimental  results  at  long  ranges.  Therefore,  since  the 
true  experimental  value  of  signal-to-noise  was  not  known  in  any 
of  the  test  cases,  it  must  be  concluded  that  an  accurate  experi- 
mental verification  of  the  theoretical  formula  has  not  been  made. 

The  standard  deviation  of  all  of  the  experimental  data 
could  be  fitted  rather  well  by  the  following  formula  (outliers 
removed) : 

aR  = 51.7  + 19.4  R + 2.68  R2  , (20) 

where  aR  is  in  yards,  and  R is  in  kiloyards. 

2 

This  fit  indicates  verification  of  the  R dependence 
of  the  range  error  which  is  portrayed  in  the  PUFFS  equation.  The 
mean  (or  bias)  error  of  PUFFS  is  assumed  to  be  zero. 

AN/BQS-4  and  AN/BQS-4  Modified.  The  data  base  used 
to  assess  the  range  accuracy  of  the  AN/BQS-4  and  the  AN/BQS-4 
Modified  was  data  taken  at  Dabob  Bay  and  at  Fleet  Operational 
Readiness  Accuracy  Check  Sites  (FORACS)  at  Guantanamo  Bay,  Cuba, 
and  San  Clemente  Island,  California.  Indications  from  the  limited 
data  available  were  that  the  present  AN/BQS-4  could  be  assigned 
a range  error  standard  deviation  of  about  47o.  Other  active  sonars, 
such  as  the  AN/SQS-23  series,  have  exhibited  a range  error  of 
about  17,  on  these  same  measurement  ranges.  These  data  are  generally 
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taken  at  ship  speed  of  3 knots  and  in  Sea  State  2 or  less  with 
the  ship  going  into  or  with  the  sea.  The  result  of  the  examina- 
tion was  that  the  AN/BQS-4  range  accuracy  was  estimated  to  be 
about  4%  and  the  AN/BQS-4  modification  was  chosen  at  27,.  The 
range  accuracy  assigned  to  the  AN/BQS-4  modified  was  based  on 
the  fact  that  the  AN/SQS-23  had  demonstrated  that  17,  accuracy 
had  been  achieved  and  that  the  modification  to  AN/BQS-4  should 
be  comparable  to  the  present  AN/SQS-23  in  measuring  range. 

Hence,  the  AN/BQS-4  was  assumed  to  have  a standard 
deviation  of  0.04R  and  a zero  mean  error  while  the  AN/BQS-4 
modified  was  assumed  to  have  a standard  deviation  of  0.02R  and 
a zero  mean  error. 

4. 4. 2. 2 Bearing  Systems 

Two  sources  of  bearing  information  are  under  con- 
sideration: The  AN/BQR-2B,  the  AN/BQR-2  (DIMUS)  and  the 

PUFFS  sonar  subsystem. 

AN/BQR-2B  and  AN/BQR-2  (DIMUS!.  The  assessment  of 
bearing  errors  supplied  by  the  AN/BQR-2B  or  the  AN/BQR-2  ( ) 
DIMUS  was  made  on  data  similar  to  that  used  for  estimation  of 
range  errors.  The  AN/BQR-2B  has  exhibited  a peak  local  mean 
error  of  about  1°  at  Dabob  Bay  and  FORACS.  This  error  includes 
a bias  error  (alignment)  and  a bearing  dependent  error.  The 
standard  deviation  about  a best  fit  curve  through  the  bearing 
dependent  bearing  errors  was  about  0.3°  for  807,  of  the  equip- 
ments tested  at  Dabob  Bay.  Therefore,  the  AN/BQR-2B  was  assumed 
to  have  a bias  of  +1°  and  a standard  deviation  of  0.3°.  The 
AN/BQR-2  (DIMUS)  data  was  so  scarce,  that  it  was  assigned  a 
bias  error  of  zero  and  a standard  deviation  of  0.25°. 

PUFFS . According  to  the  MK48  Torpedo  Simulation 
Model  Choice  IIB,  the  theoretical  standard  deviation  of  the 
bearing  error  of  the  PUFFS  subsystem  sonar  is: 
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where 


(0.25)2  + oB2 


= bearing  error  in  degrees; 


cfg  = error  due  to  system  noise  (see  PUFFS 
range  error) . 

The  bias  error  of  PUFFS  is  assumed  to  be  zero. 

Each  localization  subsystem  is  assumed  to  have  a 
S/N  threshold  level,  below  which  the  system  will  return  no 
information.  This  cutoff  is  the  S/N  level  associated  with  the 
0. 1 probability  of  detection  for  infinite  observation  time  for 
the  AN/BQR-2B  and  2 (DIMUS) ; the  S/N  level  corresponding  to 
a 1 min.  signal  integration  time  for  PUFFS;  and  the  single 
ping  0.5  probability  of  detection  S/N  level  for  the  AN/BQS-4 
and  AN/BQS-4  (MODIFIED). 

Further,  each  localization  subsystem  is  assumed  to 
operate  only  in  a certain  sector  of  a 360°  circle.  The  sector 
in  which  the  system  is  able  to  return  no  useful  information  is 
termed  the  blind  spot  of  the  system.  This  blind  spot  occurs 
due  to  the  physical  location  of  the  sonar  system  in  question 
(i.e.,  a sonar  located  on  the  nose  of  the  ship  cannot  detect 
a target  aft).  All  systems  except  PUFFS  are  assumed  blind  in 
a 30°  sector  astern. 

The  PUFFS  system  presents  a special  case,  as  it  is 
mounted  atop  the  ship  and  is  designed  to  look  abeam.  Data 
gathered  more  than  40°  from  abeam  is  extremely  error  prone 
(see  formulas  for  error  in  this  section).  Therefore,  the 
blind  spot  of  PUFFS  consists  of  the  sectors  outside  ± 40°  from 
abeam.  Also,  PUFFS  is  assumed  to  have  a range  limit  beyond  which 
it  provides  no  useful  information.  From  the  OPEVAL  data,  this 
range  was  chosen  as  18,000  yds. 
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4.4.3.  Classification  Systems 

The  classification  system  is  the  AN/BQH-2C,  which 
serves  as  an  information  collection  system  in  conjunction  with 
the  AN/BQR-2B  or  the  AN/BQR-2  (DIMUS)  sonar  systems.  The 
information  available  on  this  classification  system  is  scarce  , 
but  the  following  approach  was  selected  for  this  study. 

If  the  signal-to-noise  level  at  the  input  to  the 
BQH-2C  is  at  least  -6dB  for  a period  of  20  min. , there  is  a 
0.90  probability  of  correct  classification.  It  is  assumed 
that  if  there  is  a 0.90  probability  of  correct  classification, 
then  the  target  has  been  correctly  classified. 

4.4.4  Signal-To-Noise  Calculation 

The  calculation  of  the  signal-to-noise  ratio  at  the 
beamformer  output  is  of  fundamental  importance  to  subsystem 
performance . 

The  calculation  of  signal-to-noise  is  divided  into 
two  basic  calculations  - one  for  passive  subsystems  and  one  for 
active  systems.  The  similarity  of  the  calculations  will  be 
apparent. 

The  quantities  which  determine  the  signal  and  noise 
levels  at  the  beamformer  output  for  a passive  sonar  are  defined 
below. 

(1)  T(f)  denotes  the  spectral  density  of  the  acoustic 
energy  radiated  by  the  target  referred  to  one  yard  from  the 
target ; 

(2)  H (f ) denotes  the  spectral  propagation  loss  of 
the  acoustic  energy  which  leaves  the  target  and  travels  to  the 
receiving  array.  This  quantity  includes  the  effects  of  multi- 
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path  arrivals  and  the  beam  deviation  loss  for  each  arrival; 

(3)  N(f)  denotes  the  spectral  density  of  the  environ- 


mental noise  field; 


(4)  A(f)  denotes  the  frequency  response  of  the  system 
filters;  and 


(5)  A(f)  denotes  the  array  directivity  index  for  the 
isotropic  environmental  noise  field. 


The  noise  power,  N , at  point  A in  the  block  diagram 

3 


is  calculated  from 


fo  + B 


Na  = J A (f )N(f )A (f)  df  , 


(22) 


in  which  f is  the  lower  band  edge,  and  B is  the  receiving  band- 


width. Since  the  input  equalizer  filters  for  the  subsystems  are 
usually  designed  to  achieve  maximum  "prewhitening,"  the  inte- 
grand above  i$  approximately  constant.  This  being  the  case, 
the  result  of  the  integration  is: 

(23) 


Na  = B A(fc)N(fc)A(fc)  , 


for  which  f is  the  effective  band  center,  defined  as  C f 0 ( f + B)] 


To  obtain  the  signal  power  S at  point  A in  the  block 

3 


diagram,  it  is  necessary  to  evaluate  the  integral 


fo  + B 

Sa  = j A (f ) H (f ) T (f ) df  . 


(24) 


o 


The  integral  has  been  approximated  by 


fo  + B 


Sa  - [b  A(fc)  T(fc)j  fl  j H(£)  df 

f 


(25) 
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The  second  bracket  in  the  above  equation,  denoted  by 
Le , represents  the  effective  transmission  loss  at  the  beamformer 
output  averaged  over  the  subsystem  bandwidth.  This  quantity  is 
computed  through  a numerical  integration  over  H(f).  The  spectral 
propagation  loss  function,  H(f),  is  computed  for  a particular 
velocity  profile  as  a function  of  range,  source  depth,  and  re- 
ceiver depth,  using  a ray  tracing  propagation  loss  model  in 
which  the  directivity  pattern  of  the  subsystem  is  utilized  to 
account  for  the  beam  deviation  loss  associated  with  the  direction 
of  arrival  of  each  ray.  Examples  of  the  spectral  propagation 
loss  curves  used  in  the  analysis  are  shown  in  Figures  18  through 

20.  The  beam  patterns  were  computed  at  the  subsystem  center 

. * 

frequency,  f , using  the  actual  subsystem  array  configuration  . 
The  particular  velocity  profile  utilized  by  the  propagation  loss 
model  for  the  Norwegian  Sea  is  given  in  Figure  21. 


The  ratio  of  average  signal  power,  S , to  average 

a 

noise  power,  N , denoted  by  (S/N)  , has  been  determined  using 


the  equation: 


a’ 

Sa  _ T (f  ) L 
(S/N)  = c e 


N 


N<V  “V 


(26) 


It  is  more  convenient  to  represent  the  signal-to- 
noise  ratio  in  dB,  in  which  case  the  equation  becomes: 

10  log  (S/N)a  = 10  log  T(fc)  + log  Le 

- 10  log  N(fc)  - 10  log  *(fc). 

(27) 


SS/SSN  Sonar  Improvement  Program  Analysis  (U)  TRACOR  Document 
Number  RL/67-025-C  of  May  1967  (CONFIDENTIAL). 
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FIG.  18  - PROPAGATION  LOSS  CURVES  - BELOW  LAYER 


PATROL  - ABOVE  LAYER 
TARGET  - BELOW  LAYER 
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FIG.  19-  PROPAGATION  LOSS  CURVES- CROSS  LAYER 


BANDWIDTH  - 1-4  KC 
PATROL  - ABOVE  LAYER 
TARGET  - ABOVE  LAYER 


FIG.  20  - PROPAGATION  LOSS  CURVES  — ABOVE  LAYER 
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The  target  radiated  noise  levels  [10  log  T ( f c ) ] and 

the  subsystem  directivity  indices  T 1 0 log  L (f  )]  used  in  the 

* c 

analysis  are  listed  in  another  report  . The  environmental 
noise  levels  T 10  log  N(fQ)l  were  obtained  by  adding  the  self- 
noise values  given  in  that  report  to  the  sea  state  four  noise 
level  at 

the  depth  of  the  receiver). 

The  calculation  of 
very  similar  to  the  passive  case,  with  the  exception  that 
the  signal  energy  is  transmitted  at  essentially  a single  fre- 
quency which  is  Doppler  shifted  when  received,  but  not  particu- 
larly broadened. 

The  quantities  which  determine  the  signal  level  at 
the  beamformer  output  for  active  sonars  are  defined  below. 


f given  by  the  Knudsen  curves  (after  correction  for 


(S/N)  for  an  active  sonar  is 

cl 


(1) 


S(fc)  denotes  the  spectral  density  of  the 


acoustic  energy  emitted  by  the  active  sonar  referred  to  one 
yard  from  the  emitter; 

(2)  AT(fc)  denotes  the  array  directivity  index  for 
the  transmitted  signal; 


(3) 


TS(fc,9)  denotes  the  target  strength  as  a 


function  of  the  angle  of  arrival  9. 

The  ratio  of  the  signal  power  to  the  average  noise 
power  has  been  determined  (in  dB)  using  the  following  equation: 

10  log  (S/N)  = 10  log  S(f  ) + 10  log  AT(f  ) + 10  log  TS(f  ,9) 

- 10  log  N(f  ) - 10  log  A(fc)  + 20  log  Le  - 10  log  B. 

(28) 


ibid 
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The  source  levels  T 10  log  S ( f c ) "]  , the  transmission 
directivity  indices  T 10  log  AT(fc)1,  the  receiving  directivity 
indices  I"  10  log  &(f  )]  and  the  bandwidth  of  the  receiver  T 10  log  B] 
have  been  listed  ' for  each  subsystem.  The  target  strength  data 
r 10  log  TS(f  ,0)1  was  not  available  as  a function  of  aspect  angle  9. 
Hence,  the  target  strength  was  approximated  as  a constant  for  all 
aspect  angles: 

10  log  TS(f  ) = 15  dB  . (29) 

It  should  be  observed  that,  in  the  computer  program, 
propagation  loss  is  stored  as  a negative  number,  receiving  direct- 
ivity indices  are  negative,  and  depths  are  negative.  Further, 
there  are  three  propagation  loss  tables  stored  for  each  subsystem. 
The  tables  represent  the  cases  for  both  ships  above  the  layer, 
one  above  and  one  below,  and  both  ships  below  the  layer. 

The  above  discussion  is  based  on  the  assumption  that 
the  effect  of  the  beamformer  on  the  acoustical  signal  field  and 
the  effect  on  the  acoustical  noise  field  can  be  considered  in- 
dependently. This  is  true  for  a linear  beamformer.  In  the  case 
of  a non-linear  beamformer  (a  DIMUS  system,  for  example),  this 
independent  treatment  of  signal  field  and  noise  field  to  predict 
the  (S/N)  ratio  at  the  output  of  the  beamformer  may  not  be  valid. 
This  problem  is  currently  under  investigation. 

4.4.5  Tracking  Subsystems 

In  the  simulation  of  the  Penetrator's  course  through 
the  barrier,  it  is  assumed  that  he  will  perform  some  random 
evasive  maneuvers  about  a predetermined  straight  line  course. 
Therefore,  there  are  two  types  of  tracking  problems  for  the 
patrol  submarine.  The  first  tracking  problem  is  to  determine 
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the  base  course  about  which  the  Penetrator  is  maneuvering.  This 
is  a long  term  problem  in  which  the  Penetrator  must  be  tracked 
for  several  legs  of  his  maneuver  course. 

The  second  tracking  problem  arises  in  conjunction 
with  the  launch  of  a weapon  at  the  target.  In  this  case,  the 
patrol  submarine  is  interested  in  the  present  course  and  speed 
of  the  Penetrator,  rather  than  the  base  course.  This  second 
problem  is  a short-term  tracking  problem  in  which  the  track  is 
reinitialized  each  time  the  Penetrator  is  observed  to  turn. 

A satisfactory  approximation  to  the  resolution  of 
both  tracking  problems  is  obtained  by  making  a linear  hypothesis 
concerning  the  motion  of  the  target  for  both  tracking  situations. 
We  assume,  thereby,  that  the  target  course  in  the  barrier  co- 
ordinate system  is  the  following  (the  broken  line  in  Fig.  22 
represents  the  base  course,  the  solid  line  the  maneuver  course): 


FIG.  22.  TRANSITOR'S  COURSE  THROUGH  BARRIER 
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The  mathematical  description  of  the  problem  of  fitting 
a linear  function  to  a set  of  observable  data  is  to  postulate 
the  relation 


X = a-^  + a^t  + e 
Y = + e 


(30) 


where : 

X = the  random  variable  corresponding  to  the  X- 
coordinate  of  the  target; 

Y = the  random  variable  corresponding  to  the  Y- 
coordinate  of  the  target; 

t = non-random,  observable  variable  time; 

e = random  error  of  mean  zero. 


a-p  a^,  a^,  = unknown  parameters 

Vf 

It  is  shown  that  the  minimum- variance , unbiased  estimators  of 
the  unknown  parameters  (a-p  a. 3 , and  can  be  obtained 
from  a set  of  data  i (Xp  Yp  t^)}^  by  the  following  relations: 


where 


n 

^ (X.  - X) (t . - t) 

— • JL  1 

°2  = 7 

I (tl  - c> 

i=l 


ii  = X - a^t, 

n n 


i-1  i=l 


(31) 


Mood,  A.  M.  and  3 raybill , F.  A.,  Introduction  to  the  Theory  of 
Statistics  , McGraw-Hill  Book  Company , Inc . , New  York , 1^63 . 
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L (Y.  - Y) (t.  - t) 
1=1 


l (t.  ' t) 


where 


Y - V* 

n <‘n 

h l 4^  = ^ Yi 

i=l  i=l 


It  is  clear  that  the  assumption  that  e has  a mean 
zero  and  is  normally  distributed  implies  that  X has  a mean 
(a^  + e^t)  and  a standard  deviation  about  that  mean  of  a . The 
minimum-variance,  unbiased  estimate  of  a is  shown  in  the  same 
reference  to  be 

n 2 % 

% * Ls^2  Z (X1  - ai  - “2  'd  ] • (33) 

i=l 

Similarly,  Y has  a mean  of  (a^  + cr^  t)  and  a standard 
deviation  a which  is  found  to  be  estimated  by 

y 

n 2 % 

°y  ■ rnr2  Z <Yi  - a3  • a4  tl)  ] • (34> 

n-2 

These  estimates  of  a provide  a measure  of  the  "goodness  of  fit" 
of  the  lines  to  the  data.  It  also  should  be  noted  that  the 
physical  interpretation  of  the  a's  is  that 

= best  estimated  X-position  of  the  target  at 
the  start  of  the  track  (t^  = 0) ; 

a.2  = best  estimated  X-velocity  of  the  target  during 
the  tracking  interval; 

a ^ = best  estimated  Y-position  of  the  target  at 
the  start  of  the  track  (t^  = 0) ; 
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= best  estimated  Y-velocity  of  the  target 
during  the  tracking  interval. 

These  parameter  estimates  can  then  be  used  by  the  TACTIC 
routine  to  determine  the  estimated  course  and  speed  of  the 
target  in  the  barrier  coordinate  system. 

The  tracking  routines  are  set  up  to  provide  current 
estimates  of  target  position  and  speed.  These  estimates  are 
based  on  all  information  stored  in  the  track  and  are  found 
using  the  following  equations: 

X (x ^ + a-2  t , 

Y = a3  + t,  (35) 

V = (a  2 2 + a42)\ 

where 

X is  the  current  estimated  X-coordinate  of  the  target, 

Y is  the  current  estimated  Y-coordinate  of  the  target, 

V is  the  current  estimated  speed  of  the  target, 

t is  the  time  since  the  track  was  initialized. 

4.4.6  Computation  of  Ambush  Parameters 

To  achieve  an  ambush,  the  patrol  submarine  must 
choose  a position  in  the  barrier  at  which  he  will  attempt  to 
attack  the  transitor.  This  choice  depends  on  the  following 
constraints : 

1.  He  must  reach  the  ambush  position  before  the 
transitor  enters  his  detection  region; 

2.  He  must  not  leave  the  barrier  region. 

The  selection  of  the  ambush  course  and  speed  depends 
on  the  current  position,  course,  and  speed  of  the  transitor,  and 
upon  the  current  position  of  the  patrol  submarine.  If  the  tran- 
sitor is  assumed  to  be  a diesel -electric  submarine,  the  ambush 
is  attempted  at  the  position  where  it  is  believed  that  he  will 
next  snorkel. 
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The  diagram  (Fig.  23)  of  the  problem  establishes  the 
technique  used  in  this  program  to  determine  an  ambush  position. 

It  should  be  noted  that  the  technique  used  here  is  not  necessarily 
that  used  by  the  fleet  nor  is  it  the  best  implementation  of  the 
ambush  doctrine. 
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In  the  diagram,  the  patrol  submarine  (Barrier  Ship) 
is  labeled  with  1,  the  transitor  (Penetrator)  with  2 , p is  the 
range  vector  between  the  ships,  is  the  velocity  vector  of 
the  transitor,  and  EDR  is  the  detection  range  of  the  patrol 
submarine . 

Clearly,  the  problem  is  to  find  the  pair  (VpT) 
which  satisfies  the  criteria,  and  once  (V^,T)  is  found,  to 
compute  the  course  8 . 

If  the  transitor  is  assumed  to  be  a diesel-electric, 
an  estimated  time  until  next  snorkel  is  provided,  namely  ETS, 
and  ambush  is  attempted  at  a specified  velocity  Vm^n  using 
that  time.  The  following  inequality  must  be  satisfied  if  ambush 
is  possible  at  this  time  and  velocity: 

Vmin (ETS)  s Ip  + (V2  ETS  + EDR) 1 . (36) 

If  ambush  is  impossible,  or  if  the  predicted  ambush  point  is 
outside  the  barrier,  the  velocity  is  increased  by  a fixed  amount 
to  V-p  and  ambush  is  attempted  at  the  original  ETS. 

The  process  continues  until  either  a pair  (VpETS) 
satisfies  the  ambush  criteria,  or  until  equals  the  maximum 
velocity  of  the  ship  and  the  ambush  criteria  are  still  not 
satisfied.  In  the  latter  case,  a flag  is  set  to  indicate  no 
ambush . 

If  the  equality  is  satisfied,  the  law  of  cosines  is 
used  to  find  the  course  angle  8 : 

2 2 ^2 

(V-.  T)  + |p|  - (v9t  + EDR) 

cos  B = — • U/; 

2(vxt)  Ip! 

If  the  transitor  is  assumed  to  be  a nuclear  submarine, 
no  time  until  next  snorkel  can  be  chosen,  and  some  initial  esti- 
mate Tq  of  a time  to  intercept  must  be  determined.  This  is 
accomplished  in  the  following  manner: 
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a = tt  - cos_1[p  . V2/(|p|  |V2I)] 

= |p|  cos  a , (38) 

^1 

To  = ’ where  Vx  is  initially  Vmin, 

and  the  same  inequality  is  again  tested: 

dx  * Ip  + (V2T0  + EWR)  I . (39) 

If  the  inequality  is  not  satisfied,  the  estimate 
Tq  is  increased  (by  10  min)  and  the  coordinates  of  the  pre- 
dicted ambush  point  are  computed,  etc. 

The  ambush  parameters  (VpT)  are  finally  determined 
just  as  before,  and  the  angle  8 selected  as  in  the  previous 
discussion. 

If  no  ambush  can  occur  even  at  the  maximum  speed  of 
the  patrol  submarine,  the  program  prints  a message  and  sets  a 
flag  indicating  this  fact. 

4.4.7  Torpedo  Model 

To  estimate  the  torpedo  intercept  probability  the 
ORL  CHOICE  IIB  MK-48  torpedo  model  was  used.  Due  to  its  com- 
plexity a complete  description  cannot  be  given  here.  However, 
it  will  be  useful  to  outline  this  model  as  it  relates  to  the 
problem  under  consideration. 

The  CHOICE  IIB  model  is  essentially  a simulation  of 
the  events  which  occur  from  the  instant  of  torpedo  launch  to 
intercept,  applied  to  only  one  particular  mode  of  use;  namely, 
it  is  assumed  that  the  torpedo  is  fixed  in  the  lead  intercept 
mode.  As  designed,  the  model  takes  into  account  the  following 

*0KL  Report  - R.  Greenstone  and  J.  McCutcheon.  "CHOICE,  A 
Probabilistic  Model  of  an  Advanced  Torpedo (U)",  Tech.  Memo. 
File  No.  TM412. 2711-71,  May  1966,  Contract  NOW  65-0123-d. 
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factors : 

(a)  Target  evasion  upon  detection  of  running 
torpedoes,  including  reaction  time. 


(b)  The  errors  in  range  and  bearing  data  provided 
to  the  torpedo  from  the  submarine's  sonar  system  during  the 
wire  guided  phase. 


(c)  The  noise  characteristics  of  the  torpedo  and 

the  target. 


(d)  The  torpedo  sonar  systems  (active  and  passive). 

(e)  The  torpedo  search  and  runout  speeds  and  the 
search  pattern. 


target . 


(f)  A single  zig-zag  maneuver  on  the  part  of  the 


(g)  Command  guidance  capability.  (During  the  wire 
guided  phase  one  torpedo  guidance  change  is  allowed). 

(h)  Torpedo  fuel  limitations. 

(i)  Target  speed  and  maneuv  ring  capability. 

(j)  Wire  length. 

(k)  Target  "baffle"  region  effect. 


The  essentials  of  the  model  can  be  understood  by 
referring  to  Figure  24.  Thus,  torpedo  launch  occurs  based 
on  the  estimated  target  track  and  speed,  and  guidance  is  pro- 
vided based  upon  the  attack  submarine's  sonar  system,  where 
it  is  assumed  that  the  target  always  is  at  the  center  of  the 
position  indicated  by  the  sonar  for  these  purposes.  After 
some  elapsed  time  the  target  can  exercise  a single  turn  man- 
euver and  the  torpedo  can  be  appropriately  redirected.  Further, 
the  target  may  have  detected  running  torpedoes  and  hence  react 
as  illustrated. 


I 
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The  torpedo  heads  initially  to  a precomputed  point 
based  on  the  initially  predicted  target  intercept  position,  or 
to  a newly  computed  position  if  torpedo  redirect  occurs.  At 
that  point  the  torpedo  commences  its  search  using  either  an 
active  or  passive  sonar  system.  The  target  position  at  this 
point  in  time  is  taken  to  lie  in  a region  of  space  defined  by 
the  range  and  bearing  errors  associated  with  the  submarine  sonar, 
in  accordance  with  a normal  probability  distribution  function 
of  both  range  and  bearing.  The  consequent  intercept  probability, 
taking  into  account  this  position  distribution  and  the  other 
factors  noted  earlier,  is  then  calculated. 

The  CHOICE  IIB  model  has  been  slightly  altered  for 
the  purposes  of  this  program.  These  changes  include  the 
following: 

1.  All  signal  or  noise  calculations  not  relating 
exclusively  to  the  MK-48  torpedo  apparatus  are  handled  by  the 
signal-to-noise  routines  discussed  in  this  report. 

2.  The  calculation  of  propagation  loss,  ship  self- 
noise, ship  passive  and  active  signal,  and  signal-to-noise  ratio 
are  handled  in  routines  exterior  to  the  CHOICE  IIB  model. 

The  above  changes  produce  more  realistic  S/N  levels  than  did  the 
original  CHOICE  IIB  model,  and  hence  produce  more  realistic 
probability  of  kill  results. 

4.5  COMMAND  AND  CONTROL  (TACTIC) 

The  routine  TACTIC  is  entered  to  evaluate  the  status 
of  the  simulation  as  provided  by  the  SUBSYS  information  and  to 
make  decisions  required  by  this  information.  The  decisions  are 
made  in  accordance  with  preassigned  doctrinal  principles  for 
each  boat.  A detailed  description  of  the  doctrinal  principles 
and  their  implementation  in  the  simulation  follows. 
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4.5.1  Doctrinal  Principles 

The  doctrinal  principles  employed  in  this  computer 
simulation  were  outlined  in  discussions  with  CNO.  For  the 
patrol  submarine  the  doctrinal  principles  are: 

(1)  In  the  event  that  the  patrol  submarine  achieves 
a detection,  it  will,  subject  to  certain  constraints,  carry  out 
an  attack  on  the  enemy  submarine  in  accordance  with  the  follow- 
ing three  attack  principles,  which  are  selected  in  the  order 
presented : 

a.  "Ambush"  - the  patrol  submarine  attempts 
to  reach  a position  in  the  path  of  the 
target  submarine  in  such  a way  that  torpedo 
launch  can  be  accomplished  from  a range  of 
choice  without  being  detected. 

b.  "Tail  Chase"  - the  patrol  submarine  attempts 
to  reach  a position  in  the  "baffle"  region 
of  the  target  submarine  in  such  a way  that 
torpedo  launch  in  this  region  can  be 
accomplished  from  a range  of  choice  without 
being  detected. 

c.  "Close  and  Attack"  - the  patrol  submarine 
elects  an  approach  tactic  which  enables 
torpedo  launch  to  be  accomplished  in  minimum 
time. 

(2)  The  constraints  placed  on  the  patrol  submarine 

are : 

a.  The  patrol  submarine  shall,  to  the  extent 
possible,  avoid  detection  by  the  enemy 
prior  to  torpedo  launch. 

b.  In  the  event  that  detection  by  the  enemy  is 
believed  to  have  occurred,  the  patrol  sub- 
marine shall  employ  evasive  tactics. 
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c.  The  patrol  submarine  will  not  leave  its 
forward  area  patrol  region  (the  barrier). 

The  transitor  doctrine  is  purely  evasive  in  nature. 

Thus,  if  the  transitor  detects  the  patrol  submarine,  evasive 
tactics  are  employed. 

This  tactical  doctrine  differs  somewhat  from  that 
used  in  TRACOR's  SS/SSN  Preliminary  Cost-Effectiveness  Analysis. 

The  analytical  nature  of  that  analysis  made  it  unnecessary  to 
consider  constraints  b)  and  c)  above. 

For  attack  principles  a)  and  b) , that  analytical 
study  used  "torpedo  launch  . . . from  a position  of  choice" 
rather  than  "torpedo  launch  . . . from  a range  of  choice.”  The 
analytical  study  assumed  perfect  tracking  information.  Although 
the  computer  simulation  addressed  the  tracking  problem,  there 
are  enough  uncertainties  in  the  accuracy  of  the  treatment  to 
make  unwarranted  an  attempt  to  maneuver  into  a position  of  choice 
before  weapon  launch. 

It  should  be  emphasized  that  the  detailed  tactics 
employed  in  the  computer  simulation  are  an  interpretation  of 
the  doctrinal  principles  outlined  by  CNO.  They  are  not  necessarily 
the  tactics  employed  by  U.  S.  submarines  nor  is  it  necessarily 
recommended  that  they  be  so  employed. 

4.5.2  Implementation  of  Tactical  Doctrine  in  the  Simulation 

4. 5. 2.1  Introduction.  The  tactics  employed  by  a submarine 

during  the  course  of  time  may  be  thought  of  as  the  implementa- 
tion of  the  submarine's  tactical  doctrine  subject  to  the  cumula- 
tive information  available  to  the  decision  maker  on  board  that 
submarine.  This  particular  implementation  does  not  necessarily 
represent  that  actually  used  by  the  fleet.  For  programming  pur- 
poses it  is  convenient  to  separate  the  tactics  of  each  submarine 
into  a number  of  parts.  Each  part  is  characterized  by  the  state 
of  the  cumulative  information  available  and  the  course  of  action 


91 

CONFIDENTIAL 


A. 


I 


required  by  that  submarine's  tactical  doctrine  given  that  infor- 
mation. Each  of  these  divisions  of  the  tactics  program  is  called 
a tactical  situation.  At  each  simulated  time  point  the  informa- 
tion available  to  each  ship  is  reevaluated  and  the  decision  is 
made  either  to  remain  in  the  current  tactical  situation  or  to 
transfer  to  another  tactical  situation. 

Nine  tactical  situations  are  required  for  the  patrol 
submarine.  Two  are  required  for  the  transitor.  Tables  IV  and  V 
contain  a brief  outline  of  these  tactical  situations,  and  de- 
tailed flow  charts  are  contained  in  the  Appendix.  A detailed 
description  of  each  tactical  situation  will  now  be  given. 

4. 5. 2.2  Tactical  situations  of  patrol  submarine.  Overriding 
all  of  the  tactical  considerations  of  the  patrol  submarine  is 
the  detection  of  an  active  sonar  ping.  In  the  event  that  a 
ping  is  detected  the  decision  maker  assumes  that  the  ping  is  a 
result  of  the  use  of  active  sonar  by  an  enemy  submarine.  He 
further  assumes  that  the  enemy  would  only  employ  active  sonar 
just  prior  to  a weapon  launch.  Therefore,  in  the  event  of  ping 
detection,  the  patrol  submarine  turns  to  a direction  opposite 
that  in  which  the  ping  was  received,  dives  to  maximum  depth  and 
runs  at  maximum  speed. 

4. 5. 2. 2.1  Tactical  Situation  No.  1 - No  Detection  Has 

Occurred.  This  is  the  tactical  situation  existing  prior  to  any 
detection.  The  patrol  submarine  is  completely  unaware  of  any 
intruder  which  might  be  in  its  vicinity.  The  patrol  submarine 
proceeds  on  its  preassigned  patrol  pattern  with  its  passive 
detection  system  in  operation.  The  existing  simulation  provides 
for  two  types  of  patrol  paths:  (1)  back  and  forth  along  a 

straight  line  parallel  to  the  X-axis  (See  Fig.  10)  and  (2)  a 
circular  path  with  center  at  any  point  and  with  an  arbitrary 
radius.  On  both  of  these  paths  the  submarine  may  proceed  at 
any  fixed  depth  and  move  with  any  fixed  speed.  The  turn  at 
the  end  of  each  patrol  segment  in  the  straight  line  case  is 
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TABLE  V.  SUMMARY  OF  TACTICS  OF  TRANS ITOR  SUBMARINE 


TACTICAL  SITUATION 

1 

2 

1.  Information  State 

No  Detection 

Detection  of 
another  ship 
has  occurred. 
Contact  may 
or  may  not  be 
maintained. 

2.  Objective 

1 

Transit  barrier 
Region 

Evade 

3.  Course  of  Action 

Transit  Region 

Steer  in  a dir- 
ection opposite 
t.  that  of  con- 
tact, run  at 
maximum  depth 
and  maximum  speed. 

4.  Occurrences  which  may 
require  alteration  or 
termination  of  course 
of  action  prior  to 
reaching  objectives 
(See  Table  C-VI) 

1 

94 


CONFIDENTIAL 


' 

J 


TABLE  VI.  OCCURRENCES  REQUIRING  ALTERATION  OR  TERMIN- 
ATION OF  COURSE  OF  ACTION  PRIOR  TO  REACHING  OBJECTIVE 


OCCURRENCE  ACTION 


1.  Submarine  is  required 
to  snorkel. 

1.  Submarine  snorkels 

2.  Patrol  submarine  reaches 
boundary  of  barrier 

2a.  Submarine  takes  new 
course  so  as  to  not 
cross  boundary  while 
proceeding  to  objec- 
tive or, 

2b.  submarine  fires 
torpedo  at  target 

3.  Contact  is  lost 

3a.  Submarine  maneuvers 
in  an  effort  to  re- 
gain contact  1 

3b.  Submarine  fires  tor-  1 
pedo  at  target  if 
fire  control  solution; 
is  available  or, 

3c.  Evade  1 

i 4.  Patrol  submarine  is 

estimated  to  be  with- 
in enemy's  assumed 
detection  range  and 
weapon's  range. 

4a.  Fire  torpedo  if  a 

fire  control  solution! 
is  available  or, 

4b.  Evade 

5.  Contact  makes  an 
unexpected  course 
change 

5.  Fire  torpedo  if  a 

fire  control  solution 
is  available  or, 
patrol  if  range  is 
closing  or,  close 
and  attack  if  range 
is  opening. 
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accomplished  by  rotating  the  heading  of  the  submarine  through 
180°.  All  other  turn  effects  are  neglected. 

If  the  patrol  submarine  is  a diesel-electric  it 
snorkels  in  accordance  with  some  preassigned  snorkeling  scheme. 
In  the  current  simulation  the  submarine  snorkels  with  some  pre- 
assigned frequency  for  some  preassigned  duration.  For  example 
snorkeling  begins  at  times  tg  + nT-^  and  ends  at  times  tg  + nT^  + 
T9 ; tg  is  a reference  time,  T-^  is  the  snorkeling  period,  is 

the  snorkeling  duration  and  n = 1,2,3 More  elaborate 

snorkeling  schemes  can  easily  be  incorporated  into  the  simula- 
tion. 

Tactical  Situation  No.  1 is  terminated  upon  the 
patrol  submarine's  detection  of  an  intruder,  and  the  patrol 
submarine  enters  Tactical  Situation  No.  2. 

4 . 5 . 2 . 2 . 2 Tactical  Situation  No.  2 - Detection  has  occurred, 
gather  information  for  tactical  decision.  In  this  Tactical 
Situation,  the  objective  of  the  patrol  submarine  is  to  gather 
information  sufficient  to  allow  him  to  localize,  classify,  and 
select  the  appropriate  mode  to  attack  the  intruder.  If  an 
accurate  classification  cannot  be  made,  the  decision  maker  in 
the  patrol  submarine  will  assume  that  the  intruder  is  an  enemy 
transitor.  Thus  the  paramount  objective  of  this  Tactical 
Situation  is  to  track  the  intruder.  In  this  Tactical  Situation 
the  patrol  submarine  maneuvers  on  those  courses  and  at  those 
speeds  necessary  to  obtain  the  intruder's  track,  and  continues 
to  do  this  either  until  a track  is  obtained  and  a decision  can 
be  made,  or  until  circumstances  force  the  discontinuance  of  this 
action.  The  actions  of  the  patroller  are  at  all  times  subject 
to  the  constraints  set  forth  in  the  doctrinal  principles. 
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Initial  Information 

Upon  entering  this  Tactical  Situation  the  sonar 
systems  of  the  patrol  submarine  provide  only  two  pieces  of 
information:  (1)  a detection  has  occurred  and  (2)  the  approxi- 

mate bearing  of  the  intruder.  Both  of  these  pieces  of  informa- 
tion are  provided  by  the  detection  subsystem.  In  this  simula- 
tion, however,  the  decision  maker  in  the  patrol  submarine  is 
provided  with  additional  information.  This  additional  informa- 
tion is  an  estimate  of  both  his  and  the  enemy's  detection  range, 
an  estimate  of  the  enemy's  speed,  and  an  estimate  of  the  enemy's 
weapon  range.  These  estimators  are  input  parameters  to  the 
simulation  and  are  an  approximation  of  the  experience  and  intelli- 
gence of  the  decision  maker.  At  the  time  of  detection  the  de- 
cision maker  in  the  patrol  submarine  should  be  able  to  make  at 
least  a gross  estimate  of  the  range  of  the  contact.  If  the  de- 
cision maker  is  assumed  to  have  recent  bathymetric  data  available, 
then  he  might  be  able  to  make  a rather  accurate  estimate  of  the 
range  of  the  contact.  Thus,  the  ability  of  the  decision  maker 
to  estimate  the  range  to  the  contact  is  represented  by  the  accu- 
racy of  that  input  parameter. 

Initial  Course  and  Speed  of  Patrol  Submarine 

Since  only  a bearing  is  available  at  the  time  that 
Tactical  Situation  No.  2 is  entered,  the  decision  maker  in  the 
patrol  submarine  assumes  that  the  contact  is  an  enemy  submarine 
attempting  to  transit  the  barrier.  He  further  assumes  that  if 
the  contact  is  in  the  northward  direction  then  the  transit  is 
from  north  to  south,  and  conversely  if  the  contact  is  in  the 
southward  direction.  That  is,  it  is  assumed  that  the  supposed 
transitor  has  not  passed  the  patrol  submarine  in  the  barrier. 

The  patrol  submarine  in  the  current  simulation  is  equipped  with 
the  PUFFS  (AN/BQG-4)  passive  ranging  subsystem.  This  subsystem 
cannot  provide  accurate  range  information  unless  the  target 
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bearing  is  within  about  40°  either  side  of  the  beam  of  the  ship 
upon  which  it  is  mounted.  Therefore,  upon  entering  Tactical 
Situation  No.  2 the  patrol  submarine  takes  a course  such  that 
the  contact  bears  either  70°  or  290°  relative  in  order  to  close 
on  the  contact  at  a moderate  rate  and  to  bring  PUFFS  to  bear 
upon  the  contact.  The  speed  of  the  patrol  submarine  during 
this  phase  of  the  tactics  may  be  any  preselected  value. 


During  this  Tactical  Situation  the  estimate  of  the 
range  to  the  contact  is  made  continuously.  The  methods  of  ob- 
taining this  estimate  are,  in  order  of  preference: 


(1)  Computation  using  the  last  range  provided  by 
PUFFS  and  the  estimate  of  the  intruder's  speed. 

(2)  If  contact  has  been  maintained  for  45  minutes, 
the  estimated  range  is  assumed  to  be  707,  of  the  exact  range.  A 
range  estimate  of  this  accuracy  should  be  obtainable  with  a 
bearings  only  solution  in  45  minutes. 

(3)  Computation  using  the  estimated  detection  range 
of  the  patrol  submarine  and  the  estimated  speed  of  intruder. 

The  estimated  range  to  the  contact  is  used  frequently  in  this 
Tactical  Situation. 

The  classification  subsystem  (AN/BQH-2C)  is  con- 
tinuously interrogated  during  the  course  of  this  Tactical 
Situation.  If  this  subsystem  fails  to  provide  information  the 
intruder  is  assumed  to  be  an  enemy  transitor. 

Course  and  Speed  Alterations 

The  patrol  submarine  continues  on  the  initial  course 
at  the  initial  speed  and  attempts  to  compute  the  track  of  the 
other  submarine  until  one  of  the  following  events  occurs: 
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(1)  The  patrol  submarine  reaches  an  edge  of  the 
barrier.  If  this  occurs  the  patrol  submarine  changes  its  heading 
so  as  to  remain  in  the  barrier  region  and  to  continue  to  bring 
PUFFS  to  bear  on  the  contact. 

(2)  Contact  is  lost.  If  contact  is  lost  and  the 
estimate  of  intruder's  range  is  less  than  the  estimate  of  the 
intruder's  detection  range  and  weapon  range,  the  decision  maker 
in  the  patrol  submarine  assumed  that  the  intruder  is  making  an 
attack  approach  on  the  patrol  submarine  and  by  taking  a course  180c 
from  the  last  bearing  of  the  intruder,  going  to  maximum  speed, 

and  diving  to  maximum  depth. 

If  the  estimated  range  to  the  contact  is  not 
less  than  the  contact's  estimated  weapon  range  and  if  some 
tracking  information  is  available  then  the  patroller  will: 

(a)  Chase  on  the  last  bearing  available  if  the 
range  is  opening  (Tactical  Situation  No.  4), 

(b)  Return  to  patrol  with  modified  patrol  limits 
if  the  range  is  closing  and  the  contact  is 
believed  to  have  been  in  the  convergence 
zone  (Tactical  Situation  No.  3), 

(c)  Attempt  to  compute  an  ambush  course  and 
speed  if  the  range  is  closing.  If  ambush 
is  possible  the  patrol  submarine  proceeds 
to  the  ambush  point  (Tactical  Situation 
No.  6).  If  it  is  impossible  to  compute 
an  ambush  course  and  speed  the  patrol 
submarine  returns  to  patrol  with  modified 
patrol  limits  (Tactical  Situation  No.  4). 

If  the  estimated  range  to  the  contact  is  not  less 
than  the  contact's  estimated  weapon  range  and  no  tracking  infor- 
mation is  available,  then  the  patrol  submarine  returns  to  patrol 
with  modified  patrol  limits. 
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(3)  The  estimated  range  to  the  contact  is  less  than 
the  contact's  estimated  detection  range  and  estimated  weapon 
range.  If  this  occvirs  the  patrol  submarine  does  one  of  the 
following  three  things: 

(a)  Fire  a torpedo  if  a fire  control  solution 
is  available  (Tactical  Situation  No.  5), 

(b)  Close  and  attack  if  intruder  is  evading 
(Tactical  Situation  No.  9), 

(c)  Evade  by  turning  to  the  direction  opposite 
to  the  direction  of  the  contact,  diving 

to  maximum  depth,  and  running  at  maximum 
speed  (Tactical  Situation  No.  8). 

(4)  The  relative  bearing  of  the  contact  is  no 
longer  within  30°  either  side  of  the  beam.  If  this  occurs  the 
patrol  submarine  comes  to  a new  course  such  that  the  relative 
bearing  of  the  contact  is  either  70°  or  290°. 

When  and  if  the  approximate  track  of  the  intru- 
der is  established  the  patrol  submarines  does  one  of  the  following 

(a)  If  the  range  is  opening  the  patrol  sub- 
marine initiates  a Tail-Chase  Attack 
(Tactical  Situation  No.  7). 

(b)  If  the  range  is  closing  an  attempt  is 
made  to  calculate  an  ambush  course  and 
speed.  If  an  ambush  is  possible  the 
patrol  submarine  proceeds  toward  the 
ambush  point  (Tactical  Situation  No.  6). 

(c)  If  ambush  is  not  possible  then  the  patrol 
submarine  proceeds  to  close  and  attack 
the  intruder  (Tactical  Situation  No.  9). 
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submarine  the  electrical  energy  remaining  in  the  batteries  is 
constantly  monitored.  If  the  remaining  energy  falls  below  a 
preassigned  fraction  of  the  maximum  possible  energy  of  the 
batteries  then  the  patrol  submarine  will  execute  a quiet  snorkel 
for  preassigned  duration.  In  the  present  simulation  the  fraction 
of  the  remaining  electrical  energy  is  maintained  at  1.0  so  that 
the  patrol  submarine  never  snorkels  in  this  Tactical  Situation. 
The  snorkeling  provision  is  included  at  present  to  allow  for 
later  improvements  in  the  program. 

4. 5. 2. 2. 3 Tactical  Situation  No.  3 - Contact  lost,  patrol.  This 
Tactical  Situation  occurs  when  contact  has  been  lost,  and  when 
either:  (1)  no  tracking  information  is  available  to  the  patrol 
submarine,  or  (2)  the  detection  is  believed  to  have  occurred  in 
the  convergence  zone.  The  maneuvers  of  the  patrol  submarine 
are  similar  to  those  of  Tactical  Situation  No.  1 in  that  the 
patrol  submarine  patrols  back  and  forth  along  a straight  line 
segment  parallel  to  the  X-axis  (see  Fig.  10).  The  limits  of 
the  patrol  path  are  not  the  limits  used  in  Tactical  Situation 
No.  1,  however.  It  is  assumed  that  the  contact  is  proceeding 
through  the  barrier  on  a course  parallel  to  the  Y-axis  of  the 
barrier  and  in  a direction  approaching  the  patrol  path  of  the 
patrol  submarine.  Thus,  the  patrol  limits  are  calculated  using 
the  last  contact  bearing,  the  last  range  to  contact  estimate 
(see  Tactical  Situation  No.  2 for  the  method  of  estimating 
range  to  contact),  and  a distance  which  specifies  the  patrol 
distance  either  side  of  the  contact's  assumed  path.  The  latter 
distance  is  arbitrary  and  is  an  input  to  the  program.  The 
Y-coordinate  of  the  patrol  path  is  that  value  of  Y reached  by 
the  patrol  submarine  at  the  time  of  lost  contact  plus  the  time 
delay  associated  with  entry  into  Tactical  Situation  No.  3. 
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When  and  if  contact  is  regained  the  patrol  submarine 
returns  to  Tactical  Situation  No.  2 and  attempts  to  obtain  the 
track  of  the  intruder. 

If  the  patrol  submarine  is  a diesel-electric  and  if 
the  energy  remaining  in  the  batteries  falls  below  a preassigned 
fraction  of  the  maximum  battery  energy,  then  the  submarine  will 
execute  a quiet  snorkel  for  a preassigned  period. 

4. 5. 2. 2. 4 Tactical  Situation  No.  4 - Contact  lost,  pursue.  This 
Tactical  Situation  occurs  when  contact  has  been  lost  and  when 
circumstances  make  it  desirable  to  move  in  the  direction  of  last 
contact.  Tactical  Situation  No.  2 would  occur,  for  example, 
when  contact  was  lost  because  of  evasive  tactics  of  the  enemy. 

In  this  Tactical  Situation  the  patrol  submarine 
maintains  a course  in  the  direction  of  the  last  contact,  and 
maintains  a preassigned  speed  until  either:  (1)  the  patrol  sub- 
marine reaches  its  maneuvering  limits  within  the  barrier  or 
(2)  contact  occurs. 

If  the  patrol  submarine  reaches  its  maneuvering 
limits  inside  the  barrier,  it  takes  a new  course  parallel  to 
the  boundary  edge  and  in  the  general  direction  of  the  last 
contact . 

If  detection  occurs,  the  patrol  submarine  proceeds 
to  close  and  attack  the  contact,  the  assumption  being  that  the 
contact  is  an  enemy,  and  that  the  contact  is  attempting  to 
evade.  The  close  and  attack  course  is  computed  on  the  assump- 
tion that  the  range  to  the  contact  is  the  estimated  detection 
range  of  the  patrol  submarine,  provided  the  decision  maker  is 
the  patrol  submarine,  and  that  the  contact  is  proceeding  directly 
away  from  the  patrol  submarine.  The  close  or  attack  speed  is 
an  arbitrary  input  variable. 
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If  the  patrol  submarine  is  a diesel-electric  and  if 
the  energy  remaining  in  the  batteries  falls  below  a preassigned 
fraction  of  the  maximum  battery  energy  then  the  patrol  sub- 
marine will  execute  a quiet  snorkel  for  a preassigned  period. 

4. 5. 2. 2. 5 Tactical  Situation  No.  5 - Fire  torpedo.  This  is  the 
weapon  launch  Tactical  Situation.  It  is  designed  to  simulate 
weapon  launch  and  weapon  and  ship  control  during  the  weapon's 
run.  In  the  current  simulation  all  of  these  functions  are 
accounted  for  in  the  CHOICE  IIB  MK-48  Torpedo  Simulation.  Thus 
Tactical  Situation  No.  5 merely  initializes  and  calls  the  CHOICE 
Subroutine . 

4. 5. 2. 2. 6 Tactical  Situation  No.  6 - Ambush.  In  this  Tactical 
Situation  the  Ambush  mode  of  attack  is  carried  out.  The  Tactical 
Situation  is  divided  into  two  phases: 

a.  Travel  to  the  ambush  pivot  position 

b.  The  attempt  to  reestablish  contact  if  contact 
has  been  lost  during  the  passage  to  the  pivot 
position. 

The  actual  approach  to  the  target  is  accomplished  in  another 
Tactical  Situation  (No.  9). 

The  method  of  calculating  the  course,  speed,  and 
time  required  to  reach  the  ambush  pivot  position  has  been  des- 
cribed. Upon  entering  this  Tactical  Situation  the  patrol  sub- 
marine changes  course  and  speed  as  required  and  dives  to  a depth 
sufficient  to  place  the  submarine  below  the  thermal  layer,  if 
a layer  exists.  This  maneuver  reduces  the  chances  of  detection 
by  the  enemy  submarine. 

The  patrol  submarine  proceeds  on  the  ambush  course 
at  the  required  speed  until  the  pivot  position  is  reached  unless 
one  of  the  following  events  happens: 
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A.  The  estimated  range  to  the  contact  is  less  than 
the  target's  estimated  weapon  range  and  detection  range.  The 
method  of  estimating  the  range  to  the  target  is  the  same  as  used 
in  Tactical  Situation  No.  2.  If  this  occurs  the  patrol  sub- 
marine will : 

(1)  Fire  weapon  if  a fire  control  solution  is 
available.  (Tactical  Situation  No.  5) 

(2)  Evade  by  turning  to  direction  opposite  to 
that  of  the  target,  going  to  maximum  speed 
and  diving  to  maximum  depth.  (Tactical 
Situation  No.  8). 

B.  The  target  makes  a major  course  change  while 
contact  is  maintained.  If  this  occurs  the  patrol  submarine 
will: 

(1)  Fire  weapon  if  a fire  control  solution  is 
available.  (Tactical  Situation  No.  5). 

(2)  Close  and  attack  if  the  target  is  moving 
away  from  the  patrol  submarine.  (Tactical 
Situation  No.  9). 

(3)  Attempt  to  retrack  the  target  if  the  target 
is  moving  toward  the  patrol  submarine. 
(Tactical  Situation  No.  2). 

C.  Contact  has  been  regained  after  being  lost.  If 
this  occurs  the  patrol  submarine  will  attempt  to  retrack  the 
target.  (Tactical  Situation  No.  2). 

If,  upon  arrival  at  the  ambush  pivot  position, 
contact  has  been  maintained,  the  patrol  submarine  immediately 
begins  to  close  and  attack  the  enemy  at  minimum  speed.  (Tacti- 
cal Situation  No.  9). 

If  contact  has  been  lost,  if  the  patrol  submarine  is 
a diesel "electric , and  if  the  energy  remaining  in  the  batteries 
is  below  a preassigned  fraction  of  the  maximum  battery  energy, 
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then  the  patrol  submarine  will  execute  a quiet  snorkel.  The 
duration  of  this  snorkel  is  calculated  so  that  the  snorkel 
period  will  end  at  the  approximate  time  that  useful  detection 
information  starts  to  be  received.  The  snorkeling  is  carried 
out  as  the  ship  moves  in  a circular  path  about  the  pivot  posi- 
tion at  minimum  speed. 

Upon  completion  of  the  snorkel  period,  or  upon 
arrival  at  the  pivot  position  if  no  snorkeling  is  required  and 
contact  has  been  lost,  the  patrol  submarine  turns  so  as  to  move 
along  the  previously  computed  track  of  the  enemy  at  a pre- 
assigned speed. 

The  patrol  submarine  continues  on  this  course  at 
this  speed  until  redetection  has  occurred  or  until  the  barrier 
boundary  is  reached.  If  the  patrol  submarine  reaches  the  barrier 
boundary  he  takes  a new  course  so  as  to  move  parallel  to  that 
boundary  in  the  same  general  direction  as  the  previous  course. 

If  redetection  occurs,  he  immediately  proceeds  to  close  and 
attack  the  enemy  at  minimum  speeds.  (Tactical  Situation  No.  9). 

4. 5. 2. 2. 7 Tactical  Situation  No.  7 - Tail  chase.  In  this  Tacti- 
cal Situation,  the  Tail  chase  mode  of  attack  is  carried  out. 

The  course  taken  by  the  patrol  submarine  is  essentially  a bear- 
ing rider  course.  Instead  of  steering  toward  the  enemy,  however, 
the  patrol  submarine  steers  so  that  the  enemy  is  approximately 
60°  off  the  bow.  This  is  done  so  that  PUFFS  will  provide  ranging 
information. 

Initially  the  patrol  submarine  takes  a course  so 
that  the  enemy  bears  ±60°  off  the  bow  and  moves  at  a preassigned 
speed.  The  patrol  submarine  maintains  this  course  and  speed 
until  one  of  the  following  events  happens: 


105 

CONFIDENTIAL 


CONFIDENTIAL 


LANE 


AUSTIN  TEXAS  78721 


A.  Contact  is  lost.  If  this  occurs  the  patrol  sub- 
marine steers  in  the  direction  of  last  contact  at  a preassigned 
speed  in  an  effort  to  regain  contact.  (Tactical  Situation  No.  4). 

B.  The  patrol  submarine  reaches  the  boundary  of 
the  barrier.  If  this  occurs,  the  patrol  submarine  fires  a 
torpedo  regardless  of  the  quality  of  the  fire  control  solution. 

C.  The  enemy's  bearing  is  less  than  55°  or  more 
than  65°  either  side  of  the  bow.  If  this  occurs  the  patrol 
submarine  adjusts  course  so  that  the  enemy  bears  60°  off  the 
bow. 

D.  Tracking  information  indicates  that  the  range  is 
opening.  If  this  occurs,  the  speed  of  the  patrol  submarine  is 
increased,,  If  it  is  impossible  to  increase  speed  the  patrol 
submarine  fires  a torpedo  regardless  of  the  quality  of  the  fire 
control  solution. 

E.  The  estimated  range  to  the  enemy  is  less  than  the 
enemy's  estimated  weapon  range  and  estimated  detection  range. 

If  this  occurs,  the  patrol  submarine  fires  a torpedo  if  a fire 
control  solution  is  available.  If  a fire  control  solution  is 
not  available,  the  patrol  submarine  evades  by  turning  to  a 
direction  opposite  to  that  of  the  enemy,  diving  to  maximum 
depth,  and  running  at  maximum  speed. 

F.  A fire  control  solution  is  available  and  the 
estimated  range  to  the  target  is  less  than  the  preassigned 
maximum  firing  range.  When  this  occurs  the  patrol  submarine 
fires  a torpedo  at  the  enemy. 

The  termination  by  occurence  F is  the  ultimate 
goal  of  this  Tactical  Situation.  It  should  be  noted  that  the 
torpedo  launch  may  or  may  not  be  made  when  the  patrol  submarine 
is  in  the  baffle  region  of  the  enemy.  Uncertainties  in  the 
accuracy  of  the  treatment  of  the  tracking  problem  make  unwarranted 
any  attempt  to  maneuver  into  the  baffle  region  of  the  enemy. 
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4. 5. 2. 2. 8 Tactical  Situation  No.  8 - Evade.  This  Tactical 
Situation  occurs  when  the  decision  has  been  made  to  evade  the 
other  submarine.  Upon  entering  this  Tactical  Situation  the 
patrol  submarine  turns  to  a course  180°  away  from  the  contact, 
dives  to  maximum  depth,  and  runs  at  maximum  speed.  The  patrol 
submarine  will  continue  on  this  course  at  this  speed  and  depth 
unless  one  of  the  following  events  happens: 

A.  The  patrol  submarine  reaches  the  boundary  of 

the  barrier.  If  this  occurs  the  patrol  submarine  takes  a course 
parallel  to  the  boundary  in  a direction  away  from  the  last 
bearing  of  the  enemy. 

B.  The  patrol  submarine  obtains  the  bearing  of  the 
enemy.  He  will  then  change  course  so  that  the  enemy  bears  180° 
relative.  Ordinarily,  the  enemy  will  be  in  the  patrol  submarine' 
blind  spot. 

4. 5. 2. 2. 9 Tactical  Situation  No.  9 - Close  and  attack.  This 
Tactical  Situation  occurs  when  the  decision  has  been  made  to 
close  the  target  and  fire  as  soon  as  a fire  control  solution  is 
available  and  the  estimated  range  to  the  target  is  less  than  the 
maximum  firing  range.  Initially  the  course  taken  by  the  patrol 
submarine  is  a modified  lead  intercept  course.  That  is,  the 
patrol  submarine  takes  a lead  intercept  course  if  this  course 
will  allow  PUFFS  to  bear  on  the  enemy.  If  PUFFS  will  not  bear 
upon  the  enemy  the  course  is  adjusted  so  that  PUFFS  will  bear. 

The  speed  of  the  patrol  may  be  any  preassigned  number  greater 
than  the  speed  required  to  close  the  lead  intercept  triangle. 

The  patrol  submarine  continues  on  this  course  at 
this  speed  unless  one  of  the  following  events  happens: 
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A.  The  patrol  submarine  reaches  the  boundary  of 

the  barrier.  If  this  occurs  the  patrol  submarine  fires  a torpedo 
regardless  of  the  quality  of  the  fire  control  solution. 

B.  Contact  is  lost.  If  this  occurs,  the  patrol 
submarine  steers  in  the  direction  of  last  contact  at  a pre- 
assigned speed  in  an  effort  to  regain  contact.  (Tactical 
Situation  No.  4). 

C.  The  target  does  not  bear  between  55°  and  65° 
either  side  of  the  bow.  If  this  occurs  a new  modified  lead 
intercept  course  and  speed  is  computed  and  the  patrol  submarine 
takes  this  new  course  and  speed. 

D.  The  estimated  range  to  the  target  is  less  than 
the  estimate  of  the  target's  weapon  range  and  detection  range. 

If  this  occurs,  the  patrol  submarine  fires  at  the  target  if  a 
fire  control  solution  is  available.  If  no  fire  control  solution 
is  available,  the  patrol  submarine  evades  by  turning  to  a dir- 
ection opposite  to  that  of  the  enemy,  diving  to  maximum  depth 
and  running  at  maximum  speed.  (Tactical  Situation  No.  8). 

E.  A satisfactory  fire  control  solution  exists  and 
the  estimated  range  to  the  target  is  less  than  the  maximum 
firing  range,  '//hen  this  occurs,  the  patrol  submarine  fires  a 
torpedo  at  the  target. 

4. 5. 2. 3 Tactical  Situations  of  the  Transitor.  Overriding  all 
of  the  tactical  considerations  of  the  transitor  is  the  detection 
of  an  active  sonar  ping.  As  with  the  patrol  submarine,  if  a 
ping  is  detected,  the  transitor  immediately  evades  by  turning 
to  a direction  opposite  that  in  which  the  ping  was  received, 
diving  to  maximum  depth  and  running  at  maximum  speed. 

4 . 5 . 2 . 3 . 1 Tactical  Situation  No.  1 - No  Detection  has  occurred. 
This  is  the  tactical  situation  existing  prior  to  any  detection. 
While  in  this  tactical  situation  the  transitor  proceeds  to  pass 
through  the  barrier. 
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It  is  not  assumed  that  the  transitor  is  aware  ot  the 
precise  location  of  the  barrier  but  it  is  assumed  that  the  tran- 
sitor is  aware  that  he  is  in  a region  of  possible  enemy  activity. 
Therefore,  the  transitor  carries  out  precautionary  tactics  by 
zig-zagging.  This  is  accomplished  by  defining  a base  course, 

cd  , a maximum  deviation  from  base  course  AT,  and  a turn  time 

o 

increment.  At.  At  times  tT  + At,  where  tT  is  the  time  of  last 
turn,  a new  course  between  cpQ  - Acp  and  cpQ  + Am  is  selected  at 
random.  The  speed  of  the  transitor  is  a preselected  constant. 

If  the  transitor  is  a diesel -electric , he  snorkels 
in  accordance  with  some  preassigned  snorkeling  scheme.  In  the 
current  simulation  the  transitor  snorkels  with  some  preassigned 
frequency  for  some  preassigned  duration.  The  snorkeling  scheme 
is  the  same  as  that  employed  by  the  patrol  submarine  in  Tactical 
Situation  No.  1. 

When  the  transitor  detects  another  ship,  he  assumes 
that  the  contact  is  an  enemy  and  takes  immediate  evasive  action 
by  turning  to  the  direction  opposite  to  that  of  the  contact, 
diving  to  maximum  depth  and  running  at  maximum  speed. 

4. 5. 2. 3. 2 Tactical  Situation  No.  2 - Evade.  This  Tactical 
Situation  occurs  when  the  decision  has  been  made  to  evade  the 
contact.  The  patrol  submarine  continues  on  the  evade  course 
selected  in  Tactical  Situation  No.  1,  at  maximum  depth  and 
maintains  maximum  speed  unless  the  transitor  obtains  the  bear- 
ing of  the  enemy.  If  this  occurs,  the  transitor  will  change 
course  so  that  the  contact  bears  180°  relative.  Ordinarily 
the  enemy  will  be  in  the  transitor' s blind  spot. 
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5 . SAMPLE  RESULTS 

The  performance  description  of  the  subsystems  employed 
in  the  Forward  Area  Patrol  Mission  is  statistical  in  nature.  The 
statistical  distribution  of  the  performance  of  each  of  these  sub- 
systems has  been  accounted  for  in  the  simulation  at  the  subsystem 
level,  by  random  sampling  from  the  appropriate  distribution.  To 
demonstrate  the  effect  of  these  distributions  upon  decisions  and 
upon  the  ultimate  outcome  of  the  engagement,  a particular  geo- 
metry of  encounter  was  chosen.  The  participants  were  each  diesel- 
electric  submarines.  The  barrier  submarine  possesses  the  follow- 
ing subsystems  in  this  example:  Detection  - AN/BQR-2  (DIMUS); 

Localization  - AN/BQR-2  (DIMUS) (Bearing) , AN/BQG-4  (Range  and 
Bearing);  Classification  - AN/BQH-2C;  Weapon  - MK-48  Torpedo. 

The  transitor  submarine  possesses  the  following  subsystems  in 
this  example:  Detection  - AN/BQR-2B;  Localization  - AN/BQR-2B 

(Bearing).  By  this  choice  of  subsystems,  the  Barrier  Submarine 
is  given  a detection  superiority.  The  transitor  is  initially 
snorkeling,  on  a base  course  parallel  to  the  Y-axis  in  Figure  25. 
To  disguise  this  course  he  makes  turns  every  5 minutes.  These 
turns  are  randomly  chosen  between  ± 20°  of  his  base  course.  The 
model  was  exercised  three  times.  Each  time  a different  random 
sample  of  the  performance  of  each  subsystem  was  chosen. 

The  results  of  the  three  runs  are  somewhat  dissimilar. 

In  the  first  two  runs,  the  Barrier  Ship  detected  the  transitor, 
tracked  him  until  the  Barrier  Ship  felt  he  had  a good  estimate 
of  the  transitor 's  base  course,  computed  an  intercept  position 
based  on  the  estimated  length  of  time  before  the  transitor  would 
next  snorkel,  intercepted,  and  launched  a weapon.  The  kill 
probability  for  the  first  run  was  0.941;  the  kill  probability  for 
the  second  run  was  0.938.  Both  of  these  runs  are  displayed  in 
Figure  25.  In  the  third  run,  the  Barrier  Ship  failed  to  reacquire 
the  target  upon  reaching  the  intercept  position.  The  transitor 
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thus  escaped,  and  the  kill  probability  is  0.0,  as  no  weapon  was 
launched . 

Three  runs  are,  of  course,  insufficient  to  give  a 
complete  measure  of  the  effectiveness  of  the  subsystems  employed, 
but  the  effect  of  the  statistical  distribution  of  subsystem 
performance  is  clearly  indicated. 

It  is  now  necessary  to  mention  some  peculiar  aspects 
of  the  two  runs  illustrated  in  Fig.  25. 

First,  it  will  be  observed  that  the  Barrier  Submarine 
does  not  snorkel  throughout  the  encounter.  This  is  because  it 
is  assumed  that  his  reserves  are  sufficient  to  sustain  him 
without  snorkeling.  In  the  absence  of  a contact  the  barrier 
ship  would  snorkel  normally. 

Second,  it  is  noted  that  the  ambush  position  is  some 
distance  from  the  base  course  of  the  transitor.  This  testifies 
to  the  success  of  the  transitor  in  disguising  his  true  course. 

It  also  suggests  that  the  tactics  used  by  the  Barrier  Ship  to 
decipher  the  transitor 's  base  course  require  some  improvement, 
if  not  in  real  life,  then  at  least  in  the  model. 

The  last  position  for  each  ship  shown  in  the  diagram 
is  that  at  weapon  launch.  From  that  point,  the  torpedo  model 
Choice  IIB  completes  the  simulation,  but  does  not  now  permit 
any  recovery  of  the  participant's  positions  as  a function  of 
time . 

One  point  peculiar  to  both  runs  should  be  mentioned. 

The  Barrier  Ship  redetected  the  transitor  after  reaching  the 
intercept  position  and  attempted  to  close  with  the  transitor. 

Due  to  a discrepancy  in  the  current  model,  the  Barrier  Ship 
increased  his  speed  to  such  a point  that  the  ensuing  increase 
in  self-noise  prevented  his  obtaining  any  range  data  for  a 
short  time.  At  close  range,  and  prior  to  detection  by  the 
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transitor,  the  Barrier  Ship  was  able  to  obtain  range  data, 
establish  a fire  control  track,  and  launch  his  weapon. 
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6.  IMPROVEMENT  AND  EXTENSION  OF  PRESENT 
SIMULATION  MODEL 

The  two-submarine  simulation  model  described  in  this 
Technical  Memorandum  is  designed  to  be  both  a model  of  the 
Forward  Area  Patrol  Mission  and  the  initial  phase  of  a multi- 
participant ASW  Combat  Simulation  Model.  The  purpose  of  this 
chapter  is  to  discuss  means  by  which  the  current  model  can  be 
made  to  simulate  more  realistically  the  Forward  Area  Patrol 
Mission  and  to  point  out  the  problems  which  must  be  solved 
before  a multi-participant  model  can  be  constructed. 

The  present  model  can  more  realistically  simulate  the 
performance  of  each  participant  in  ASW  Combat  by  removing  the 
deficiencies  in  each  of  the  major  areas: 

1.  Subsystem  Performance, 

2.  Weapon  Performance, 

3.  Information  and  Noise, 

4.  Assessment,  Decision,  and  Action  Processes,  and 

5.  Vehicular  Motion. 

The  deficiencies  stem  partially  from  insufficient  information, 
and  partially  from  time  allowed  for  construction  of  the  current 
model . 

6 . 1 SUBSYSTEM  PERFORMANCE 

The  most  critical  deficiencies  in  the  simulation  are 
in  subsystem  and  weapon  performance.  This  is  true  from  the 
point  of  view  of  both  the  Forward  Area  Patrol  Mission  Simulation 
and  the  general  multi-participant  simulation. 

Serious  deficiencies  exist  in  all  categories  of  sonar 
subsystem  performance.  The  simulation  of  passive  sonar  detection 
of  a single  contact  for  the  dynamic  situation  is  believed  to  be 
rather  accurate.  However,  the  technique  lacks  experimental 
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validation.  At  present,  the  effect  of  multiple  contacts  on 
passive  detection  is  not  accurately  known.  The  same  is  true 
for  active  sonar  detection.  The  effect  of  mutual  interference 
on  active  sonar  detection  is  not  accurately  known. 

The  deficiencies  in  the  simulation  of  existing  local- 
ization subsystems  result  primarily  from  a paucity  of  experimental 
data  of  sufficient  quality.  Theoretical  error  analyses  exist 
for  many  localization  subsystems;  however,  these  analyses  do  not 
consider  all  of  the  factors  involved  in  subsystem  performance. 

For  example,  the  effect  of  own  ship  motion  is  almost  always 
neglected.  Furthermore,  these  analyses  consider  only  one 
contact  even  though  the  performance  of  the  localization  subsystem 
in  the  presence  of  multiple  contacts  and  mutual  interference  is 
a matter  of  extreme  importance.  Localization  subsystems  can 
be  accurately  simulated  only  after  more  and  better  performance 
data  has  been  obtained  and  analyzed. 

The  classification  problem  requires  a large  amount  of 
study  before  it  can  be  accurately  simulated. 

The  final  major  simulation  deficiencies  in  the  area 
of  subsystem  performance  involve  the  fire  control  and  tracking 
subsystems.  In  the  present  simulation,  the  tracking  problem 
is  treated  by  a least  squares  linear  fit  to  position  estimates 
of  the  contact  provided  by  the  range  and  bearing  subsystem  and 
by  assuming  perfect  knowledge  of  own  ship  position.  To  be 
realistic,  the  model  should  be  capable  of  accurately  simulating 
the  performance  of  automatic  trackers,  fire  control  computers, 
and  plotting  parties.  Time  limitations  prevented  an  attempt  to 
do  this  for  the  present  model.  If  adequate  experimental  data  is 
available,  these  simulations  will  be  included  in  the  general 
model . 
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6 .2  WEAPON  PERFORMANCE 

The  present  computer  model  relies  on  the  Choice  IIB 
torpedo  model  to  simulate  weapon  performance . This  torpedo 
model  was  developed  by  Ordnance  Research  Laboratory  at  Penn  State 
for  use  in  the  design  of  the  MK-48  torpedo.  Though  useful  as  a 
design  tool,  Choice  IIB  is  too  inflexible  to  be  satisfactory  as 
the  torpedo  simulation  for  the  computer  model  for  the  following 
reasons : 

1)  Choice  IIB  models  only  the  MK-48  torpedo  while 
in  reality  each  ship  may  possess  a variety  of  torpedoes. 

2)  Choice  IIB  permits  the  interaction  of  only  two 
ships  - one  a firing  ship,  the  other  a target  ship  - while  in 
reality,  each  ship  could  play  both  roles,  simultaneously. 

3)  Choice  IIB  permits  each  ship  very  limited  interro- 
gation of  the  model  once  the  weapon  has  been  launched,  thus 
making  decisions  on  unrealistic  information. 

4)  Choice  IIB  simulates  the  MK-48  in  conjunction 
with  the  PUFFS  sonar  system,  when  in  reality,  the  MK-48  or  any 
torpedo  can  be  used  in  conjunction  with  various  sonars. 

To  serve  the  ASW  Combat  Simulation  Model,  a torpedo 
simulation  must  be  developed  to  satisfy  the  following  require- 
ments : 

1)  The  simulation  must  be  of  a modular  design.  To 
model  different  torpedoes  should  merely  require  replacing 
component  subprograms. 

2)  The  simulation  must  permit  the  torpedo  to  become 
a participant  in  the  model  with  time-stepped  motion,  command 
and  control  function,  etc. 

3)  Each  ship  should  be  able  to  launch  a torpedo 
while  himself  under  attack. 

4)  The  simulation  must  give  each  ship  the  same 
ability  to  interrogate  the  model  both  before  and  after  launch  of 
weapon . 
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5)  The  simulation  must  be  compatible  with  any  weapon 
system  of  any  ship. 

Since  the  Forward  Area  Patrol  Simulation  model  required 
the  simulation  of  torpedo  performance  only,  the  problem  of 
simulating  the  performance  of  other  weapons  was  not  addressed. 

It  seems  reasonable  to  assume  that  the  performance  of  all 
required  weapons  can  be  adequately  simulated  provided  adequate 
knowledge  of  this  performance  is  available. 

6 . 3 THE  MEDIUM  AND  INFORMATION  FLOW 

Insufficient  target  strength  data  necessitated 
the  assumption  of  a 15  dB  target  strength.  More  and  better 
data  is  required  to  correct  this  deficiency. 

Inaccuracies  in  the  simulation  of  both  radiated  noise 
and  self-noise  also  resulted  from  insufficient  experimental 
data.  All  noise  levels  simulated  in  the  present  model  are 
based  on  experimental  data;  but,  in  general,  only  a few  data 
points  are  available.  Thus,  interpolation  inaccuracies  are 
introduced.  Furthermore,  the  experimental  data  is  such  that 
isotropic  noise  must  be  assumed.  In  fact,  noise  is  non- 
isotropic and  this  should  be  accounted  for  in  the  simulation. 

The  simulation  of  reverberation  must  be  included  in  future 
improvements  of  the  model. 

Propagation  loss  for  a particular  oceanographic 
condition  and  for  each  receiving  sensor  is  stored  as  a table  of 
source  depth,  receiver  depth,  and  horizontal  range  between 
source  and  receiver.  Currently,  only  three  combinations  of 
source  and  receiver  depths  should  be  included.  Propagation 
loss  curves  for  other  oceanographic  conditions  should  also  be 
included  if  oceanographic  conditions  are  expected  to  change 
during  the  course  of  the  combat. 
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In  reality,  propagation  loss  is  a function  of  time, 
due  to  the  rapid,  local  variation  of  oceanographic  microstruc- 
ture, as  well  as  to  the  longer  term,  general  variation  of 
oceanographic  macrostructure.  A technique  of  including 
propagation  loss  as  a function  of  time  should  be  included  in 
the  simulation. 

A more  accurate  technique  of  computing  propagation 
loss  would  increase  the  accuracy  of  the  simulation. 

6.4  ASSESSMENT,  DECISION,  AND  ACTION  PROCESSES 

To  better  simulate  the  assessment  and  decision 
processes  effected  by  the  decision  maker,  it  is  necessary  to 
determine  in  more  detail  what  criteria  are  used  by  the 
decision  maker  to  formulate  action.  A more  detailed  description 
of  tactical  doctrine,  and  a complete  idea  of  the  role  of 
knowledge  and  experience  are  both  necessary  to  a realistic 
simulation  of  the  command  and  control  function. 

In  addition  to  these  areas  of  human  judgment,  a 
knowledge  of  how  battery  level  is  a function  of  ship  activity, 
and  the  length  of  time  between  decision  and  action  for  each 
action  are  also  needed. 

6.5  VEHICULAR  MOTION 

The  current  model  permits  a participant  three 
possible  types  of  motion: 

1.  Straight  line,  constant  speed; 

2.  Circular  course,  constant  speed; 

3.  Straight  line,  constant  acceleration. 

There  are  three  dificiencies  in  this  description  of  motion: 

1.  Acceleration  is  constant. 

2.  The  equations  of  motion  are  integrated  over  a 
complete  time  step. 
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3.  Only  three  degrees  of  freedom  are  allowed  each 
participant . 

A time -dependent  acceleration  is  required  to 
accurately  simulate  both  linear  acceleration  and  turning  motion. 
If  the  appropriate  accelerations  were  available,  the  command  and 
control  simulation  for  each  participant  would  specify  the  new 
course  desired  and  rudder  angle,  the  new  depth  desired  and  rate 
of  descent,  or  the  new  speed  desired,  and  the  simulation  could 
be  programmed  to  accurately  simulate  the  resulting  motion. 

The  finite  time  increment  technique  requires  that  the 
change  in  position  of  each  participant  occuring  during  the  time 
step  be  evaluated  at  the  end  of  the  time  step.  There  is  no 
requirement  that  the  equations  of  motion  be  integrated  over  the 
entire  time  step.  Thus,  the  accuracy  of  motion  simulation  can 
be  increased  by  allowing  velocity  changes  to  occur  during  a 
time  step.  This  can  be  accomplished  without  difficulty.  It  was 
felt,  however,  that  the  absence  of  a time -dependent  acceleration 
made  such  a feature  superfluous. 

The  third  motion  deficiency  is  the  restriction  of 
each  participant  to  three  degrees  of  freedom.  This  restriction 
implies  that  each  participant  is  a point  mass  while  in  reality 
each  is  a rigid  body.  To  be  realistic  each  participant  should 
be  allowed  six  degrees  of  freedom.  This  would  permit  the 
simulation  of  roll,  pitch,  and  yaw  effects.  These  quantities 
affect  the  performance  of  sonar  sensors  and  own  ship  motion 
sensors . The  six  degrees  of  freedom  problem  can  be  handled  in 
the  simulation  if  all  forces  and  torques  acting  on  each 
participant  are  known. 
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APPENDIX 

FLOW  CHART  DESCRIPTION  OF  FORWARD  AREA 
PATROL  SIMULATION  MODEL 

The  purpose  of  this  appendix  is  to  provide  detailed 
information  on  each  routine  in  the  Forward  Area  Patrol  Simulation 
Model.  The  theory  upon  which  the  routines  are  based  is  contained 
in  the  body  of  the  report. 

This  appendix  is  conveniently  divided  into  four  parts, 
each  part  describing  a major  section  of  the  simulation  model.  The 
four  parts  are  CONTOL  (Control  and  I/O  Routine),  POSIT  (Position 
Routine),  SUBSYS  (Routines  associated  with  Subsystems),  and  TACTIC 
(Routines  associated  with  command  decisions).  The  format  for  each 
part  is  the  same.  Each  separate  routine  is  identified  by  a brief 
description  of  its  purpose.  A description  of  the  input,  output 
and  internal  variables  follows.  The  documentation  of  each  routine 
concludes  with  a flow  chart  of  the  routine. 

The  torpedo  model  used  in  the  Forward  Area  Patrol 

1 2 

Simulation  is  described  in  separate  documents  ’ , and  a description 
will  not  be  included  in  this  report. 


Greenstone,  R. , and  J.  W.  McCutcheon,  "Choice)  A Probabi' istic 
Model  of  an  Advanced  Torpedo,"  Tech.  Memo.  TM  412.2711-71,  18  May 
1966,  (Contract  N0w65-0123-d) . 

O 

Greenstone,  R. , and  J.  J.  Tully,  "The  Effect  of  Torpedo  Radiated 
Noise  on  the  MK48  MOD  0 Torpedo  Weapon  System,"  Part  5,  Tech.  Memo 
TM  412.2711-67,  6 May  1966,  (Contract  N0w65-0123-d) . 
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A . 1 CONTOL 

PURPOSE:  This  is  the  driver  routine  for  the  simulation  model. 

The  routines  POSIT,  SUBSYS,  and  TACTIC  are  called  in  that  order 
by  the  driver  routine  at  each  time  step. 

In  addition  to  being  the  driver  routine  CONTOL  is  the 
only  input  and  principal  output  routine.  The  input  and  output 
roles  of  CONTOL  will  be  considered  in  detail  in  this  section. 

Under  the  description  of  input,  each  card  is  labeled 
with  the  variables,  and  a definition  is  included  for  each  variable, 
A sample  ot  the  labeled  input  as  printed  by  the  program  is  in- 
cluded after  all  input  variables  have  been  defined. 

In  the  description  of  output,  each  column  heading  is 
defined,  and  a copy  of  sample  output  included  for  reference. 
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INPUT 


FORMATS: 


INTEGER 
FLOATING 
LABELS 
CARD  IMAGE 


(1615) 

(4E20.10) 

I 8(A6,4X)] 
(13A6.A2) 


READ  - (INTEGER)  - N BOAT , r ITYPE  (I),  1=1,  NBOAT] 


NBOAT  - Number  of  boats  participating  in  model=2 
ITYPE  - Defines  MISSION  for  each  boat 


0 = Barrier  Boat 

1 = Penetrator 


(3  numbers)  (1  card) 

READ  - (INTEGER)  - IPRNT1,  IPRNT2,  IPRNT3 , IPRNT4 , IPRNT 5 , 

IPRNT6,  IPRNT7 

IPRNT1  - Tag  for  printing  positions  of  both  ships 

IPRNT 2 - Tag  for  printing  subsystem  values  of  Barrier 
Ship 

IPRNT 3 - Tag  for  printing  subsystem  values  of  Penetrator 
Ship 

IPRNT4  - Tag  for  printing  tracking  information  of  Barrier 
Ship 

IPRNT 5 - Tag  for  printing  fire  control  information  of 
Barrier  Ship 

IPRNT6  - Tag  for  plotting  information  to  be  placed  on 
Fastrand 

IPRNT 7 - Tag  for  PRINT1  call  at  every  time  step 

Prints  all  pertinent  information  for  all  options: 

0 = no 

1 = yes 


A 4 

UNCLASSIFIED 


UNCLASSIFIED 


LANE 


AUSTIN,  TEXAS  78721 


IF  PLOT  OPTION  IS  SELECTED,  THE  FOLLOWING  READS  ARE  NOW  NECESSARY. 

READ  - (LABELS)  - r LABVER (1) , 1 = 1,  11] 

LABVER  - Vertical  labels  for  plot  (11  labels) (2  cards) 

READ  - (LABELS)  - r LABHOR(I) , 1=1,  21] 

LABHOR  - Horizontal  labels  (21  labels) (3  cards) 

READ  - (CARD  IMAGE)  - 2 cards  to  identify  plot 

READ  - (LABELS)  - Master,  Name 

Master  - Master  Name  for  Fastrand  file 

Name  - Version  (or  subname)  for  Fastrand  file 

READ  - (FLOATING)  - XMINP,  XMAXP , YMINP , YMAXP 
XMINP  - Minimum  value  of  X for  plot, 

X-direction  is  vertical,  minimum  at  top. 

XMAXP  - Maximum  value  of  X for  plot, 

X-direction  is  vertical,  maximum  at  bottom. 

YMINP  - Minimum  value  of  Y for  plot, 

Y-direction  is  horizontal,  minimum  at  left. 

YMAX P - Maximum  value  of  Y for  plot, 

Y-direction  is  horizontal,  maximum  at  right. 

READ  - (INTEGER)  - INCP 

INCP  - Printing  and  drum  tag  for  number  of  time  steps  to 
be  printed  and  read  out  to  drum  (e.g.,  INCP=5 
means  every  fifth  time  step  will  be  printed) 

READ  - (INTEGER)  - [KEY  (I),  1=1,  7] 

KEY  - Random  number  generator  key  for  probability  of 

detection  for  Barrier  Boat  (7  numbers)  (one  card) 
KEY1  - Random  number  generator  key  for  probability  of 

detection  for  Penetrator  Boat  (7  numbers)  (one  card) 

READ  - (INTEGER)  - [ JEY  (I),  1=1,  7] 

JEY  - Random  number  generator  key  for  random  ZIG-ZAG  course 
for  Penetrator  (7  numbers)  (one  card) 
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READ  - (FLOATING)  - r HEY  (I),  1=1,  7] 

HEY  - Random  number  generator  key  for  normal  distribution 
curves  used  in  SUBSYS.  (7  numbers)  (two  cards) 

READ  - (FLOATING)  - TO,  DELT 

TO  - Initial  value  of  time  variable  (min) 

DELT  - Increment  for  time  step  (min) 

READ  - (FLOATING)  - XMAX,  YMAX , ZMAX 

XMAX  - Value  of  X-coordinate  which  determines  barrier 
boundaries  (yds) 

YMAX  - Value  of  Y-coordinate  which  determines  barrier 
boundaries  (yds) 

ZMAX  - Maximum  depth  allowed  under  any  circumstances 
(negative)  (yds) 

READ  - (INTEGER)  - NSYSTM  (I),  1=1,  5] 

NSYSTM  - Array  of  indexes  for  output  which  defines  which 

system  is  being  used  for  this  run  for  the  following 
subscripts 
NSYSTM  (1) 

NSYSTM  (2) 

NSYSTM  (3) 

NSYSTM  (4) 

NSYSTM  (3) 

For  all  indexes: 

1 = BQS  4 

2 = BQS  4 

3 = BQR  2B 

4 = ROVER 

5 = PUFFS 

6 = BQR  2 (DIMUS) 

7 = BQH  2C  CLASSIFIER 
(5  numbers)  (one  card) 


= Barrier  Detection  Index 
= Barrier  Passive  Ranging  System  Index 
= Barrier  Passive  Bearing  System  Index 
= Barrier  Active  Ranging  System  Index 
= Barrier  Classification  System  Index 

MODIFIED  ACTIVE 
ACTIVE 
PASSIVE 
PASSIVE 
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READ  - 


READ  - 


READ 


READ 


(INTEGER)  - ISYSTT , ISYSRR 

ISYSTT  - Flag  for  tracking  system  bearing  equipment 
ISYSRR  - Flag  for  tracking  system  ranging  equipment 


1 = BQS  4 

2 = BQS  4 

3 = BQR  2B 

4 = ROVER 

5 = PUFFS 


MODIFIED 


ACTIVE 

ACTIVE 

PASSIVE 

PASSIVE 


READ  - 


6 = BQR  2 (DIMUS) 

7 = BQH  2C  CLASSIFIER 

(FLOATING)  - TDETLM  (I),  1=1,  2],  ATTENP  (1) 

DETLM  - Detection  thresholds  (dB) 

ATTENP  - The  absorption  coefficient  used  in  the  Knudsen 
depth  correction  (dB/Yard) 

(INTEGER)  - ISEAS 
ISEAS  - Sea  State  number 

(FLOATING)  - [ F (I , J) , DF(I,J),  DI(I,J),  J=7],  1=1,  2^ 

£ - Center  frequency  for  each  system  of  each  boat  (Hz) 

DF  - Bandwidth  for  each  system  of  each  boat  (Hz) 

PI  - Directivity  Index  for  each  system  of  each  boat  (dB) 

Note:  2 sets  of  7 cards  each,  3 numbers  per  card 

(FLOATING)  - XMIN , YMIN,  ZMIN 

XMIN  - Value  of  X-coordinate  which  determines  barrier 
boundaries  (yds) 

YMIN  - Value  of  Y-coordinate  which  determines  barrier 
boundaries  (yds) 

ZMIN  - Value  of  Z-coordinate  which  determines  barrier 
boundaries  (yds) 

(3  numbers,  1 card) 
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READ  - 


READ  - 


READ  - (FLOATING)  - BLIND (1,1),  BLIND (2, I),  1=1,  NBOAT] 

BLIND  - For  each  boat,  blind  spot  definition  relative  to 
ships  headings  (smallest  value  first) 

Note:  2 cards  each  with  2 numbers 

READ  - (INTEGER)  - [NPASSD(I),  1=1,  2J 

NPASSD  - Number  of  information  points  for  sliding  average 
to  be  used  in  passive  detection  for  each  boat 
(2  numbers)  (one  card) 

READ  - (INTEGER)  - [MOTION(I),  1=1,  2] 

MOTION  - For  each  boat  an  index  defining  type  of  motion: 

1 = Straight  line  path,  constant  velocity 

2 = Straight  line  path,  constant  acceleration 

3 = Circular  path,  constant  speed 
(2  numbers)  (one  card) 

READ  - ( INTEGER)  - IKIND(I),  1=1,  2] 

IKIND  - Kind  of  boat: 

0 = Nuclear 

1 = Diesel -Electric 
(2  numbers)  (1  card) 

READ  - (FLOATING)  - STIME 

STIME  - Stop  Time,  whenever  time  reaches  this,  value 
program  stops  (min) 

READ  - (INTEGER)  - [ITACT(I),  1=1,  2] 

ITACT  - Defines  current  tactic  employed  by  each  boat 

FOR  BOAT  NO.  1: 

1 = No  detection  has  occurred 

2 = Gather  data  on  other  boat 

3 = Contact  Lost-Patrol 

4 = Contact  Lost-Chase 

5 = Fire  Weapon 

6 = Ambush 

7 = Tail-Chase 

8 = Evade 

9 = Close  and  Attack 
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FOR  BOAT  NO.  2: 

1 = No  detection  has  occured 

2 = Evade 

(2  numbers)  (i  card) 

READ  - ( INTEGER)  - ISNORK(I),  1=1,  NBOAT] 

ISNORK  - Defines  snorkeling  status  of  each  boat: 

0 = Not  snorkeling 

1 = Snorkeling 

(2  numbers,  1 card) 

READ  - (FLOATING)  - ESCAPE,  TTURN 

ESCAPE  - The  distance  from  the  boundary  at  which  the 

Penetrator  has  definitely  escaped  the  Barrier 
Boat  (yds) 

TTURN  - Time  for  Penetrator  to  execute  a course  change 
(min  ) 

HEAD  - (FLOATING)  - [EDR(I,J,K),  K=l,  2;  J=l,  2;  1=1,  2] 

EDR  - Set  of  Estimates  of  Detection  Range  in  the  absence 
of  any  information  (yds) 

For  each  boat,  set  consists  of  2 cards,  2 numbers. 

On  each  - The  order  of  data  is  non-snorkeling  estimate 
then  snorkeling  estimate  on  each  card. 

(8  numbers)  (4  cards) 

READ  - (FLOATING)  - [EWR(I,J),  J=l,  2;  1=1,  2] 

EWR  - Set  of  Estimates  of  Weapon  Range  (yds) 

For  each  boat,  set  consists  of  1 card  and  2 numbers; 
order  is: 

IBoat's  estimate  of  JBoat's  Weapon  Range 
(4  numbers)  (2  cards) 

READ  - (FLOATING)  - [ESP(I,J),  J=l,  2;  1=1,  2] 

ESP  - Set  of  Estimates  of  Speed  in  the  absence  of  all 
information  (kts) 
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For  each  boat,  set  consists  of  1 card  with  2 numbers; 
order  is: 

IBoat's  estimate  of  JBoat's  speed 
(4  numbers)  (2  cards) 

READ  - (FLOATING)  - [EDR1(I,J),  J=l,  2;  1=1,  2] 

EDRl  - Set  of  Estimates  of  Detection  Range  for  use  in 
time  to  fire  determination  (yds) 

For  each  boat,  set  consists  of  1 card  with  2 numbers; 
order  is : 

IBoat's  estimate  or  JBoat's  detection  range 
(4  numbers)  (2  cards) 

READ  - (FLOATING)  - [EWR1(I,J),  J=l,  2;  1=1,  2] 

EWRl  - Set  of  Estimates  of  Weapon  Range  for  use  in  time 
to  fire  determination  (yds) 

For  each  boat,  set  consists  of  1 card,  2 numbers; 
order  is : 

IBoat's  estimate  of  JBoat's  Weapon  range 
(4  numbers)  (2  cards) 

(FLOATING)  - LZT(J,I),  J=l,  9;  1=1,  2] 

ZT  - Normal  depth  of  each  boat  during  each  tactical 
situation  (negative)  (yds) 

For  each  boat,  set  of  9 numbers,  (3  cards),  (4  numbers 
per  first  2 cards,  1 number  3rd  card)  / 

(18  numbers)  (6  cards) 

(FLOATING)  - [ZS(I),  1=1,  2] 

ZS  - The  depth  at  which  each  boat  snorkels  (negative) 
(yds),  (2  numbers)  (1  card) 

(FLOATING)  - [ ZLIM(J , I) , J=l,  2;  1=1,  2] 

ZL1M  - For  each  boat,  a two  component  vector  which 

defines  the  maximum  and  minimum  operating  depth, 
(max  value  first)  (negative)  (yds) 

(4  numbers)  (2  cards) 
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READ  - (FLOATING)  - . XMAN(I),  1=1,  2J 

XMAN  - Minimum  and  maximum  values  of  X inside  of  which 
Boat  1 must  maneuver  (small  value  first)  (yds) 

(?  numbers) 

READ  - (FLOATING)  - YMAN(I) , 1=1,  2] 

YMAN  - Minimum  and  maximum  value  of  Y inside  of  which  Boat 

1 must  maneuver  (small  value  first)  (yds)  (2  numbers) 

READ  - (FLOATING)  - WRMAX(I),  1=1,  2] 

WRMAX  - Maximum  range  of  each  boat's  weapons  (yds) 

(2  numbers) 


READ  - (FLOATING)  - LXPAT(I),  1=1,  2J 

XPAT  - X coordinate  of  minimum  and  maximum  patrol  points 
for  Boat  No.  1 (small  value  first)  (yds) 

(2  numbers) 


READ  - (FLOATING)  - 'CT(J.I),  J-l , 9;  1=1,  2] 

Cl  - Set  of  normal  speeds  during  tactical  situations  (kts) 

For  each  boat,  set  consists  of  9 numbers,  3 cards,  (4 

numbers  each  for  first  2 cards,  3rd  card,  1 number) 
(18  numbers)  (6  cards) 

READ  - (FLOATING)-  CS(J,I),  J=l,  9;  1=1,  2] 

CS  - Set  of  speeds  for  each  boat  during  each  tactical 
situation  while  snorkeling  (kts) 

For  each  boat,  set  consists  of  9 numbers,  3 cards  (4 
numbers  each,  first  2 cards,  3rd  card,  1 number) 

(18  numbers)  (6  cards) 

READ  - (FLOATING)  - [CMIN(I),  1=1,  2] 

CM IN  - For  each  boat  minimum  speed  required  to  maintain 
steerage  (kts) 

( 2 numbe  rs ) 

READ  - (FLOATING)  - iCMAX(I),  1=1,  2] 

CMAX  - For  each  boat  maximum  speed  (kts) 

(2  numbers) 

UNCLASSIFIED 


READ  - (FLOATING)  - [PHI02(I),  1=1,  2] 

PHIQ2  - Basic  course  of  IBOAT  (degrees) 

(2  numbers)  f 

READ  - (FLOATING)  - [ RFIRE(I) , 1=1,  2] 

RFIRE  - Maximum  firing  range  of  each  boat  (yds) 

(2  numbers) 

READ  - (FLOATING)  - [TDLY(J,I),  J=l,  9;  1=1,  2] 

TDLY  - For  each  boat  a vector  which  specifies  a set  of 
delay  times  required  to  react  to  each  tactical 
situation  change  (min) 

For  each  boat,  set  consists  of  (9  numbers,  3 cards) 

(4  numbers  each  first  2 cards,  3rd  card,  1 number) 

(18  numbers)  (6  cards) 

READ  - (FLOATING)  - [TSS(I),  1=1,  2] 

TSS  - For  each  boat,  time  to  start  next  snorkel  period 
(min) 

(2  numbers) 

READ  - (FLOATING)  - [TCS(I),  1=1,  2] 

TCS  - For  each  boat,  time  to  stop  next  snorkel  period 
(min) 

(2  numbers) 

READ  - (FLOATING)  - [ TLSNK(I) , 1=1,  2] 

TLSNK  - For  each  boat,  time  that  boat  completed  last 
snorkel  (min) 

(2  numbers) 

READ  - (FLOATING)  - [SNBC(J,I),  J=l,  9;  1=1,  2] 

SNBC  - For  each  boat,  a fraction  representing  the  battery 
power  level  requiring  snorkeling  during  tactical 
situations 


For  each  boat,  set  consists  of  9 numbers,  3 cards  (4 

numbers  each,  first  2 cards;  third  card,  1 number) 
(18  numbers)  (6  cards) 
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READ  - (FLOATING)  - [ETS(I),  1=1,  2] 

ETS  - Estimate  of  snorkel  frequency  (min) 

(2  numbers) 

READ  - (FLOATING)  - [BSO(J,I),  J-l,  7;  R0(1,I),  R0(2,I),  R0(3,I), 

1=1,  2] 

BSO  - For  each  boat,  a seven  component  vector  which  con- 
tains the  motion  parameters  at  time  TO 
BS0(1,NB0AT)  = Horizontal  linear  speed  (kts) 
BS0(2,NB0AT)  = Heading  (degrees) 

BS0(3,NB0AT)  = Horizontal  circular  speed  (kts) 
BS0(4,NB0AT)  = Radius  of  circle  of  motion  (yds) 
BS0(5,NB0AT)  = X-coordinate  of  center  of  circle  (yds) 
BS0(6,NB0AT)  = Y-coordinate  of  center  of  circle  (yds) 
BS0(7,NB0AT)  = Horizontal  linear  acceleration  (kts/min) 
RO  - For  each  boat,  a three  component  vector  which  defines 
the  initial  positions  (X,Y,Z)  respectively  (yds) 

For  each  boat,  set  consists  of  7 motion  parameters,  then 
3 position  parameters 

(1st  card  - 4 numbers,  2nd  card  - 3 numbers, 

3rd  card  - 3 numbers) 

(20  numbers)  (6  cards) 
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OUTPUT 

1.  SHIP  POSITIONS  (For  both  Barrier  Boat  and  Penetrator 

Boat: 

TIME  = Time  (min) 

^ HEADING  = Angle  from  North  to  boat's  axis  (degrees) 

SPEED  = Speed  of  boat  (kts) 

^ (X,Y,Z)  = Coordinates  of  boat  in  barrier 

coordinate  system  (yds) 

RHO  = True  range  between  boats  (yds) 

?.  SUBSYSTEM  VALUES  (BARRIER  BOAT  VALUES) 

TIME  = Time  (min) 

TACTIC  = Tactical  situation  for  Barrier  Boat 

SNORK  = Flat  to  indicate  if  Barrier  Boat  is 

snorkeling  (=1)  or  not  (=0) 

AVE  S/N  = The  maximum  average  signal-to-noise 

ratio  for  the  detection  subsystem  at 
the  time  (dB) 

S/N  = Signal-to-noise  level  computed  at  the 

time  by  the  S to  N routine  as  seen  at 
the  detection  subsystem  (dB) 

PROB  D = Probability  of  detection 

D.  BEARING*  = Relative  bearing  reported  by  the 
AN/BQR-2(DIMUS)passive  bearing 
subsystem  (degrees) 

D.  B ERR  = Error  in  the  bearing  reported  by  the 
DIMUS  subsystem  (degrees) 

P.  BEAR*  = Relative  bearing  reported  by  PUFFS 
(AN/BQG-4)  subsystem  (degrees) 

P.B  ERR  = Error  in  the  bearing  reported  by 
PUFFS  (degrees) 


Negative  value  means  signal  too  low  for  subsystem  to  provide 
information. 
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P.  RANGE* 


P.  R ERR 


= Range  reported  by  the  PUTTS  subsystem 
(yds) 

= Error  in  range  reported  by  PUFFS  (yds) 


SUBSYSTEM  VALUES  (PENETRATOR  BOAT  VALUES) 

TIME  = Time  (min) 

TACTIC  = Tactical  situation  for  Penetrator 

Boat 

SNORK  = Flag  to  indicate  if  Penetrator  is  snorkeling 

(=1)  or  not  (=0) 

AVE  S/N  = Maximum  average  signal- to-noise  ratio 
for  the  detection  subsystem  at  the 
time  (dB) 

S/N  = Signal- to-noise  ratio  currently  available 

to  the  detection  subsystem  of  the 
Penetrator  (dB) 

PROB  D = Probability  of  detection 

BEARING*  = Relative  bearing  reported  by  Penetrator' s 

bearing  subsystem  (degrees) 

BEAR  ERR  = Error  in  error  reported  by  bearing 
subsystem  (degrees) 

RHO  = Actual  range  between  boats  (yds) 

SUBSYSTEM  VALUES  (Barrier  Boat,  Tracking  Information) 

TIME  = Time  (min) 

X(EST.)  = Current  estimated  X-coordinate  of 

target  on  base  course  (yds) 

Y(EST.)  = Current  estimated  Y-coordinate  of 

target  on  base  course  (yds) 

SPEED  (EST.)=  Current  estimate  of  target's  speed 
along  base  course  (kts) 


Negative  value  indicates  signal  level  too  low  for  subsystems 
to  provide  information. 
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= Number  of  points  stored  in  base 
course  track 

= Flags  indicating  that  X or  Y motion 
estimates  satisfy  (=1)  criteria  or 
don't  (=0) 

= Estimate  of  X-velocity  of  target  (kts) 

= Estimate  of  Y-velocity  of  target  (kts) 

= Standard  deviation  of  X-data 
about  least  squares  line  (yds) 

= Standard  deviation  of  Y-data  about 
least  squares  line  (yds) 

SUBSYSTEM  VALUES  (Barrier  Boat,  Fire  Control  Information) 
(Same  as  TRACKING  INFORMATION  DEFINITIONS  EXCEPT 
current  course  must  be  substituted  for  Base  course) 
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YVEL(EST) 
SIGMA  X 
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A. 2 POSIT 


PURPOSE:  This  subroutine  computes  new  position  and  velocity 

coordinates  for  each  boat  as  well  as  computing  the  range,  true 
bearing,  heading,  relative  bearing,  and  other  data  used  in  the 
program.  The  subroutine  first  integrates  the  equations  of  motion 
for  each  boat,  then  computes  the  required  data. 

INPUT: 

TIME  = Time 

DELT  = Length  of  time  step  (min) 

NBOAT  = Number  of  boats  in  the  simulation 

MOTION  = Array  describing  the  type  of  motion  which  each 

boat  is  effecting  1 = Straight  line  path,  constant 

velocity 

2 = Straight  line  path,  constant 

acceleration 

3 = Circular  path,  constant 

speed 

BS(I,IBOAT)  = Array  of  motion  parameters  for  each  boat  IBOAT 


1=1  Horizontal  linear  speed  (kts) 

1=2  Heading  (degrees) 

1=3  Horizontal  circular  speed  (kts) 

1=4  Radius  of  circular  motion  (yds) 

1=5  X-coordinate  of  center  of  circle  (yds) 
1=6  Y-coordinate  of  center  of  circle  (yds) 
1=7  Horizontal  linear  acceleration  (kts/min) 


OUTPUT: 

RDOT(J , IBOAT)  = J-component  of  IBOAT' s velocity  vector 

(kts) 

R(J, IBOAT)  = A three  component  vector  which  defines  the 

current  position  (X,Y,Z)  of  boat  IBOAT  (yds) 
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RDDOT(J , IBOAT)  = Three  component  vector  defining  the  linear 

acceleration  of  boat  IBOAT  (kts/min) 

RHO  = Range  between  IBOAT  and  JBOAT  (yds) 

PHIRHO  = True  bearing  of  JBOAT  from  IBOAT  (radians) 

RELB  = Relative  bearing  of  JBOAT  from  IBOAT 

(radians) 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

PSI  = Angle  which  radius  vector  to  IBOAT  makes 

with  North,  in  a clockwise  sense,  at  the 
beginning  of  the  time  step  (radians) 

ARGC  = Angle  which  radius  vector  to  IBOAT  makes 

with  North,  in  a clockwise  sense,  at  the 
end  of  the  time  step  (radians) 

Conversion  from  knots  to  yards  per  minute  = 33.783 
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A.  3 SUBSYS 


PURPOSE:  This  subroutine  evaluates  the  various  subsystems  for 

each  of  the  boats  in  the  simulation  and  stores  the  outputs  in 
labeled  common  for  reference  by  the  TACTIC  subroutine  (i.e.,  for 
use  by  the  boats  to  determine  their  future  course  of  action) . 

INPUT: 

TIME  = Time  (min) 

NBOAT  = Number  of  boats  in  the  simulation 
ISYSON  = Array  of  flags  indicating  whether  a subsystem  on 
a particular  boat  is  on  (=1)  or  off  (=0) 

ITYPE  = Array  indicating  the  role  of  IBOAT  in  the  simulation 

0 = Barrier  Boat 

1 = Penetrator 

NSS  = Array  containing  total  number  of  subsystems 

available  to  each  boat  according  to  the  definitions 
of  the  first  argument  listed  in  control 
IQUIP  = Tag  describing  equipment  available  on  a boat 

0 = does  not  possess  equipment 

1 = does  possess  equipment 

IDETEC  = Flag  denoting  if  boat  has  detected  anyone 

0 = no 

1 = yes 

DLSN  = Threshold  of  useful  information  for  the  detection 
subsystem  for  IBOAT  (dB) 

RELB  = Array  of  relative  bearing  of  JBOAT  from  IBOAT, 
where  IBOAT  is  the  first  index  (radians) 

BLIND  = Array  of  blind  spot  angles  for  IBOAT  (radians) 

ISTART  = Flag  to  instruct  each  boat  whether  or  not  to  clear 
passive  detection  information  field  and  rebuild 
(0  = rebuild).  If  non-zero,  it  indicates  the 
index  of  the  oldest  S/N  in  the  probability  of 
detection  table 
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NPASSD 

DETLM 

LST 

DELT 

RHO 

DI 

HEY 

PHIRHO 

R 

ISYSRR 

ISY3TT 

OUTPUT : 

STN 

PDF 

PD 

SNAV 

NTIME 


Array  listing  the  number  of  points  required  to 
fill  the  S/N  table  used  for  passive  detection 
Lower  limit  for  probability  of  detection  below 
which  there  is  assumed  to  be  no  chance  that 
detection  can  occur 

Length  of  time  for  which  S/N  is  less  than  thres- 
hold before  detection  is  assumed  to  be  lost  (min) 
Length  of  each  time  step  (min) 

Actual  distance  between  IBOAT  and  JBOAT  (yds) 
Receiving  directivity  index  for  system  NSYS 
(negative)  (dB) 

Array  of  random  numbers  used  as  keys  for  random 
number  generator 

True  bearing  of  JBOAT  with  respect  to  IBOAT 
(radians) 

Vector  of  the  true  X,  Y,  and  Z components  of 
IBOAT  in  the  Earth -fixed  barrier  coordinate  system 
Number  of  subsystem  providing  range  information 
to  FIRCON  and  TRACK 

Number  of  subsystem  providing  bearing  information 
to  FIRCON  and  TRACK 

Array  of  signal- to-noise  ratio's  indexed  by  IBOAT, 
JBOAT  and  NSYS  (dB) 

Us^d  to  store  array  of  S/N  for  each  averaging 
process  for  each  time  step 

Array  used  to  store  maximum  probability  of 
detection  for  all  averaging  processes  at  each 
time  step.  The  maximum  value  is  stored  in  this 
array  upon  exit  from  the  routine 
Array  of  the  maximum  average  S/N  computed  at  the 
current  time  step 

Number  of  points  currently  in  the  table  of  S/N 

ratios  used  to  determine  detection,  up  to  the 

maximum  allowable  number 
A 62 
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= Number  of  consecutive  time  steps  for  which 
detection  has  been  lost 

= Distance  between  IBOAT  and  JBOAT  reported  by 
ranging  system  NSYS  (yds) 

= Actual  error  in  DR  (yds) 

= Relative  bearing  of  JBOAT  from  IBOAT  reported 
by  system  NSYS  (radians) 

= Actual  error  in  THET  (radians) 

= Time  since  the  classification  subsystem  was 
updated  (min) 

= Length  of  time  for  which  classification 
subsystem  has  been  continuously  gathering 
information  (min) 

= Flag  indicating  that  IBOAT  has  classified  JBOAT 
with  0.90  probability 

= Flag  indicating  that  JBOAT  has  detected  IBOAT' s 
ping 

= Tracking  flags  indicating  the  status  of  the 
track  for  the  fire  control  system 
= Number  of  points  in  current  fire  control  track 
= Number  of  points  in  current  track  to  determine 
target's  base  course 

= Array  of  a -parameters  determined  in  FIRCON 
(yds  and  kts) 

= Array  of  standard  deviations  determined  in 
FIRCON  (yds) 

= Estimated  current  X-position  of  target  in  Earth- 
fixed  barrier  coordinate  system  from  FIRCON  (yds) 
= Estimated  current  Y-position  of  target  in  Earth- 
fixed  barrier  coordinate  system  from  FIRCON  (yds) 
= Time  track  was  last  updated  in  FIRCON  (min) 

= Array  of  a-parameters  determined  in  TRACK  (yds 
and  kts) 

= Array  of  standard  deviations  determined  in  TRACK 
(yds) 
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MFLGTK  = Flags  indicating  whether  or  not  TRACK  has 

accurately  determined  the  base  course  of  the 
target 

XESTTK  = Estimated  current  X-position  of  target  on  base 
course  in  Earth- fixed  barrier  coordinate  system 
(yds) 

YESTTK  = Estimated  current  Y-position  of  target  on  base 
course  in  Earth- fixed  barrier  coordinate  system 
(yds) 

TLSTTK  = Time  track  was  last  updated  in  TRACK  (min) 


INTERNAL  VARIABLES  AND  CONSTANTS: 

IBOAT  = Tag  indicating  the  boat  on  which  the  subsystem 
is  located 

JBOAT  = Tag  indicating  the  boat  about  which  information 
is  being  gathered 

NSYS  = Tag  indicating  the  number  of  the  subsystem 

under  consideration 

NSYSMX  = Upper  limit  of  do  loop  for  each  type  of  subsystem 

IDDETl  = Number  of  time  steps  for  which  no  useful 

information  was  accumulated  prior  to  detection 
by  Barrier  Boat 

IDDET2  = Serves  same  function  for  Penetrator  that  IDDETl 
serves  for  Barrier  Boat 

INI  = Least  number  of  time  steps  before  a probability 

of  detection  is  computed 

II  = Number  of  S/N  that  will  be  averaged 

SUM  = Average  S/N  for  the  interval  characterized  by  II 

PDD  = Probability  of  detection  returned  by  PROBD 

NRAN  = Random  number  between  1 and  2000. 

PDT  = 1000  times  the  maximum  probability  of  detection 

for  the  current  time  step 

Increment  by  which  DELTM  is  increased  after  each  test  for 
detection  is  0.05 
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SNT  = Signal- to-noise  ratio  at  JBOAT  when  IBOAT  pings 

with  an  active  sonar 

IXFIR  = Flat  determining  whether  or  not  a fire  control 
quality  track  has  been  achieved  1 = no 

2 = yes 

X = True  bearing  of  JBOAT  with  respect  to  IBOAT 

reported  by  bearing  subsystem,  including  error 
(radians ) 
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ATERPO 

PURPOSE:  The  function  ATERPO  is  a routine  for  straight  line 

interpolation  using  the  slope-Y-intercept  form  for  a straight 
1 ine . 

INPUTS:  There  is  one  formal  argument  to  the  routine,  ATERPO(B). 

B is  the  Y-intercept  for  the  interpolation.  All  other  inputs 
are  in  the  form  of  a labeled  common  block  called  SNTAB2 , the 
variables  of  which  are  X,  S,  T,  D. 

X is  the  X-value  of  the  point  to  be  interpolated 

S and  T are  the  two  functional  values  used  in  the 
interpolation 

D is  the  difference  in  abscissas  of  the  two  functional 
values 

OUTPUT:  The  output  is  the  interpolated  value 
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BQRANG 


PURPOSE:  This  subroutine  returns  the  relative  bearing 

information  determined  by  the  AN/BQR-2B  and  the  AN/BQR-2  (DIMUS) . 
This  is  accomplished  by  taking  the  actual  relative  bearing  between 
the  two  boats  in  question  and  adding  to  it  a random  error 
characteristic  of  the  bearing  errors  of  the  2B  and  2 (DIMUS) 
equipment . 

INPUT: 


R 

TIJ 

SN 

HEY 

NSYS 


OUTPUT: 


= Actual  range  between  the  two  boats  (yds) 

= Actual  relative  bearing  between  the  two 
boats  (radians) 

= Signal- to-noise  ratio  at  the  output  of  the 
beamformer  (dB) 

= A dimensioned  array  of  random  numbers 
= A flag  for  determining  whether  the  2B  or 

the  2 (DIMUS)  is  to  be  used  to  return  information 
3 = 2B 

6=2  (DIMUS) 


BIJ  = Relative  bearing  returned  by  2B  or  2 (DIMUS) 

including  error  (radians) 

DB  = Bearing  error  characteristic  of  the  2B  or 

2 (DIMUS)  (radians) 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

DR  = Conversion  factor  from  degrees  to  radians 

RD  = Recognition  threshold  below  which  equipment 

is  inoperable  (dB) 

Blind  area  is  a sector  from  165°  to  195° 
bearing  in  which  the  equipment  can  return  no  useful  information. 
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ENTER 


NO  INFORMATION 
AVAILABLE 

B I J ■ -~1  ,~ 

DB-0.0 


RETURN 


TARGET 

IN 

3LIN0  SPOT 


XBAR-DR 
SIGMA- . 3*DR 


XBAR-0 

SIGMA*. 25*DR 


_CALLl  RANDN J 

r ANOOMLY  SELECT 
DB  FROM 
fi(XBAR, SIGMA) 


RETURN 
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BS4RAN 


PURPOSE:  The  purpose  of  this  subroutine  is  to  return  the 

range  determined  by  the  BQS-4  and  the  BQS-4  (MOD  2).  This  is 
accomplished  by  adding  to  the  actual  range  a random  range  error 
which  is  characteristic  of  the  particular  system  being  used. 


INPUT: 


R 

= 

Actual  range  between  the  two  boats  (yds) 

TIJ 

Actual  relative  bearing  as  seen  by  own  ship 
(radians ) 

SN 

= 

Signal- to-noise  ratio  at  own  ship,  measured 
at  output  of  the  beamformer  (dB) 

HEY 

= 

A dimensioned  array  of  random  numbers 

TJI 

— 

Actual  relative  bearing  as  seen  by  target  ship 
(radians ) 

I BOAT 

= 

Tag  identifying  own  ship  (ship  gathering 
information) 

JBOAT 

= 

Tag  identifying  target  ship  (ship  which  is 
object  of  information  being  gathered) 

BLIND 

= 

A dimensioned  array  of  the  blind  spot  angular 
limits 

IMOD 

OUTPUT: 

Tag  indicating  equipment  to  be  modeled 

1 = BQS-4 

2 = BQS-4  (MOD  2) 

I FLAG 

= 

Flag  indicating  whether  JBOAT  has  been  alerted 
by  the  ping  of  IBOAT  (1  = alert,  0 = not  alert) 

SNT 

= 

Signal- to-noise  ratio  at  JBOAT  as  a result  of 
IBOAT' s ping  (dB) 

RIJ 

= 

Range  including  error  returned  by  the  active 
system  (yds) 

DR 

= 

Range  error  characteristic  of  the  active  system 
used  (yds) 
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INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 


HALPD 


RD 


IM0D1 


Signal- to-noise  level  for  which  the  acti.ve  sonars 
have  a 0.50  probability  of  detection  (dB) 
Recognition  threshold  for  JBOAT  to  detect  IBOAT's 
ping  (dB) 

Tag  used  by  subprogram  STON  to  calculate  the 
term  SNT 
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CLASA 


PURPOSE:  The  purpose  of  this  subroutine  is  to  model  the 

BQH-2C  classification  subsystem. 


INPUT: 

SN 


DELT 


= The  signal- to-noise  ratio  at  the  classification 
subsystem  (dB) 

= The  time  increment  at  which  the  classification 
subsystem  is  updated  (min) 


OUTPUT : 

TLOSS 

TCLASS 

ICLASS 


= The  time  since  the  classification  subsystem  was 
updated  (min) 

= The  time  over  which  useful  classification  in- 
formation has  been  accumulated  (min) 

= Flag  which  is  set  when  the  probability  of 
classification  is  0.90  or  higher 


INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 
THRESHOLD  FOR  CLASSIFIER  = -6  dB 

UPPER  LIMIT  FOR  TLOSS  = 1 min 

LOWER  LIMIT  FOR  CLASSIFICATION  = 20  min 
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INPUT 


INCREMENT 

, NO 

IS 

SN> 

\ YES 

INCREMENT 

TLOSS 

I 

THRESHOLD  OF 

.classifier^ 

^ ^ 

T CL ASS 

' TLOSS> 
UPPER  LIMIT 
FOR  TLOSS^ 


Rt I URN 


CLASSIFICATION 
MUST  BE 
RESTARTED 

TLOSS^O 

TCLASS:0 


TCLASS> 

LOWER  LIMIT 
FOR  CLASSIFICATION, 


CLASSIFICATION 

ACHIEVED 

I CLASS = I 
TLOSS-O 
TCLASS=0 


RETURN 


FIG  A-6-CLASA  FLOW  CHART 
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FNOISE 


PURPOSE:  Function  subprogram  FNOISE  determines  the  total 

masking,  background  noise  seen  by  the  receiving  sonar. 


INPUTS: 

I SEAS 
IKINDI 

FI ISON 
ATTENP 

RI 

BSIBT 

DIISON 

DFIISN 

JSNORK 


Sea  State  number 

Type  of  boat  0 = Nuclear 

1 = Diesel-Electric 

Operating  frequency  (center  band)  (Hz) 
Absorption  coefficient  for  propagation  of  sound 
through  sea  water  (dB/yd) 

Depth  of  IBOAT  (negative  yds) 

Velocity  of  IBOAT  (kys) 

Receiving  directivity  index  of  sonar  (negative) 
(dB) 

Receiving  bandwidth  of  sonar  (Hz) 

Flag  for  snorkel  0 = No 

1 = Yes 


OUTPUT: 

The  value  of  FNOISE  is  the  only  output  of  the  routine 
INTERNAL: 

SEANOI  and  SEFNOI  are  called  internally. 
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FIG.  A-7  - FNOISE  FLOW  CHART 
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NOISY 


PURPOSE:  The  subroutine  NOISY  is  a three-dimensional  inter- 

pola  ion  /extrapolation  routine  to  find  a Y-value  for  a given  X 
and  Z.  The  Y-value  is  determined  by  interpolation  in  or  extrapola- 
tion of  a set  of  tabular  values. 


INPUTS:  The  inputs  to  NOISY  are:  3 formal  arguments,  and  a 

common  block  SNTAB1  with  7 variables.  The  formal  arguments  are: 

1.  Starting  addr*ess  of  two-dimensional  array  of 
tabular  values.  The  tabular  values  are  Y-values. 

2.  Starting  address  of  table  of  X-values  which 
correspond  to  the  first  dimension  of  the  Y-tabular 
values . 

3.  Size  of  first  dimension  of  Y-tabular  values,  i.e., 
number  of  X-values. 

The  common  block  variables  are: 

1.  X-value  for  interpolating/extrapolating 

2.  Z-value  for  interprolating/extrapolating 

3.  Index  for  start  of  values  in  X-array  and  first 
dimension  of  Y-arrays 

4.  Index  for  last  of  values  in  X-array  and  first 
dimension  of  Y-arrays 

5.  Array  (dimensioned  4)  of  Z-values  corresponding  to 
2nd  dimension  of  Y-tabular  values  array.  Z-values 
are  step  widths  (or  increments)  between  tables. 

6.  Number  of  Z-values  for  Z-array  and  2nd  dimension 
of  Y-tabular  values  array 

7.  Variable  for  output. 
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OUTPUTS:  The  only  output  from  NOISY  is  that  of  the  interpolated/ 

extrapolated  Y-value  which  is  returned  as  the  7th  variable  in 
the  common  block  SNTABl. 

INTERNAL:  NOISY  calls  the  function  ATERPO  to  do  all  the  inter- 
polation/extrapolation. The  inputs  for  ATERPO  are:  1 formal 

argument,  and  a common  block  SNTAB2  with  4 variables. 

The  formal  argument  is  the  Y- intercept  value  for  the 
straight  line  interpolation. 

The  common  block  variables  are: 

1.  X-value  to  be  interpolated  for, 

2&3.  Two  tabular  values  to  interpolate  or  extrapolate 

from, 

4.  Difference  between  X-values  of  two  tabular  values. 
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PURPOSE:  This  function  subprogram  determines  the  propagation 

loss  from  a source  to  a receiver  if  the  sonar  mode  is  passive, 


and  from  the  source  to  target  to  receiver  if  the  sonar  mode  is 
active.  This  is  accomplished  by  a table  look-up  technique,  and  the 
values  stored  in  the  tables  are  negative.  The  values  stored  in 
the  tables  are  those  for  the  path  between  the  two  boats  which 
permits  the  highest  intensity  signal  to  reach  the  receiver. 

INPUT: 

ISONDX  = Index  indicating  the  type  of  sonar  equipment 
RI  = Depth  of  IBOAT  (negative)  (yds) 

RJ  = Depth  of  JBOAT  (negative)  (yds) 

RHOIJ  = Range  of  JBOAT  from  IBOAT  (yds) 

INVERS  = Flag  for  inverse  interpolation  0 = No 

1 = Yes 

OUTPUT: 

The  value  of  PRLOSS  is  the  only  output  of  the  routine 
INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

CDEPTH  = Layer  Depth  = -33.3  (yds) 

The  variables  of  common  block  SNTABl  must  be  set  before 
calling  NOISY. 

FREQUY  = Range  in  kyds 

VELOCY  = Not  used  =0.0 

NFIRST  = Starting  index  in  PROPAGATION  LOSS  TABLE 
NLAST  * Last  index  in  PROPAGATION  LOSS  TABLE 
STEPWD(l)=  0.0  (not  used) 

SOUND  = Output  from  NOISY 

For  INVERSE  INTERPOLATION,  the  input  variable  RHOIJ  must  be 
the  propagation  loss  (multiplied  by  1000.0),  the  call  to  NOISY 
has  the  X and  Y arrays  reversed. 

NOISY  has  three  formal  arguments: 

1.  Starting  address  of  X-array 

2.  Starting  address  of  Y-array 

3.  Size  of  first  dimension  of  X-array 
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ENTER 

PRLOSS) 


both  above 

LAYER 


DIB  = -33  3-  R1 
DJB--333-RJ 


BOTH  BELOW 
LAYER 


IS  = 1SONDX 


: ONE  ABOVE  LAYER  -° 
1 ONE  BELOW 


ISONDX-6 


[SONDX-7 


ISONPX 


PRINT 

ERROR 

MESSAGE 


N - NSTART(IS) 

I PRE0UY=  RHOIJ/IOOOO 
I VEL0CY  = 0.0 

NFIRST=I 

NLA  iT  = NS  TART  ( IS  -tl)-N 
STEPWD(I)  -O.C 


F RE  QU Y - RHOI J / IOOO  0 
VELOCY-  0 0 
ntable=i 

NFIRST  = I 


ISONDX-IC 


ISONDX-I 


NTABLE=I 
NSIZE  = 
NLAST-H 


NLAST=3I 

NSIZEr3l 

N=3l 


NLAST=7 

NSIZE=7 

N=l 


INVERS. 


/ CALL  NOISY 
(STRAIGHT 
\ INTERPOLATION)/ 


INVERS 


W NOISY  \ 
((INTEPOLATION) 
\ FOR  / 
Xtorpedo/ 


PRLOSS- 

SOUND 


/ CALL  N OISY 

(INVERSE 
, INTERPOLATION) 


CALL  NOISY  \ 
(INVERSE 

\ interpolation 
\ FOR  TORPEDO)/ 


RETURN ' 
PRLOSS / 


FIG.  A-9  - PRLOSS  FLOW  CHART 
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PROBD 


PURPOSE:  The  subroutine  PROBD  is  a routine  to  determine  the 

probability  of  detection  given  a signal- to-noise  ratio.  If  the 
maximum  amount  of  usable  detection  data  is  not  available  then  the 
signal- to-noise  ratio  is  reduced.  This  reduction  accounts  for  the 
fact  that  the  tabular  values  are  based  on  the  maximum  amount  of 
detection  data. 

INPUT:  There  are  7 formal  arguments  to  PROBD: 


1. 

ISONDX  - = 

3 - Use  tables  for  BQR-2B 
6 - Use  tables  for  BQR-2  (DIMUS) 

else  - error  halt 

2. 

DLSN 

not  used 

3. 

DELT 

not  used 

4. 

NTIME 

The  time  since  the  beginning  of  the 
simulation  (min) 

5. 

SNAV 

Average  S/N  computed  at  the  current 
time  step  (dB) 

6 . 

NTMAX 

Maximum  amount  of  usable  detection  data 

7. 

PDD 

Output 

OUTPUT: 

The  output 

is  PDD  - the  probability  of  detection 

INTERNAL: 

F(x)  - A function  for  the  reduction  of  S/N  when 

the  maximum  amount  of  usable  detection  data  is  not  available. 
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PUFANG 

PURPOSE:  This  subroutine  returns  the  relative  bearing  deter- 

mined by  the  AN/BQG-4  (PUFFS).  This  is  accomplished  by  taking  the 
actual  relative  bearing  between  the  two  boats  in  question  and 
adding  to  it  a random  error  which  is  characteristic  of  the  bearing 
error  of  the  PUFFS  equipment. 

INPUT: 

R = Range  between  the  two  boats  (yds) 

TIJ  - Actual  relative  bearing  between  the  two  boats 

(radians ) 

SN  = Signal- to-noise  ratio  at  the  output  of  the  beam- 

former 

HEY  = A dimensioned  array  of  random  numbers 

DI  = Directivity  Index  of  the  PUFFS  sonar  (currently 

unused) 

OUTPUT: 

BIJ  = Relative  bearing  returned  by  PUFFS  including 

error  (radians) 

DB  = Bearing  error  characteristic  of  the  PUFFS 

equipment  (radians) 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

RD  = Recognition  threshold  below  which  PUFFS  is 

inoperable  (dB) 

Maximum  range  beyond  which  PUFFS  furnishes 
no  useful  information  = 18000  yds.  Angle  sectors  outside  of  which 
information  from  PUFFS  is  disregarded  are  + 40°  of  abeam. 

T = Signal  integration  time  for  PUFFS  correlator 
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PUFRAN 


PURPOSE:  This  subroutine  returns  the  range  between  the  two 

boats  determined  by  the  AN/BQG-4  (PUFFS).  This  is  accomplished 
by  taking  the  actual  range  between  the  boats  and  adding  to  it  a 
random  error  which  is  characteristic  of  the  range  error  of  the 
PUFFS  equipment. 


Actual  range  between  the  two  boats  (yds) 

Actual  relative  bearing  between  the  two  boats 
(radians) 

Signal- to-noise  ratio  at  the  output  of  the 
beamformer  (dB) 

A dimensioned  array  of  random  numbers 

Range  returned  by  PUFFS  including  error  (yds) 
Range  error  characteristic  of  the  PUFFS 
equipment  (yds) 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

RD  = Recognition  threshold  below  which  PUFFS  is 

inoperable  (dB) 

Maximum  range  beyond  which  PUFFS  furnishes  no 
useful  information  = 18000  yds. 

Angle  sectors  outside  of  which  information  from 
PUFFS  is  disregarded  are  + 40°  of  abeam. 


INPUT: 

R 

TIJ 

SN 

HEY 

OUTPUT: 

RIJ 

DR 
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AS  DR 


RIJ  - R + DR 


RETURN 
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RADSIG 


PURPOSE:  The  purpose  of  this  function  subprogram  is  to 

compute  the  passive  signal  received  by  IBOAT  from  JBOAT.  This  is 
accomplished  by  a table  look-up  of  the  spectrum  level  noise 
radiated  at  the  center  frequency  of  IBOAT' s sonar,  adjusting 
this  value  for  the  propagation  loss  between  JBOAT  and  IBOAT  and 
adding  a term  to  account  for  the  bandwidth  of  the  receiver  on 
IBOAT. 

INPUT: 

ISONDX  = Index  indicating  the  type  of  sonar  equipment 
IKINDJ  = Type  of  boat 

0 = Nuclear 


1 = Diesel-Electric 
JSNORK  = Flag  for  snorkel 

0 = No 

1 = Yes 

FIISON  = Geometric  or  center  frequency  of  sonar  ISONDX  (Hz) 

RJ  = Depth  of  JBOAT  (negative) (yds ) 

RI  = Depth  of  IBOAT  (negative) (yds ) 

RHOIJ  = Range  of  JBOAT  from  IBOAT  (yds) 

DFIISN  = Receiving  bandwidth  of  sonar  JSONDX  (Hz) 

BSJBT  = Velocity  of  JBOAT  (kts) 

OUTPUT:  The  value  of  RADSIG  is  the  only  output  of  the  program 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

The  variables  in  common  block  SNTABI  must  be  set  prior  to 
calling  NOISY. 

FREQUY  = Center  frequency  of  sonar 
VELOCY  = Velocity  of  JBOAT 

NFIRST  = Starting  index  in  radiated  signal  tables 

NLAST  = Last  index  in  radiated  signal  tables 
STEPWD  = Array  of  stepwidths  between  tables,  the  first 
one  being  the  velocity  for  the  first  table 
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NTABLE  = Number  of  tables  of  radiated  signals 

SOUND  = Output  from  subroutine  NOISY 

NOISY  has  three  formal  arguments: 

1.  Starting  address  of  table  of  radiated  signals 

2.  Starting  address  of  center  frequencies 

3.  Size  of  the  first  dimension  on  the  array  of  signals 


A 101 

UNCLASSIFIED 


UNCLASSIFIED 


ENTER 

RADSIG 


ERROR  STOP 


PRINT 

roono 

MESSAGE 

1 1 T rt - O 

I KINDJ-2 


ITYPE= 

l-JSNORK 


TYPE  = 4 


FREQUY=  FIISON 
VELOCY=  BSJBT 
> IT  =ITYPE 
NFIRST  = NSTART(IT) 
N=  I 


J NO 

I=NTABLE^ 

^ES 

NSIZE= 

116 


NLAST=  N START  (IT  + 1 )-l 
NTABLE=  NUM  (ITYPE) 
STEPWD(I)  =CUTOFF  (ITYPE) 
ITABr  NSTI(ITYPE) 


, STEPWD(I)  = 
CHAG  (1-2  + 1 TAB) 


/ CALL  NOISY 
TEMP=SOUND/ 


CALL  PR LOSS 
RAD  =TEMP  + PRLOSS 


RADSIG 


RETURN 

RADSIG 


FIG  A- 13  - RADSIG  FLOW  CHART 
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REFSIG 


PURPOSE:  This  function  subprogram  determines  the  active 

signal  received  by  IBOAT,  after  being  adjusted  for  JBOAT's  target 
strength  and  two-way  propagation  loss.  The  program  also  determines 
the  signal  received  at  JBOAT  due  to  IBOAT's  active  ping,  reduced 
by  the  one-way  propagation  loss. 

INPUT: 

ISONDX  = Sonar  index 

RI  = Depth  of  IBOAT  (negative) (yds) 

RJ  = Depth  of  JBOAT  (negative) (yds) 

RHOIJ  = Range  of  JBOAT  from  IBOAT  (yds) 

OUTPUT: 

Value  of  REFSIG 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS 

SOJ  = Source  level  adjusted  for  the  sending 

directivity  index  of  the  sonar 
Target  strength  for  the  active  sonar  is  assumed  to  be  +15 
dB  for  all  aspect  angles. 
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FIG.  A-14  - REFSIG  FLOW  CHART 
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SEANOI 


PURPOSE:  Function  SEANOI  determines  the  background  noise 

at  the  receiver  due  to  Sea  State.  This  is  accomplished  by  taking 
Sea  State  and  computing  the  Knudsen  value  for  the  spectrum  level 
of  the  sound  at  the  sonar  frequency.  This  Knudsen  value  is 
adjusted  for  the  depth  of  the  receiver,  the  directivity  index  of 
the  receiver,  and  the  bandwidth  of  the  receiver. 

INPUT: 

ISEAS  = Sea  State  number 

FIISON  = Geometric  or  center  frequency  of  sonar  (Hz) 

ATTENP  = Absorption  coefficient  for  the  propagation  of 
sound  through  sea  water  (dB/yd) 

RI  = Depth  of  IBOAT  (negative) (yds ) 

DIISON  = Receiving  directivity  index  of  sonar  (negative) 
(dB) 

DFIISN  = Receiving  bandwidth  of  sonar  (Hz) 

OUTPUT: 

The  value  of  SEANOI  is  the  only  output  of  the  routine 
INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

B(ISEAS+1)  - Array  of  Y-intercepts  for  Knudsen  curves 
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FIG.  A-15  — SEANOI  FLOW  CHART 
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SEFNOI 


PURPOSE:  This  function  subprogram  computes  the  self-noise 

seen  by  the  sonar  receiver  on  IBOAT.  This  self-noise  is  stored 
in  a table,  and  the  tabular  value  is  adjusted  for  the  receiving 
directivity  index  and  bandwidth  of  the  receiver. 

INPUT: 


IKINDI 

FIISON 

BSIBT 

DIISON 

DFIISN 

JSNORK 


Type  of  Boat  0 = Nuclear 

1 = Diesel-Electric 

Geometric  or  center  frequency  of  sonar  (I-’z) 
Velocity  of  IBOAT 

Receiving  directivity  index  of  sonar  (negative) 
(dB) 

Receiving  bandwidth  of  sonar  (Hz) 

Flag  for  snorkel  0 = No 

1 = Yes 


OUTPUT: 

The  value  of  SEFNOI  is  the  only  output  of  the  routine 
INTERNAL: 

Before  calling  NOISY  (TABLE  INTERPOLATION  ROUTINE),  the 
variables  of  common  block  SNTAB1  must  be  set. 


FREQUY 

= 

FIISON 

VELOCY 

= 

BSIBT 

NFIR3T 

= 

Index  of  start  in  noise  array 

NLA  ST 

= 

Index  of  last  in  noise  array 

STEPWD 

— 

1st  in  array  is  velocity  for  first 
are  stepwidths  in  velocity  between 

table,  rest 
tables 

NTABLE 

= 

Number  of  tables  = Z 

SOUND 

= 

Output  from  NOISY 

NOISY  also 
1. 
2. 

has  three  formal  arguments: 
Array  of  self-noise  values 
Array  of  operation  frequencies 

3.  Size  of  1ST  dimension  of  self-noise  values 
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FIG  A- 16  — SEFNOI  FLOW  CHART 
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STON 


PURPOSE:  The  purpose  of  this  function  subprogram  is  the  return 

the  s igna 1- to-noise  ratio  received  at  IBOAT  from  JBOAT,  when  IBOAT 
is  using  the  sonar  equipment  indexed  ISONDX. 

INPUT: 


IBOAT 

JBOAT 

ISONDX 

ISEAS 

ATTENP 

F 

DF 

DI 

ISNORK 


Index  of  boat  on  which  the  sonar  equipment  is 
located  (own  boat) 

Index  of  boat  which  is  the  target 

Index  indicating  the  type  of  sonar  equipment 

for  IBOAT 

Sea  State  number 

Array  of  absorption  coefficient  for  the 
propagation  of  sound  through  sea  water  (dB/yd) 
Geometric  or  center  frequency  of  sonar 
ISONDX  (Hz) 

Received  bandwidth  of  sonar  ISONDX  (Hz) 
Receiving  directivity  index  of  sonar  ISONDX 
(negative) (dB) 

Flag  for  snorkel 
0 = No 


1 = Yes 

R(3,JBQAT)=  Depth  of  JBOAT  (negative) (yds) 

R(3,IB0AT)=  Depth  of  IBOAT  (negative) (yds) 

BS(1,IB0AT)=  Velocity  of  IBOAT  (kts) 

RHO  = Actual  range  between  IBOAT  and  JBOAT 

IKIND  = Type  of  boat  0 = Nuclear 

1 = Diesel-Electric 

OUTPUT:  The  value  of  STON  is  the  only  output  of  the  program 

INTERNAL  VARIABLES  AND  CONSTANTS: 

The  value  +8  dB  has  been  added  to  the  S/N  ratio  calculated  from 
the  tables  to  correspond  with  the  SS/SSN  Preliminary  Cost  Effectiveness 
Analysis  noise  curves. 
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A. 4 TACTIC 

PURPOSE:  This  subroutine  evaluates  the  information  provided 

by  the  subsystems  of  each  boat,  and  makes  tactical  decisions 
based  on  the  tactical  doctrine  under  which  the  boat  is  oper- 
ating. 

As  the  simulation  currently  involves  only  two  boats, 
the  boat  making  the  decisions  is  termed  IBOAT  and  the  other 
boat  is  termed  JBOAT  for  purposes  of  describing  the  variables 
used  in  the  program. 

Certain  variables  are  generated  in  TACTIC  at  one  time 
step  for  use  by  TACTIC  at  a later  time  step.  As  such,  these 
variables  are  both  input  and  output  variables  of  tactics.  In 
addition,  there  are  certain  variables  which  are  either  pure  in- 
put or  pure  output.  For  this  routine,  all  variables  are  lumped 
into  one  general  category. 

VARIABLES  USED  IN  TACTIC: 


TIME 

= Current 

t ime  (min ) 

XMAX 

= Maximum 

X-coordinate 

value 

for 

barrier 

(yds) 

YMAX 

= Maximum 

Y-coordinate 

value 

for 

barrier 

(yds) 

ZMAX 

= Maximum 

Z-coordinate 

value 

for 

barrier 

(yds) 

XMIN 

= Minimum 

X-coordinate 

value 

for 

barrier 

(yds) 

YMIN 

= Minimum 

Y-coordinate 

value 

for 

barrier 

(yds) 

ZMIN 

= Minimum 

Z-coordinate 

value 

for 

barrier 

(yds) 
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R (3, IBOAT)  = Array  of  three  component  vectors 

(X,Y,Z)  specifying  current  posi- 
tion of  IBOAT  (yds) 

IKIND(IBOAT)  = Array  defining  the  type  of  IBOAT 

0 = Nuclear 

1 = Diesel-Electric 


ISNORK( IBOAT) 


TSS (IBOAT) 

TCS (IBOAT) 

DELT 

IPFLAG (IBOAT) 


= Array  specifying  the  current 
snorkeling  status  of  IBOAT 

0 = Not  Snorkeling 

1 = Snorkeling 

= Array  containing  the  next  time 

to  start  snorkeling  for  IBOAT  (min) 

= Array  containing  the  next  time 

to  cease  snorkeling  for  IBOAT  (min) 

= Time  step  for  the  simulation  (min) 

= Flag  indicating  that  IBOAT  has 
reacted  to  a ping  by  JBOAT 

0 = No 

1 = Yes 


IPING (JBOAT)  = Flag  indicating  that  IBOAT  has 

detected  JBOAT' s active  sonar  ping 

0 = No 

1 = Yes 

TOCCUR(I , IBOAT)  = Time  when  IBOAT  decides  to  enter 

tactical  situation  I (min) 


ITACT (IBOAT)  = Current  tactic  employed  by  IBOAT 

for  Boat  No.  1; 

1 = No  detection  has  occurred 

2 = Gather  data 

3 = Contact  Lost-Patrol 
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ITBR  (I, IBOAT) 


TDLY ( I , IBOAT) 
IDETEC( IBOAT) 

XPAT 

ZS( IBOAT) 
BS(I, IBOAT) 

CS(I, IBOAT) 
CT(I, IBOAT) 
TLSNK( IBOAT) 


4 = Contact  Lost -Chase 

5 = Launch  Weapon 

6 = Attempt  Ambush 

7 = Attempt  Tail-Chase 

8 = Evade 

9 = Close  and  Attack 
For  Boat  No.  2: 

1 = No  detection  has  occurred 

2 = Evade 

= Flag  indicating  whether  course 
and  speed  changes  for  tactical 
situation  I have  been  carried 
out  by  IBOAT  0 - yes 
1 = no 

= Time  delay  before  IBOAT  can 

enter  tactical  situation  I (min) 

= Array  flag  denoting  whether  IBOAT 
has  detected  anyone  0 = no 

1 = yes 

= Array  of  X-coordinates  for  the 
minimum  and  maximum  patrol  points 
for  Barrier  Boat  (small  number 
first,  yds) 

= Array  of  depths  at  which  IBOAT 
snorkels  (yds) 

= A seven  component  vector  which 
contains  the  motion  parameters  for 
each  boat  (see  SUBSYS  for  definition) 
= Snorkel  speed  of  IBOAT  during 
tactical  situation  I (kts) 

= Normal  speed  of  IBOAT  during 
tactical  situation  I (kts) 

= Array  time  IBOAT  last  completed 
snorkel  (min) 
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PHIRHO( IBOAT, JBOAT)  = True  bearing  of  JBOAT  relative  to 

IBOAT  (radians) 


DTHET (I BOAT, JBOAT, NSYS) 

MOTION (IBOAT) 

THET( IBOAT, JBOAT, NSYS) 

ITRK 

EDR ( I BOAT , JBOAT , I SNORK) 


IMP 

ESP (IBOAT, JBOAT) 
MFLGTK 

EWR( IBOAT, JBOAT) 
ZT (I , IBOAT) 


Array  of  errors  associated  with 
bearing  of  JBOAT  from  the  bearing 
subsystem  NSYS  on  IBOAT  (radians) 
Array  of  tag  for  type  of  motion 
undergone  by  IBOAT  (see  CONTOL) 

Array  of  bearing  of  JBOAT  from 
the  bearing  subsystem  NSYS  on 
IBOAT  (radians) 

Tag  indicating  number  of  points 
in  the  current  base  course  track 
IBOAT' s estimate  of  JBOAT1 s detection 
range  in  the  absence  of  any 
information  when  JBOAT  is  described 
by  I SNORK  (snorkeling  or 
non- snorke ling) 

Flag  indicating  whether  JBOAT  is 
snorkeling  (=2)  or  not  snorkeling 
(=1)  when  estimate  of  detection 
range  is  made 

IBOAT' s estimate  of  JBOAT 's  speed 
in  the  absence  of  all  information 
(kts) 

Array  of  flags  indicating  the 
existence  of  a satisfactory  (=1) 
or  unsatisfactory  (=0)  base 
course  track 

IBOAT 's  estimate  of  JBOAT 's  weapon 
range  (yds) 

Normal  depth  of  IBOAT  during 
tactical  situation  I (yds) 
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SNBC( I, I BOAT) 


PI2 

PI 

RD0T(2, IBOAT) 
CMAX(  IBOAT) 

ETS( IBOAT) 

ALFTRK 


XESTTK 


YESTTK 


AUSTIN,  TEXAS  78721 

= Array  of  minimum  and  maximum  values 
of  the  X-coorciinate  permitted  the 
Barrier  Boat  within  the  barrier 
(smaller  value  first,  yds) 

= Array  of  minimum  and  maximum  values 
of  the  Y-coordinate  permitted 
the  Barrier  Boat  within  the  barrier 
(smaller  value  first,  yds) 

= A fraction  for  IBOAT  representing 
the  battery  power  level  requiring 
snorkeling  during  tactical  situation 
I 

= 2n 

= Velocity  vector  for  IBOAT  (kts) 

= Array  of  maximum  velocities  for 
IBOAT  (kts) 

= IBOAT's  estimate  of  target's 
snorkel  frequency  (min) 

= Array  of  four  parameters  determined 
in  TRACK 

ALFTRK(2)  = estimated  X-velocity 
of  target  (kts) 

ALFTRK(4)  = estimated  Y-velocity 
of  target  (kts) 

= Estimate  of  current  X-coordinate 
of  target  on  his  base  course  in 
barrier  coordinate  system  (yds) 

= Estimate  of  current  Y-coordinate 
of  target  on  his  base  course  in 
barrier  coordinate  system  (yds) 
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ICLASS (IBOAT) 


MFLGFR 


DR ( IBOAT , JBOAT , NSYS ) 


WRMAX(IBOAT) 

RELB (IBOAT, JBOAT) 

I SEAS 

RHO ( I BOAT , J BOAT ) 

CMIN( I , IBOAT) 

ZLIM ( I , JBOAT) 

STN( IBOAT, JBOAT, NSYS) 

IFIR 

XESTFR 


YESTFR 


= Flag  indicating  that  IBOAT  has 
(=1)  or  has  not  (=0)  classified 
the  target 

= Array  of  tracking  flags  for  fire 
control  system  1 = accurate  track 

0 = inaccurate  track 

= Range  from  JBOAT  to  IBOAT  re- 
ported by  ranging  system  NSYS  on 
JBOAT  (yds) 

= Maximum  weapon  range  of  IBOAT  (yds) 

= Actual  relative  bearing  of  JBOAT 
with  respect  to  IBOAT  (radians) 

= Sea  State 

= Range  of  JBOAT  from  IBOAT  (yds) 

= Array  of  minimum  speeds  for  each 
boat  (kts) 

= A two  component  vector  for  JBOAT 
defining  operating  depths 
1=1  Minimum  depth  (yds) 

1=2  Maximum  depth  (yds) 

= Signal- to-noise  ratio  at  IBOAT 
with  system  NSYS  from  target 
JBOAT  (dB) 

= Number  of  points  in  fire  control 
track 

= Fire  control  estimate  of  X- 

coordinate  of  target  on  his  current 
course  in  the  barrier  coordinate 
system  (yds) 

= Fire  control  estimate  of  Y- 

coordinate  of  target  on  his  current 
course  in  the  barrier  coordinate 
system 
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EDR1 ( IBOAT, JBOAT) 

EWR1 ( I BOAT , JBOAT ) 

RFIRE(IBOAT) 
BLIND(I, IBOAT) 
ESCAPE 


PHI02( IBOAT) 
TTURN 

JEY 

BATLVI 

PHI2 


DTSAVE 

RSAVE2 

T SAVE 2 
RTEST 

IXG02 


= IBOAT' s estimate  of  JBOAT* s 

detection  radius  for  use  in  time 
to  shoot  determination  (yds) 

= IBOAT’s  estimate  of  JBOAT's  weapon 
range  for  use  in  time  to  shoot 
determination  (yds) 

= Maximum  range  at  which  IBOAT  will 
launch  his  weapon  (yds) 

= Array  for  IBOAT  of  blind  spot 
angles  (radians) 

= The  distance  from  the  boundary 
beyond  which  the  Penetrator  has 
definitely  escaped  the  Barrier 
Boat  (yds) 

= Basic  heading  of  IBOAT  (radians) 

= Time  for  Penetrator  to  execute 
next  turn  (min) 

= Array  of  random  number  keys 
= Current  battery  level  of  Boat  1 
= True  bearing  (with  error)  reported 
by  bearing  subsystem  used  in  tactical 
situation  2 (TS2) 

= Last  error  in  bearing  reported 
by  bearing  subsystem  (radians) 

= Last  range  reported  by  ranging 
subsystem  for  TS2  (yds) 

= Time  tag  for  RSAVE2 
= Estimate  of  Range  to  target  boat 
(yds) 

= Tag  used  to  indicate  whether 
Barrier  Boat  hits  boundary  (=2) 
in  tactical  situation  2 or  not 

(=D 
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VC4 

VC6 


ALF 
BETA 
TWA  IT 

PHI6 

RSAVE6 

TSAVE6 

PHI99 

RSAVE9 

TSAVE9 

PHI7 


= Velocity  for  Barrier  Boat  in  TS4 
(kts) 

= Velocity  for  Barrier  Boat  in 
proceeding  to  ambush  as  returned 
by  SSAMB  or  SSNAMB.  As  an  input 
to  these  routines  it  is  the 
minimum  speed  at  which  the  boats 
proceed  to  ambush  (kts) 

= Heading  Barrier  Boat  assumes 
when  he  reaches  ambush  position 
(radians) 

= Angle  through  which  Barrier  Boat 
must  turn  to  proceed  to  ambush 
(radians) 

= Length  of  time  Barrier  Boat  can 
wait  at  ambush  position  before 
target  enters  his  detection 
radius  (min) 

= Last  true  bearing  (with  error) 
from  the  bearing  subsystem  used  in 
TS  6 

= Last  range  from  the  ranging  sub- 
system used  in  TS  6 

= Time  tag  for  PHI6  and  RSAVE6 

= Last  true  bearing  of  target  from 
bearing  subsystem  (with  error) 
used  in  TS  9 

= Last  range  from  ranging  subsystem 
used  in  TS  9 

= Time  tag  for  PHI99  and  RSAVE9 

= Last  true  bearing  of  target  (with 
error)  from  bearing  subsystem 
used  in  TS  7 
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RSAVE7 

T SAVE 7 
R12 

XPAT3 

CSAVE6 

IXG06 

SNAV (IBOAT, JBOAT, I) 

PD(IBOAT, JBOAT, I) 
ISTART( IBOAT) 

NTIME( IBOAT) 

DETLM( IBOAT) 

NLOSS( IBOAT) 

IMCC( IBOAT) 


Last  range  reported  by  the  ranging 
subsystem  used  in  TS7 
Time  tag  for  PHI7  and  RSAVE9 
Current  range  to  target  used  in 
TS5  (yds) 

Array  of  patrol  limits  on  the 
X-coordinate  of  the  Barrier  Boat 
used  in  TS3  (smallest  value 
first,  yds) 

Estimate  of  target  velocity  from 
TRACK  used  in  TS6  (kts) 

Tag  used  to  indicate  whether 
Barrier  Boat  has  lost  contact  (=2) 
or  not  (=1)  with  target  in  TS6 
For  each  boat  array  of  average 
signal- to-noise  used  in  probability 
of  detection  portion  of  SUBSYS 
Probability  that  IBOAT  with 
system  I detects  JBOAT 
For  each  boat  flag  to  indicate 
whether  IBOAT  must  restart  detection 
process 

For  each  boat  number  of  points 
currently  in  S/N  table  for 
detection 

For  each  boat  lower  limit  for 
probability  of  detection 
For  each  boat  number  of  consecutive 
time  steps  for  which  detection 
has  been  lost 

Flag  indicating  that  IBOAT  has 
completed  a major  course  change 


A 
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MCC(IBOAT) 

JTS6 

IXG08 

JTACT 

BATLV2 

ISYSON(IBOAT, NSYS) 

Y1 

Y2 

PX2 

TX2 

THX2 

RSAVE3 

XLIM3 

PH14 


= Number  of  time  steps  after  IBOAT's 
course  change  that  the  other  boat 
realizes  a major  course  change 
has  occurred 

= Flag  to  indicate  possible  action 
upon  arrival  at  ambush  position 
by  Barrier  Boat 

= Tag  used  to  indicate  whether 
Barrier  Boat  hits  boundary  (=2) 
or  not  (=1)  in  TS8 
= Flag  indicating  that  simulation 
is  completed 

= Current  battery  level  of  Boat  2 
= Array  for  each  boat  and  subsystem 
indicating  status  of  subsystem 

0 = Off 

1 = On 

= Dummy  variable 
= Dummy  variable 

= Range  used  to  calculate  estimated 
range  to  target  (yds) 

= Time  used  to  calculate  estimated 
range  to  target  (min) 

= Last  bearing  from  the  bearing  sub- 
system in  TS2 

= Estimate  of  target's  range  when 
contact  was  lost  (yds) 

--  Used  to  readjust  patrol  path  as 
Barrier  Boat  moves  a distance 
XLIM3  from  the  line  he  estimates 
the  other  boat  is  taking  (yds) 

= Last  true  bearing  to  be  used  in 
TS4 
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RHOVEC 


RVEC 


WEC 


TTESTT 


CMAX6 


JAMB 


TSTOP 


PH19 

C9 

PHIR9 


THX7 


= Vector  of  target's  coordinates 
(X,Y,Z)  with  respect  to  own  ship 
= Vector  of  target's  coordinates 

(X,Y,Z)  in  barrier  coordinate  system 
= Estimate  of  target's  velocity 
vector  (V^jVyjV^)  in  knots 
= Estimated  time  from  present 
when  target  will  snorkel  (min) 

= Maximum  speed  of  Barrier  Boat 
proceeding  to  ambush  (kts) 

= Flag  indicating  ambush  can  (=0) 
or  cannot  (=1)  be  made 
= Time  from  beginning  of  intercept 
when  Barrier  Boat  reaches  ambush 
position  (min) 

= Estimate  of  current  X-coordinate 
of  target  used  in  TS9 
= Estimate  of  current  Y-coordinate 
of  target  used  in  TS9 
= Estimate  of  current  velocity  of 
target  used  in  TS9 
= Heading  of  other  boat 
= Normal  speed  of  Barrier  Boat 

during  close  and  attack,  TS9  (kts) 

= Last  true  bearing  of  target  from 
bearing  subsystem  (with  error) 
for  use  in  TS9 

= Probability  of  Kill  in  TS5  (only 
zero  value  is  contained  in  this 
program  as  CHOICE  IIB  calculates 
PK) 

= Last  bearing  from  the  bearing 
subsystem  in  TS7 
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CTEST7  = (Current  estimate  of  target's 

velocity  along  his  base  course  (kts 
MTK  = Redetection  flag 

Threshold  for  detection  with  AN/BQR-2  (DIMUS)  is  a 
signal- to-noise  of  -33.3  dB 
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BEGIN  TACT  I CAL 
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FOR  BOAT 
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s'  HAS 
"COURSE  BEEN 
INITIALIZED 
\ FOR  TS5^ 


'"RI2>WRMAX(D 

(TARGET  IN 
^RANGEV* 


PRINT  ERROR 
MESSAGE 

PK  O 

Vjja ct^_j  / 

\return  / 


INITIALIZE 
VARIABLES 
FOR  THE 
TORPEDO 
PROGRAM 


CALL 
CHO I CE 
(TORPEDO) 


JT  ACT  = 1 
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SAVE  BEARING  INFO/ 

SAVE  PUFFS 

FROM  BQR-2 ( ) IF 
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ATCS 


PURPOSE:  This  function  subprogram  computes  the  next  time  to 

cease  snorkeling  for  each  boat. 

INPUT: 

T = Time 

I = Number  of  boat  for  which  ATCS  is  to  be  computed 

OUTPUT: 

Value  of  the  function  ATCS  = Time  to  cease  snorkeling 
INTERNAL  CONSTANTS: 

Boat  1 stops  snorkeling  again  in  20  min 
Boat  2 stops  snorkeling  again  in  35  min 
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ATSS 


PURPOSE:  This  function  subprogram  computes  the  next  time  to 

start  snorkeling  for  each  boat. 


INPUT: 


OUTPUT: 


T = Time 

I = Number  of  boat  for  which  ATSS  is  to  be  computed 


Value  of  the  function  ATSS  = Time  to  start  snorkeling 


INTERNAL  CONSTANTS: 


Boat  1 starts  snorkeling  again  in  20  min. 
Boat  2 starts  snorkeling  again  in  30  min. 
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BATLl 

■ 

PURPOSE:  This  routine  is  designed  to  compute  the  battery  level  at 
each  time  step  for  the  Barrier  Boat.  The  routine  has  a branch 
determined  by  whether  or  not  the  Barrier  Boat  is  snorkeling. 

Currently  the  routine  is  a dummy;  the  battery  level  is  set  to 
1.0  each  time. 
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BATL2 

PURPOSE:  This  routine  computes  the  battery  level  for  the 

Penetrator  Boat  at  each  time  step,  branching  on  whether  that  boat 
is  or  is  not  snorkeling.  The  routine  is  currently  a dummy;  the 
battery  level  is  always  set  to  1.0  at  each  time  step. 
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INPUT: 

X = Current  X-coordinate  of  Barrier  Boat 
Y = Current  Y-coordinate  of  Barrier  Boat 
XL  = Lower  limit  for  X-coordinate  inside  barrier 
XH  = Upper  limit  for  X-coordinate  inside  barrier 
YL  = Lower  limit  for  Y-coordinate  inside  barrier 


YH  = Upper  limit  for  Y-coordinate  inside  barrier 
PKIT  = True  bearing  of  target  from  Barrier  Boat 
PHI  = Current  heading  of  Barrier  Boat 


OUTPUT: 

PHI  = New  course  for  Barrier  Boat 
X = New  X-position  of  Barrier  Boat 

Y = New  Y-position  of  Barrier  Boat 
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PURPOSE:  This  routine  handles  the  case  in  which  the  Barrier 

Boat  strikes  the  boundary  of  the  barrier  and  runs  parallel  to 

it. 


INPUT: 

X 

Y 

XL 

XH 

YL 

YH 

PH  IT 

PHI 

OUTPUT: 

PHI 

X 

Y 


Current  X-coordinate  of  Barrier  Boat 
Current  Y-coordinate  of  Barrier  Boat 
Lower  limit  for  X-coordinate  of  boundary 
Upper  limit  for  X-coordinate  of  boundary 
Lower  limit  for  Y-coordinate  of  boundary 
Upper  limit  for  Y-coordinate  of  boundary 
True  bearing  of  target  from  Barrier  Boat 
Current  heading  of  Barrier  Boat 

New  course  for  Barrier  Boat 
New  X-position  of  Barrier  Boat 
New  Y-position  of  Barrier  Boat 
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BTURN 


PURPOSE:  This  subroutine  is  called  when  the  Barrier  Boat  is 
between  its  maneuvering  limits  and  the  sides  of  the  barrier.  It 
is  only  called  from  TS1,  and  returns  the  heading  which  will  move 
die  boat  away  from  the  nearest  barrier  boundary. 

INPUT: 


X 
Y 

XI 
X2 
Y1 
Y2 

OLTP  IT: 

PHI 


Current  X-coordinate  of  Barrier  Boat 
Current  Y-coordinate  of  Barrier  Boat 
Minimum  value  of  X-coordinate  of  maneuvering  limit 
Maximum  value  of  X-coordinate  of  maneuvering  limit 
Minimum  value  of  Y-coordinate  of  maneuvering  limit 
Maximum  value  of  Y-coordinate  of  maneuvering  limit 

New  heading  for  Barrier  Boat 
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CCHG 

PURPOSE:  This  function  computes  the  course  required  for  the 

Barrier  Boat  so  that  a contact  at  the  apparent  true  bearing  (A+E) 
hears  +D  relative.  A is  the  true  bearing  and  E is  the  true  bearing 
error . 

A = True  bearing  of  target  (without  error) 

E = Error  in  bearing  reported  by  bearing  subsystem 
D = Angle  target  is  to  make  from  bow  on  new  course  of 
Barrier  Boat 

OUTPUT: 

Value  of  the  function  CCHG  = New  heading  of  Barrier  Boat 
INTERNAL  VARIABLES: 

T = True  bearing  returned  by  the  bearing  subsystem  (with  error) 
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CNC 


PURPOSE:  This  function  computes  the  new  heading  for  the 

Penetrator  by  adding  a randomly  chosen  angle  within  given  limits 
:.o  that  boat's  basic  heading. 

INPUT: 

PHI  = Base  course  heading  of  ship  (radians) 

key  = Array  of  random  numbers 

OUTPUT: 

Value  of  the  function  CNC  = New  heading  of  boat 
INTERNAL  CONSTANTS: 

The  current  limits  on  the  maneuver  are  that  it  must  be 
between  +20°  of  the  base  course. 
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PURPOSE:  This  function  subprogram  computes  the  next  time  to 

turn  for  the  Penetrator  on  his  evasive  course. 

INPUT: 

TIME  = Time  (min) 

OUTPUT: 

Value  of  the  function  CTT  = Time  to  turn 
INTERNAL  CONSTANTS: 

Penetrator  maneuver  frequency  = 5 minutes 
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FIRCON 

PURPOSE:  This  program  attempts  to  determine  a fire  control 

quality  track  of  the  target.  This  is  accomplished  by  fitting  a 
Least  squares  line  to  the  X and  Y motion  of  the  target,  independ- 
ently. This  routine  assumes  constant  velocity  motion  by  the  target. 

INPUT : 

RHO  = Range  from  the  ranging  subsystem  (yds) 

THETA  = True  bearing  of  other  boat  from  the  bearing 
subsystem  (radians) 

XO  = X-coordinate  of  own  ship's  estimate  of  own 

position  in  an  Earth-fixed  coordinate  system  (yds) 
YO  = Y-coordinate  of  own  ship's  estimate  of  own  position 

in  an  Earth- fixed  coordinate  system  (yds) 

= Number  of  points  in  the  current  track 
= Time  since  the  beginning  of  the  simulation  (min) 


TP INTI 
TIME 
OUTPUT: 

ALPHA 


SIGMA 
I FLAG 
XEST 
YE  ST 
TLA  ST 


Array  of  four  parameters  to  characterize  motion 
of  target 


X = a i + a 2^ 


Y = a.^  + a^t 

Array  of  standard  deviations  of  X(ct^)  and  Y(<J2)  ' 
values  ibout  the  least  squares  lines  (yds) 

Array  of  two  tracking  flags  to  indicate  whether 
track  is  of  sufficient  quality  to  launch  a weapon 
Current  estimate  of  target's  X-position  in  Earth- 
fixed  coordinate  system 

Current  estimate  of  target's  Y-position  in  Earth- 

fixed  coordinate  system 

Time  when  track  was  last  updated  (min) 
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FIG.  A- 26  - CTT  FLOW  CHART 
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INTERNAL  VARIABLES  AND  INTERNAL  CONSTANTS: 

IPINT2  = Internal  flag  used  to  determine  if  track  is  to  be 
updated 

XKNYPM  = Conversion  factor  from  knots  to  yards  per  min 

K = Dimensioned  array  of  variables  against  which 

SIGMA  is  compared  to  determine  if  tracking  flag 
IFLAG  should  be  set 


XUMXSQ 

= 

*i2 

SUMXT 

= 

E 

Vi 

SUMX 

= 

E 

Xi 

SUMT 

= 

E 

Ci 

SUMTSQ 

= 

E 

Ci2 

SUMYT 

= 

E 

Vi 

SUMY 

= 

E 

Yi 

SUMYSQ 

= 

E 

Yi2 

X.^  = XQ  + p sinG 
= XQ  + p cos0 
t±  = TIME  - T START 


N 

(All 

B 

CXCA 

ALF 

BETA 

T2 


= Number  of  points  in  current  track 
sums  are  taken  from  i = 1 to  i = N.) 

= Range  between  boats 

= Velocity  necessary  to  intercept  target  boat  on  a 
]_  course  (kts) 

= Interior  angle  of  close  and  attack  triangle  at 
target  position 

= Interior  angle  of  close  and  attack  triangle  at 
own  ship  position 

= Intermediate  variable  used  to  compute  final 

heading  which  assures  that  the  Barrier  Boat  always 
leads  the  target. 
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ENTER 

FIRCON 


INITIALIZE  INTERNAL 
VARIABLES  TO  BEGIN 
A NEW  TRACK 


IP  I NT  1 


RETURN 


"I  RACK  CONTAIN 
AT  LEAST  THREE 
\ POINTS  ^ 


SUM  XT  + a?2  SUM  TSO  + 2a,a, 
SUM  YT  + a 2 SUM  TSQ  +2a,a, 


ALL  errs 
BEEN  TESTED 


SET  IFLAG 


FIG.  A-27  - FIRCON  FLOW  CHART 
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PURPOSE:  This  routine  computes  the  heading  for  the  Barrier 

Boat  which  permits  him  to  intercept  the  enemy  in  the  close 
and  attack  situation  (TS9).  The  heading  is  also  chosen  subject 
to  the  constraint  that  PUFFS  be  able  to  effectively  bear  on  the 
target . 

INPUT: 


(XI, Yl)  = 


(X2,Y2)  = 


u2 

PHI2 
PH  IT 


Current  position  of  Barrier  Boat  in  the  barrier 
coordinate  system 

Speed  of  Barrier  Boat  (returned  as  speed  with 
which  boat  proceeds  to  intercept)  in  knots 
Current  position  of  target  in  the  barrier 
coordinate  system 
Speed  of  target  boat  (kts) 

Heading  of  target 

Last  true  bearing  of  target  reported  by 
bearing  subsystem 


OUTPUT: 


PHICA  = Heading  of  Barrier  Boat  to  effect  intercept 
(radians) 

Cl  = Speed  to  intercept  (kts) 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

X = Relative  X-separation  of  boats 

Y = Relative  Y- separation  of  boats 

A = The  dot  product  of  the  relative  position 

vector  and  the  target  velocity  vector 
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PURPOSE:  This  routine  is  called  when  it  is  time  for  any  boat 

to  snorkel.  The  routine  adjusts  the  depth  and  initiates  the 
snorkeling.  The  routine  calls  ATSS  to  set  the  next  snorkel 
time. 


INPUT: 

T 

Z2 

C2 

12 


Time 

Depth  at  which  boat  snorkels  (yds) 
Speed  at  which  boat  snorkels  (kts) 
Index  number  of  boat 


OUTPUT: 

Z1 

Cl 

II 


T3 


Current  Depth  of  boat 
Current  speed  of  boat 

Flag  to  indicate  whether  or  not  boat  is 
snorkeling  0 =•  No 

J.  = Yes 

Next  time  to  start  snorkeling 
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PURPOSE:  This  routine  is  called  when  it  is  time  for  a boat 

to  cease  snorkeling.  It  adjusts  the  depth  and  speed  to  those 
read  in  as  appropriate  to  the  boat,  and  sets  the  next  time  to 
cease  snorkeling  by  calling  ATCS. 


INPUT: 

T 

= 

Time 

Z2 

= 

Depth  for  boat  after  snorkeling  (yds) 

C2 

= 

Speed  for  boat  after  snorkeling  (kts) 

12 

= 

Index  number  of  boat 

OUTPUT: 

Z1 

= 

Current  depth  of  boat 

Cl 

= 

Current  speed  of  boat 

11 

Flag  to  indicate  whether  boat  is  snorkeling 

0 = No 

1 = Yes 

Tl 

= 

Next  time  for  boat  to  cease  snorkeling 

T2 

= 

Time  when  boat  last  completed  snorkeling 
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II  =0 
Zl  =Z2 
C I =C2 
T2=T 


&l=ATCS(T,I2)y 


RETURN 


FIG.  A-30  -SN0RK2  FLOW  CHART 
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PURPOSE:  This  subroutine  determines  if  the  Barrier  Boat 

can  ambush  an  enemy  Diesel-Electric  boat  based  on  the  information 
currently  available  to  the  Barrier  Boat.  If  the  ambush  can  be 
consumated,  the  routine  returns  the  parameters  to  guide  the 
Barrier  Boat  to  the  ambush  position. 


INPUT: 

R 

V 

T 

VC 


DWR 

TEST 

CMAX 

RHO 


XMAN 

YMAN 


OUTPUT: 

ALPH 

BETA 


Current  estimate  of  the  target  position  vector 
in  barrier  coordinates 

Current  estimate  of  the  target  velocity  vector 
in  barrier  coordinates  (kts) 

Time 

Minimum  speed  at  which  Barrier  Boat  will 
attempt  to  ambush  (returned  as  speed  at  which 
Barrier  Boat  will  actually  proceed  to  ambush) 
in  kts 

Barrier  Boat  detection  range  (yds) 

Estimated  time  from  present  when  target  will 
snorkel  (min) 

Maximum  speed  of  Barrier  Boat  (kts) 

Current  estimate  of  the  relative  position  vector 
from  Barrier  Boat  to  target  in  barrier 
coordinates  (yds) 

Minimum  and  maximum  values  of  X-coordinate  inside 
which  Barrier  Boat  must  maneuver 
Minimum  and  maximum  values  of  Y-coordinate 
inside  which  Barrier  Boat  must  maneuver 

Heading  to  be  assumed  upon  reaching  ambush  pivot 
point  (radians) 

Angle  between  present  heading  and  ambush  heading 
for  Barrier  Boat  (radians) 
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TSTOP 
TWA  IT 


JFLGSS 


Time  from  beginning  of  simulation  when  Barrier  Boat 
will  reach  ambush  pivot  point  (min) 

Time  from  beginning  of  simulation  when  target 
will  enter  Barrier  Boat's  detection  range  EWR 
after  Barrier  Boat  reaches  ambush  pivot  point  (min) 
Flag  for  ambush  0 = can  ambush 

1 = cannot  ambush 


INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

XKNYPM  = Conversion  factor  from  knots  to  yards  per  minute 

DELT  = Time  step  for  incrementing  ambush  calculation 

D2  = Length  of  target  leg  of  triangle  for  time  TEST 

D1  = Length  Barrier  Boat  leg  of  triangle  for  time  TEST 

V2  = Estimated  speed  of  target  boat  (yards  per  minute) 

(XE,YE)  = Vector  of  estimated  target  position  at  time  TEST 

from  present 

DTEST  = Length  of  the  sum  of  the  vector  RHO  and  a vector 

of  length  D2  with  direction  of  V 
RHOMAG  = Magnitude  of  RHO 

Cl  = X-component  of  vector  along  V with  magnitude  D2 

C2  = Y-component  of  vector  along  V with  magnitude  D2 
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PURPOSE:  This  subroutine  determines  if  the  Barrier  Boat  can 

ambush  an  enemy  Nuclear  boat  based  on  the  information  currently 
available  to  the  Barrier  Boat.  If  the  ambush  can  be  consumated, 
the  routine,  returns  the  parameters  to  guide  the  Barrier  Boat 
to  the  ambush  position. 


INPUT: 

R 

V 

T 

VC 


DWR 

CMAX 

RHO 


XMAN 

YMAN 


OUTPUT: 

ALPH 

BETA 

TSTOP 


Current  estimate  of  the  target  position  vector 
in  barrier  coordinates  (yds) 

Current  estimate  of  the  target  velocity  vector 
in  barrier  coordinates  (kts) 

Time 

Minimum  speed  at  which  Barrier  Boat  will  attempt 
to  ambush  (returned  as  speed  at  which  Barrier 
Boat  will  actually  proceed  to  ambush)  (kts) 

Barrier  Boat  detection  range  (yds) 

Maximum  speed  of  Barrier  Boat  (kts) 

Current  estimate  of  the  relative  position  vector 
from  Barrier  Boat  to  target  in  barrier  coordinates 
(yds) 

Minimum  and  maximum  values  of  X-coordinate  inside 
which  Barrier  Boat  must  maneuver 

Minimum  and  maximum  values  of  Y- coordinate  inside 
which  Barrier  Boat  must  maneuver 

Heading  to  be  assumed  upon  reaching  ambush  pivot 
point  (radians) 

Angle  between  present  heading  and  ambush 
heading  for  Barrier  Boat  (radians) 

Time  from  beginning  of  simulation  when  Barrier  Boat 
will  reach  ambush  pivot  point  (min) 
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TWAIT  = Time  from  beginning  of  simulation  when  target 
will  enter  Barrier  Boat's  detection  range  EWR 
after  Barrier  Boat  reaches  ambush  pivot  point  (min) 
JFLGSN  = Flag  for  ambush  0 = can  ambush 

1 = cannot  ambush 

INTERNAL  VARIABLES  AND  ARBITRARY  CONSTANTS: 

XKNYPM  = Conversion  factor  from  knots  to  yards  per  minute 
DELT  = Time  step  for  incrementing  ambush  calculations 

D2  = Length  of  target  leg  of  triangle  for  time  TEST 

Dl  = Length  of  Barrier  Boat  leg  of  triangle  for 

time  TEST 

V2  = Estimated  speed  of  target  boat  (yds  per  min) 

(XE,YE)  = Vector  of  estimated  target  position  at  time  TEST 

from  present 

D = Length  of  the  sum  of  the  vector  RHO  and  a vector 

of  length  D2  with  direction  of  V 
RHOMAG  = Magnitude  of  RHO 

Cl  = X-component  of  vector  along  V with  magnitude  unity 

C2  = Y-component  of  vector  along  V with  magnitude  unity 

ALPH  = Internal  angle  in  triangle  opposite  Barrier  Boat 
leg 

C = Horizontal  range  between  boats 

TO  = Time  from  present  for  which  boat  motions  are 

extrapolated 

D2PRM  = Initial  estimate  for  length  of  target  leg  of 
triangle  (yds) 
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TRACK 


PURPOSE:  This  program  attempts  to  determine  a base  course  for 

a maneuvering  target.  This  is  accomplished  by  fitting  a least 
squares  line  to  the  X and  Y motion  of  the  target,  independently. 
The  routine  assumes  constant  velocity  motion  in  both  X and  Y 
by  the  target. 


INPUT: 

RHO 

THETA 

XO 

YO 

IP  INTI 
TIME 


Range  from  the  ranging  subsystem  (yds) 

True  bearing  of  other  boat  from  the  bearing 
subsystem  (radians) 

X-coordinate  of  own  ship's  estimate  of  own 
position  in  an  Earth-fixed  coordinate  system  (yds) 
Y-coordinate  of  own  ship's  estimate  of  own 
position  in  an  Earth-fixed  coordinate  system  (yds) 
Number  of  points  in  the  current  track 
Time  since  the  beginning  of  the  simulation  (min) 


OUTPUT: 

ALPHA 


SIGMA 

IFLAG 


XEST 
YE  ST 


Array  of  four  parameters  to  characterize  motion  of 
target 

X = a i + a 2^ 

Y = ctg  + a^t 

Array  of  standard  deviations  of  X(c^)  and  Y(o2) 
values  about  the  least  squares  lines  (yds) 

Array  of  two  tracking  flags  to  indicate  whether 
track  is  of  sufficiently  good  quality  to  launch 
a weapon 

Current  estimate  of  target's  X-position  in  Earth- 
fixed  coordinate  system 

Current  estimate  of  target's  Y-position  in  Earth- 
fixed  coordinate  system 
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TLAST  = Time  when  track  was  last  updated  (rain) 


INTERNAL  VARIABLES  AND  INTERNAL  CONSTANTS: 

IPINT2  = Internal  flag  used  to  determine  if  track  is  to  be 
updated 

XKNYPM  = Conversion  factor  from  knots  to  yards  per  min 
K = Dimensioned  array  of  variables  against  which 

SIGMA  is  compared  to  determine  if  tracking  flag 
IFLAG  should  be  set 

TSTART  = Time  initial  data  on  current  track  was  accrued 
SUMXSQ  = £ X? 

SUMXT  = E Xiti 


SUMX  = V xi 

SUMT  = E t± 

SUMTSQ  = r t±2 


X^  = XQ  + p sinG 
= XQ  + p cos9 


SUMYT  = I Yiti 


t±  = TIME  - TSTART 


SUMY  = V Y^ 

SUMYSQ  = r Yt2 


N = Number  of  points  in  current  track 

(All  sums  are  taken  from  i = 1 to  i = N) 

TTRACK  = Length  of  time  for  which  data  must  be  accumulated 
before  a base  course  for  the  target  can  be 
predicted  (min) 
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Washington,  D.C.  20360 


In  the  front  portion  of  this  Technical  Memorandum,  the 
general  approach  to  the  complete  digital  computer  simulation  of  Anti- 
Submarine  Warfare  (ASW)  Combat  is  formulated  and  discussed.  The 
discussion  centers  upon  the  three  basic  components  of  combat — command 
and  control  functions,  machines  and  the  environment — and  identifies 
the  general  parameters  and  dynamic  variables  required  to  simulate 
these  components.  (While  the  primary  concern  is  ASW  combat,  the 
general  princit  es  are  applicable  to  any  type  of  combat^  The  second! 
portion  of  this  memorandum  presents  a detailed  discription  of  the  sim- 
ulation of  combat  between  two  submarines.  In  this  simulation  one  sub4 
marine  is  deployed  on  a Forward  Area  Patrol  Mission.  The  mission  of  the 
second  submarine  is  to  transit  the  patrol  area  of  the  first  submarine [ 
This  digital  computer  simulation  contains  many  features  discussed  in  jhe 
general  approach  and  demonatrates  the  feasibilty  of  constructing  a 
general  ASW  Combat  Simulation  for  a digital  computer 
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